MONTHLY NOTICES 


OF THE 


ROYAL ASTRONOMICAL SOCIETY 


Vol. 114 No. 4 1954 


Published and Sold by the 


ROYAL ASTRONOMICAL SOCIETY 
BURLINGTON HOUSE 
LONDON, W.1 


Price Thirteen Shillings and Sixpence 





NOTICE TO AUTHORS 


1. Communications.—Papers must be communicated to the Society by a Fellow. They 
should be accompanied by a summary at the beginning of the paper conveying briefly the 
content of the paper, and drawing attention to important new information and to the main 
conclusions. The summary should be intelligible in itself, without reference to the paper, 
to a reader with some knowledge of the subject; it should not normally exceed 200 words 
in length. Authors are requested to submit MSS. in duplicate. These should be 
typed using double spacing and leaving a margin of not less than one inch on the 
left-hand side. Corrections to the MSS. should be made in the text and not in the 
margin. Unless a paper reaches the Secretaries more than seven days before a Council 
meeting it will not normally be considered at that meeting. By Council decision, MSS. of 
accepted papers are retained by the Society for one year after publication; unless their 
return is then requested by the author, they are destroyed. 


2. Presentation.—Authors are allowed considerable latitude, but they are requested to 
follow the general style and arrangement of Monthly Notices. References to literature 
should be given in the standard form, including a date, for printing either as footnotes or in 
a numbered list at the end of the paper. Each reference should give the name and 
initials of the author cited, irrespective of the occurrence of the name in the text (some 
latitude being permissible, however, in the case of an author referring to his own work). 
The following examples indicate the style of reference appropriate for a paper and a book, 
respectively :— 

A. Corlin, Zs. f. Astrophys., 15, 239, 1938. 

H. Jeffreys, Theory of Probability, 2nd edn., section 5.45, p. 258, Oxford, 1948. 


3. Notation.—For technical astronomical terms, authors should conform closely to the 
recommendations of Commission 3 of the International Astronomical Union (77ans. 
1.A.U., Vol. VI, p. 345, 1938). Council has decided to adopt the I.A.U. 3-letter abbrevi- 
ations for constellations where contraction is desirable (Vol. IV, p. 221, 1932). In general 
matters, authors should follow the recommendations in Symbols, Signs and Abbreviations 
(London: Royal Society, 1951) except where these conflict with I.A.U. practice. 


4. Diagrams.—These should be designed to appear upright on the page, drawn 
dine twice the size required in print and prepared for direct photographic 
reproduction except for the lettering, which should be inserted in pencil. 
Legends should be given in the manuscript indicating where in the text the 
figure should appear. Blocks are retained by the Society for 10 years; unless the author 
requires them before the end of this period they are then destroyed. 


5. Tables.—These should be arranged so that they can be printed upright on 
the page. 


6. Proofs.—Costs of alterations exceeding 5 per cent of composition must be borne by 
the author. Fellows are warned that such costs have risen sharply in recent years, and it 
is in their own and the Society’s interests to seek the maximum conciseness and simplifi- 
cation of symbols and equations consistent with clarity. 


7. Revised Manuscripts ——When papers are submitted in revised form it is especially 
requested that they be accompanied by the original MS. 


Reading of Papers at Meetings 
8. When submitting papers authors are requested to indicate whether they will be 


willing and able to read the paper at the next or some subsequent meeting, and approxi- 
mately how long they would like to be allotted for speaking. 


g. Postcards giving the programme of each meeting are issued some days before the 
meeting concerned. Fellows wishing to receive such cards whether Ordinary 
Meetings or for the Geophysical Discussions or both should notify the Assistant Secretary. 
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MEETING OF 1954 MARCH 12 


Dr J. Jackson, President, in the Chair 
The election by the Council of the following Fellows was duly confirmed :— 


Arthur Thomas Goodchild, 3. Murley Road, Winton, Bournemouth, Hamp- 
shire (proposed by G. Merton); 

Francis James Hooper, Nineteen Steps, Orchard Hill, Bideford, Devon- 
shire (proposed by F. Molesworth); 

*Harry Hughes, Ph.D., 27 Yewlands Avenue, Leyland, Preston, Lancashire 

(proposed by F. Birch); 

Charles Lloyd, 11 Benson Avenue, Goldthorn Park, Wolverhampton, 
Staffordshire (proposed by M. Johnson); and 

David Kenneth Northrup, F.R.S.A., 62 Ballater Road, London, S.W.2 
(proposed by A. W. Lane Hall). 


The election by the Council of the following Junior Members was duly 
confirmed :— 
+Thomas William Morris, 21D Berryland Road, Moreton Wirral, Cheshire 
(proposed by J. Kershaw); and 
Raymond Henry Preston, 4 South Road, West Kirby, Cheshire (proposed 
by J. Kershaw). 


Sixty-six presents were announced as having been received since the last 
Meeting, including :— 
Max Waldmeier, Radiowellen aus dem Weltram (presented by the author); 


and 
Wolfgang Gleissberg, Die Haufigkeit der Sonnenflecken (presented by the 
author). 


MEETING OF 1954 APRIL 9 
Dr J. Jackson, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 


George Edward William Gosnell, 75 Woodmansterne Road, Coulsdon, 
Surrey (proposed by F. M. Holborn); 


* Transferred from Junior Membership. 
+ Election later declared to be null and void. 
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Guillermo Haro, ‘Tonanzintla Observatory, Ap. Postal 216, Puebla, Mexico 
(proposed by G. Miinch); 
*Garry Ross Nankivell, 45 Gayhurst Road, Dallington, N.E.1, Christchurch, 
New Zealand (proposed by A. F. A. L. Jones); 
Sibte Nabi Naqvi, Department of Meteorology, Government of Pakistan, 
(proposed by A. Hunter); 
Domingo Taboada, Av. Reforma, 2111, Puebla, Mexico (proposed by G. 
Miinch); 
John Graham Whitaker, B.Sc., 26 Woodfield Lane, Ashtead, Surrey 
(proposed by T. G. Cowling); and 
Wilhelm Levi Whittle, Parkstone, 706 Doncaster Road, Foulby, Wakefield, 
Yorkshire (proposed by G. P. Meredith). 


The President announced that the Council had elected the following 
Associates of the Society :-— 


Jean Dufay, Professor of Astronomy and Director of the Observatory, 
Université de Lyon, St Genis-Laval, Rhone, France; 

Max Waldmeier, Professor of Astronomy and Director of the Eidgenossische 
Sternwarte, Ziirich, Switzerland; and 

Herman Zanstra, Professor of Astronomy in the University of Amsterdam 


and Director of the Astronomical Institute of the University, Amster- 
dam, Holland. 


One hundred and one presents were announced as having been received 
since the last Meeting, including :-— 


American Association for the Advancement of Science, Astronomical 
photoelectric photometry (presented by the American Association for 
the Advancement of Science); 

Astronomische Gesellschaft, Portrdtgallerie, 1931 (presented by F. Gonda- 
latsch) ; 

Hamburg Observatory, Bergedorfer Spektraldurchmusterung, Band V; 
and Zweiter Katalog der Astronomischen Gesellschaft, Bande 5-9 (pre- 
sented by the Hamburg Observatory); 

International Geodetic and Geophysical Union, La deuxiéme opération 
internationale des longitudes, Octobre-Novembre, 1933 (presented by the 
International Geodetic and Geophysical Union); and 

Bernhard Sticker, Fiinfstelige Tafel der trigonometrischen Funktionen, 
Ausgabe A (presented by the author). 


* Transferred from Junior Membership. 














WHOLE-SKY STATISTICS OF CELESTIAL OBJECTS 
C. W. Allen 
(Received 1954 February 11) 


Summary 


The advantages of statistical representations based on only those celestial 
objects with greatest apparent brightness have been enumerated. ‘The 
system of observations and results that has been devised can be applied 
to almost any recognizable category of celestial object. The observable 
parameters are m,(10), g, M, oy and 5, representing respectively magnitude of 
the tenth object, change of number with magnitude, absolute magnitude and 
its dispersion, and mean galactic latitude. The equations and numerical data 
are given for conversion of these (in most cases) into the descriptive quantities 
Mb), «, Do and y, representing respectively the absolute magnitude and its 
dispersion, the space density and scale height perpendicular to the galactic 
plane. Some alternatives to My and Dy, are given which are not so sensitive 
to the evaluation of o. The mean absorption coefficient and optical 
thickness of the galaxy must be known. 

The analysis has been exemplified by application to Wolf-Rayet stars, 
@ stars, flare stars, all stars, white dwarfs, globular clusters, open clusters, 
extra-galactic nebulae, novae and A stars. 





1. Introduction.—The first step towards a systematic study of the many 
objects found in the sky is to segregate them into various classes, groups, or cate- 
gories. One can then assume a certain uniformity within each class and proceed 
to investigate its intrinsic and statistical characteristics. Observations made 
on individual objects may serve the dual purpose of revealing intrinsic character- 
istics and also assigning each object to its appropriate class. Once sufficient 
objects are classified into a particular category one may proceed with statistical 
studies from which to determine the number density of the objects, their number 
distribution and their motions. (‘The statistics in the present paper are con- 
cerned with numbers only—not motions.) 

The ideal treatment for any particular category would aim at a complete 
analysis of number density as a function of distance. ‘The data would be obtained 
by counts of the object per square degree and per magnitude range for a large 
part of the sky. Such an analysis may be possible for a few well-defined cate- 
gories, but more usually there is too much difficulty in recognizing the objects to 
to be counted. Hence the present study is based entirely on the nearer and 
brighter objects of each category. The objects actually used are the ten or twenty 
brightest objects found anywhere in the sky. 

The advantages of basing a statistical analysis on this selection are as follows: 

(a) There is no arbitrariness in the choice of object since those of lowest 
apparent magnitude are used exclusively. Results are therefore as general as 
possible. 

26* 
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(6) Those objects used will most probably be the best-recognized members 
of their class, and the physical measurements are likely to be the most accurate 
and complete. 

(c) Working with small numbers involves very little labour and justifies 
obtaining a result even if high accuracy cannot be expected. The small list of 
objects may be constantly revised, and the effects of changing values can be 
readily assessed. 

(d) The procedure can be applied consistently for all categories of object. 

(e) Since the objects used are randomly taken from any part of the sky it is 
legitimate to use a mean value for interstellar absorption. 

The whole-sky analysis of celestial objects to be described aims at exploiting 
these advantages. 

2. Parameters describing the true distribution.—Since the statistics are to be 
obtained from a small number of objects (of the order of 10), it is quite evident 
that one cannot reveal very much statistical detail. Simple assumptions on the 
form of the distributions must be made, and a small number of parameters intro- 
duced to describe them quantitatively. 

If the number density of the objects in the neighbourhood of the Sun is D,, 
the number per unit magnitude range may be assumed to follow the Gaussian 
form D, sini 

™ —(M—M,}'/20" 
L(M ) ae oat e ’ (1) 


where M, is the most probable absolute magnitude and o the dispersion in 
absolute magnitude. We have thus introduced the three parameters M,, Do, oc. 
Visual magnitudes are employed. 

The fall of density with distance z from the galactic plane is approximately 
exponential (1), and can be represented by the fourth parameter, the scale height 


feng D,=D, e-*", 
A,=D,/D,=¢e7*". (2) 


The solar neighbourhood is regarded as being on the galactic plane, and z is 
measured positively in either direction. 

It should be possible to introduce a fifth parameter to represent the rate of 
density decrease as one moves from the galactic centre through the solar neigh- 
bourhood towards the rim of the galaxy. Probably the most satisfactory one- 
parameter formula to represent the change of density with R, the radial distance 
from the galactic centre, is 

R?+e 
Dp=Do rw , (3) 


where R, is the solar distance from the centre. However, no attempt has been 
made to determine the fifth parameter « in the present paper. 

When a is large or unknown M, and D, become rather unsatisfactory para- 
meters for describing the distribution. The objects actually used for the analysis, 
which are those of greatest apparent brightness, tend to be intrinsically brighter 
than M,. Hence the evaluation of M, (and with it Dy) becomes an extrapolation 
towards fainter objects. It would be preferable if one could find a system of 
parameters that are analogous to M, and D, but which give more weight to 
intrinsically brighter objects. Unfortunately no system satisfies all purposes, 
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and it has been found desirable to introduce two further systems represented 
by M, and D, and by M, and D, respectively. The first is somewhat closer to 
observations, but the second is needed for calculations. They are very similar 
and it is easy to convert from one to the other. 

M, is the absolute magnitude of objects giving the greatest emission of visible 
radiation per cubic parsec, and D, is the emission per cubic parsec expressed as 
the equivalent number of stars of magnitude M,. M, and D, therefore have a 
physical meaning that is just as clear as for M, and Dp. 

M, is the most probable absolute magnitude of objects with a luminosity 
function (1) and seen to a finite limiting apparent magnitude m, when the objects 
are distributed uniformly in non-absorbing space. The number seen would 
be D, times the volume of the sphere whose radius is given by the modulus 
m,—M,. This serves to define D,. 

The relations with M5, Dy and o are 

M,= M, —0°4 o?/logio e, 
M,= M,—0°6 o?/logyy e, 
D, =D, exp {—} (0°4 a/logio €)*}, 
Dz, = Dy exp {—4 (0°6 a/logis €)*}. 

These relations may be seen from Fig. 1, where log,, number of, brightness 
of, or volume occupied by objects is plotted against absolute magnitude. The 
curve A represents the number of objects per cubic parsec and per unit of M 
from (1). Dg, is made equal to 0/27 for convenience. 

Curve B represents the luminosity B per object in terms of the luminosity 
M,, and hence dlog,, B/dM=-—o-4. The curve C, obtained by adding the 
ordinates of curves A and B, represents the distribution of emitted light expressed 
as the number of objects of brightness M,. The dispersion of the curve C is 
the same as A, and C has a maximum at M, = M,—0-4 o*/logye. If the C curve 
is displaced to represent distribution of emitted light in terms of M, instead 
of M, we obtain curve D. The total emission D, in terms of M, is given by 
logy9 D = logy) Dp — 0-0807/log,)e as in (4). 

The argument for M, and D, is similar. If V is the volume of non-absorbing 
space within which uniformly distributed objects brighter than apparent magni- 
tude m, can be seen, then dlog,,V/dM=-—o-6. When m,=M, and the unit 
of V isasphere of 10 pc radius, V is represented by curve E. By adding ordinates 
of A and E we obtain F with a maximum at M, = M, —0°6 o?/log,,e. Displacing 
the curve F to represent the case of m, = M, instead of m,= M, we obtain curve 
G, and find D, given by log D,=log D,—0-18 o?/logy, e as in (4). 

If the slope of the line B or E had been —gq it is evident we would obtain a 
maximum at M,= M,—g o?/log,) e. This is used in Section 4. 

3. Observable parameters.—It is intended that the observable parameters 
should be derived from the objects with smallest apparent magnitude. The 
desirable data for each object are the apparent and absolute visual magnitudes 
m, and M,, and the galactic coordinates / and 6. We must reduce these indi- 
vidual observations to suitable mean values that may represent the class of 
objects and can be analysed to give the parameters of Section 2. 

The suggested procedure is to list the brightest objects in order of apparent 
magnitude. About 20 objects are listed if available. If N is the order number 
on this list we obtain two parameters by plotting log N against m,. It may 
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be assumed that the curve should be linear and that departures from linearity 
are due to statistical sampling. ‘The two parameters chosen to describe the 
curve are 

(a) m,(10) the brightness of the tenth object (smoothed), i.e. m, where the 
curve crosses log,, V=1, and 

(b) the slope g=d logy) N/dm,. 

One sees from the definition of m, (10) that results are effectively standardized 
to the tenth object, and that if the list contains more than ten objects the others 
are used for statistical stability only. ‘The theory of Section 4 however could 
be applied equally well to any number of objects. 

Since so few objects are used it is advisable to study the individual points 
on the curves and correct irregularities. Least-squares solutions are not suitable. 
Examples in Section 7 show that g may be influenced by the irregularities but 
that m,(10) is a fairly stable measurement. 

For determining the appropriate mean of the twenty (say) absolute magnitudes 
in the list we must consider that the theory to which comparison will be made 
averages the objects brighter than the tenth. We therefore wish to determine 
mean M, for the brightest ten, but also to let the results be influenced by 
the next ten in order to give greater statistical stability. We must decide 
whether there is, or should be, any systematic trend in mean M, as N varies. 
The number of objects seen within the range M to M+dM is proportional 
to two factors, firstly L(M) the number per unit volume, and secondly the 
volume V through which objects can be seen to the required limiting apparent 
magnitude m, V is a function of the modulus (m,—™), and therefore 
d logy, V/dM = —dlog,, V/dm,= —q. We saw at the end of Section 2 that such 
a distribution has a maximum at M,=M,—qo*/log,,e, which is independent 
of m. Mean M will therefore have no strong tendency to change with N, and 
we use all available values of M to determine the mean M= M,. However, if 
for any reason certain objects have been discarded when determining m, (10) 
and q, it is advisable to omit these objects from the estimation of M. If the 
distribution is Gaussian M=M,= My, — qo?/logye. 

The same absolute magnitudes may also be used for the determination of 
the observed dispersion in M which we call oy, but the values so obtained are 
not very reliable. 

The tabulated values of 6 taken as positive may be used to obtain the mean 
value 5 and the median value 4 of the galactic latitude. In this case there is often 
a distinct tendency to vary with m, for which allowance should be made. Suppose 
that the average is obtained from N objects instead of the standard 10, then the 
limiting m, is m,(10)+(log,, N—1)/¢. The correction Ab(=b,)—4y) may be 
determined from Fig. 2 and is found to depend very little on y or o. With 
sufficient accuracy the correction may be written 


Ab = p|m.(.N)— m,(10)| =p (logy) .V—1)/q, (5) 
where p has the values 


fe Cc 


borb 1 2 5 2... — 
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p ° 3 o°:6 I a7 2 ‘4 2°-8 2°"9 2 5 I 5 o 
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‘he most accurate estimates of b and 6 can be obtained by using all the (non- 
discarded) objects in the list, and, if necessary, making a small correction from 
this table. 
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_ The most appropriate measure of galactic central concentration is the mean 
/ or median / of |(/—J,)|, where /, is the longitude of the galactic centre = 328° 
and longitude differences |(/—J,)| are less than 180°. As mentioned earlier 
there is no analysis for using this observation as yet. _ eae 

The reason for duplicating the observations 5, 6 and /, / is that sometimes 
one figure is more reliable than the other and one should provide the opportunity 
of using either or both as the occasion demands. 

It can be seen that there is a close correspondence between the descriptive 
and the observable parameters as follows: 





Descriptive | M,, M;, Mz; Do, D,, Ds 


o € | 
Observable m,(10); gq; M 7 


If absolute magnitudes are available there are at least as many observable para- 
meters as independent descriptive ones, and it should be possible to obtain an 
estimate of the latter. If there is no information at all on absolute magnitude 
there is, of course, no possibility of determining number density. However, 
the observations still pave the way to the descriptive solution in the sense that if 
either M, or y be assumed the other can be determined. 

4. Equations for calculations—The problem is to compute the observable 
parameters of Section 3 from the descriptive parameters of Section 2. To do 
so one must make an assumption on the form of the interstellar absorption. 
Following Parenago (2) we assume it to be distributed exponentially from the 
galactic plane. Then the absorption in magnitudes per parsec can be expressed 


bb om i, 











a=a,e~*"*, 
where z is the distance from the galactic plane and £ the absorption scale height. 
The total absorption in magnitudes to a distance r in parsecs is 
a(r)=a,B(1 —e~" ™"/*)/sin b, 
where b is the galactic latitude. 
As in Seares’ method of analysis (3) we define ry by 
5 logigto9=5 logyyrt+a(r)=m+5—M. (6) 
The number of stars within solid angle w and unit apparent magnitude range is 
(e.g. Smart (4)) sd 
A(m) =e | 792 Ao(r—)LAm + 5 —5 l0gs0 70)drw (7) 
where A,(7) =A(r)/G(r) expressed in terms of 7. 


A(r) = D,/Dy = density of objects at the distance r relative to density of 
solar neighbourhood. 


afesils = 3a(7) 
G0 = (1+ Spee 3%) (FHiee 3) 
The luminosity function L(M) is given in (1) and it is assumed that the dis- 


persion o is invariable. Changing the variable in (7) to M by (6) and expressing 
A, in terms of the modulus m— M we have 


@W 
5 logy e 





A(m) = 


i 8) 
«o 


10% 10% |" Ax(m—M)L(M) 10-%™ dM. 
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In order to convert to the constants M, and D, of Section 2 we note from Fig. 1 
(curve F) that 
L(M)10- 9%" — Mo) — 7 ¢0°180*/logi0 € a exp[ — (M- M,)*/20"}, 
giving 
4(m) = ogre exp[ —(M-— M,)*/202]dM 
5 logy € s = 
(8) 
The advantages of using M, and D, are that o does not appear among the constants 
outside the integral and that the integral can be easily evaluated with suitable 
accuracy. Calling the integral A, the following simple approximation has been 


used. 
A,(m — Mg) =0-8A,(m — M,) + 0° 1A9(m — My + 20) + 0°1A4(m — M, — 20). (9) 


Silas - 
19ee™m— Ms) FT), | A,(m—M) - 
J —@ 


\/270 


The total number of objects brighter than m, is 
10°w 
5 logy € 


10%w 
= 3 10° 6") DS, 


m™, 
N(m) = 190-06": P), | A,(m — M,)10°%" dm 


where 


|  A,(m—M,)10°%™ din 
S=)-2 





m™, 
{ 106" dm 
= mean value of A,(m— M,) weighted by 10°”. 


Values of S may be determined from these formulae for selected values of 4, y, 
(m,—M,), o. To determine values for the whole sky S must be expressed as a 
function of 5 and integrations carried out. ‘The mean and median values of 6, 
which are denoted b and 4, are given by 


~ en /2 m/2 
b= | Sb cos db/| S cos b db 
Jo 


~0 


-b 7/2 
| S cos 6 db= [__ S cos b db. 
0 J% 


The number of objects in the whole sky becomes, from (10), 


10° { +2 


/2 
N= 7 19°'6™—Ma)D), 27S cos b db. (11) 


—a/2 


For the special case of objects uniformly distributed throughout non-absorbing 
space S=1, and we recover the expected result which defined D, in Section 2: 


N= . 10? 1974") J), (12) 
If we write 


+2/2 
[" amS cos b db, 


/ —2f2 
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we obtain the general resulting formula in a convenient form 
N 3 = 10°10 6" MD, S. ( I 3) 


Comparison of (12) and (13) shows that S can be regarded as the factor by which 
absorption and change of space density affect the number of visible objects. 
Since N is proportional to 10°” S, the slope q becomes 


q=d logy, N/dm, =0-6 + d logy) S/dm,. 


5. Calculations.—Before the calculations can proceed it is necessary to choose 
values for a, and 8. ‘The values adopted for the main calculations are 

ay = 2°5 x 10-* mag/pe (vis.), 
B =150 pe, 
giving an optical half-thickness of the galaxy of a)8 =0-375 mag (vis.). 

Values usually quoted for mean absorption in the galactic plane are about 
0-8 x 10°% mag/pe (§), but van de Hulst has pointed out (6) that for absorption 
by interstellar clouds the amount of absorption along a random sight-line would 
be two or three times greater than this. It is the larger absorption that is 
required for our present purpose. ‘The question of interstellar absorption 
has been studied by Parenago (2), who obtained a mean random absorption a, 
of 3:4 10-8 mag/pe (pg). These two points of view are in approximate agree- 
ment and lead to the adopted value of ay. 

The most useful data for obtaining the optical half-thickness of the galaxy, 
a)B, are counts of extra-galactic nebulae. Hubble’s counts (7) gave a value of 
0°25 mag (pg), but re-analysis of his data by Parenago (2) gave 0-34 mag (pg). 
Other surveys give a greater gradient of nebular number with 6 and indicate a 
greater value of a8. The data by Shapley (8), Shapley and Ames (9), and 
Neyman, Scott and Shane (10) have been plotted against b to obtain the adopted 
value of a8. Parenago’s value of f is 100 pe (2). 

The values used for other constants were 

sin b=0-0, O-I, 0-5, 1°0. 
y =20, 50, 500, oo pe. 
o=0, I, 2 mag. 

m—M from —5 to about + 20. 

For each selection of 6 the quantities r, log,,r, exp(—rsin6/B), «(r), G(r), 
5 logyo7%, and m—M were tabulated with intervals of o-1 in logyr. For 
each y and 4 further tabulations were made of exp(—rsin6/y)=A(r), and 
exp(—rsinb/y)/G(r)=A,(r). By comparison of A,(r) with 5 log,)7) we obtain 
Ag(7o), or Ag(m—M). The argument for further tabulations was m—M, 
or m,—M, in intervals of 1mag. A,(m—M) was integrated by the approxi- 
mation (9) to give A,(m—M,) for o=o0, 1, 2mag. The integration of (10) 
to give S was carried out in two parts. For lower values of m the expression 
A,(m— M,)10°*” was approximately logarithmic and the integration was made 
numerically; for higher m the form changed and the integration was made 
graphically. Finally S was plotted against 5 so that S, 6 and 6 could be found by 
graphical methods. Only four values of 6 were available for this graph and 
some inaccuracy was introduced when there were rapid changes of S with 6. 
This sometimes occurred for the higher values of m,—M,. However, the 
accuracy throughout appears to be adequate for the application of the analysis. 





No. 4, 1954 Whole-sky statistics of celestial objects 395 


The results of the calculations for o=1 are given in Figs. 2 and 3 and in 
Table I. The data for y=100 and 200 pc were obtained by interpolation. 

The values for =o and 2 magnitudes are not given completely, but the 
effect of varying o between o and 2 can be seen from Tables II and III in the 
case of y=20, 50 and 500. ‘The effect is negligible for y= «. 





T T T T T 











10 
Modulus m,—M, 


Fic. 2.—Curves for mean galactic latitude. 











6 








Modulus m,—M, 


Fic. 3.—Factor by which absorption and space density change affects 
the number of visible objects. 

There can be no certainty that the best choice has been made for a, and p. 
‘The computational labour prevents repetition of the whole calculations for other 
values but the general trend of a change in the constants can be seen from com- 
parisons within Tables II and III. In the tables the y=50 values have been 
repeated for (a) ay= 1-6 X 10-3 mag (vis)/pce, B= 150 pe, and (6) ag=2-5 x 10-* mag 
(vis)/pe,8=100pe. From these calculations it is seen that logy) S is not greatly 
affected by a change in a, or 8. Reduction of a, from 2-5 to 1-6 mag/kpc has the 
effect of lowering 6 by 2° or 3° when 6<15°. Reduction of 8 from 150 pc to 
100 pc (to agree with Parenago) has almost negligible effect on 5. 

6. Interpretation.—The aim of the calculations was to enable one to determine 
the descriptive parameters from the observable parameters, and we now study 
Figs. 2 and 3 with this in mind. 
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Suppose in the ideal case that m,(10), g, M, o,, 6, and 6 are all available from 
measurements and are standardized to the roth object. We may assume 
a =ay (see below). M will equal M, as in Section 3, whence from (4) 


M,=M+ go? logy) e= M+ 2°3qg0", ] 


M, = M +(q—0-4)o®/logy) e = M + 2:3(q—0°4)o”, (14) 


M,= M+ (q—0°6)o"/log,,)e = M+ 2°3(q—0°6)o". | 
M is always between M, and M,. Using m,(10)— M, with bin Fig. 2 we obtain 
an estimate of y. Fig. 2 may be modified by the data given in Tables I and 
II if (a) b is used instead of b, (b) o departs seriously from 1 mag, or (c) the adopted 
constants a, and f are not considered suitable. With y so determined one may 
obtain S from Fig. 3 or Tables I and III. 
Putting N= 10 in (13) and using (4) and (14) one obtains 
logo Do = logy Dz + 0°41 50° o _ 
= — 2:62 —0-6m,(10) + 0°6M + 1-38(q — 0°3)o” — logy, S, 
logy) D, =logy) Dz + 0-230? “ 
= — 2-62 —0-6m,(10) +0°6M + 1-38(q—0-43)o?— log, S, [ 
log) Do = — 2°62 — 0°6(m,(10) — M) — logy S 
= — 2:62 —0-6m,(10) + 0-6M + 1-38(q—0-6)o? — logy, S. 





‘Thus all the expected descriptive parameters are obtainable. However, if the 
modulus m,(10)— M, is less than about 3, 6 and 6 become insensitive to y and 
hence y cannot be determined. In such cases y is not needed for the determination 
of S. Acheck on the measured value of g may be obtained since it should equal 
0-6 +dlog,,S/dm, which may be found from Fig. 3. ‘The measured value may 
be amended if necessary and the analysis repeated. 

It is reasonable to put o=o,,.. Although accidental errors in M would make 
ay >c there will be a tendency for individual estimates of M to be somewhat 
dependent. This will reduce o,, and may counterbalance the accidental errors. 
In any case the constants o, and o can seldom be determined with suitable 
accuracy because the distributions may not be Gaussian. For this reason it is 
desirable to choose from (14) and (15) the form for which the coefficients of o? 
are the smallest. It is found for various objects that g can vary over the whole 
range from o-o to 0-6, but it often has intermediate values about 0-3 and 0-4. 
The representations in M, and D, are the least sensitive to errors in o*, and 
as their values have also a clear meaning they are regarded as the main 
results. 

For some categories it will not be possible to obtain all of the observable 
parameters. Some estimate of m,(10), g and 6 will nearly always be possible 
if the category is recognizable at all. If o cannot be estimated it is reasonably 
safe to assume o is between 1 and 2. Ifo is greater than 2 (as for “‘all stars’’) 
errors introduced can be large. The most serious likely omission from the 
observables is M. Fortunately, even though individual values of M for the 
objects concerned may not be available, it often happens that a general estimate 
for the category has been made and can be used in place of VM. If M is quite 
unknown but there are some estimates of y, a reverse procedure may be used 
to determine M. 
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7. Examples.—The examples given in this section are intended as preliminary 
illustrations of the foregoing analysis. Any final studies should be accompanied 
by complete lists of the objects and values used, but the lists used for the present 
examples have not been compiled with the accuracy and completeness that would 
justify their publication. ‘The assortment of objects chosen is varied enough to 
show up the main features and difficulties of the whole-sky statistics. 

Fig. 4 gives the plots of N against m, (or m,, in two cases) from which m,(10) 
and g are obtained. ‘Table IV shows the observable measurements on the left 
and the resulting descriptive parameters on the right. The calculated value of 
9 (Geate= 0°6 + d logy, S/dm) is also added. 

Wolf-Rayet stars.—Data are mostly from Beals (11, 12), but a brighter value 
of —4 is used for M in order to give a more acceptable result for y. Even with 
M= —4 the value of 25 pe derived for y is smaller than might be expected. 

O stars.—The number of O stars is sensitive to the criteria defining the division 
between O and B stars. ‘The data in Fig. 4 and Table IV represent O stars from 
the Henry Draper Catalogue (13). ‘The Og-5 stars of the Yerkes Classification 
(14) are mostly Bo on the HD system. Among the bright stars these Og-5 sub- 
types are more numerous than all other O stars. We might conclude therefore 
that if they had been included the number density of O stars would have been at 
least doubled. 

Flare stars.—This is an example of a category of which less than ten are 
known (15). Since they are near stars the slope q should be o-6, and this has 
been adopted in Fig. 4 and Table IV, although the plotted points would have 
given a much smaller value. ‘The solution in Fig. 4 is equivalent to making the 
assumption that very few of the flare stars with m,>12 have been detected but 
that most of those with m,<12 have been. The two faintest stars on the diagram 
would not then be among the ten brightest and are discarded for M, 6 and oy,. 
This does not leave enough stars to provide a usable value of 5, which, in any case, 
could not give a value of y because m,(10)— M, is too small. 

All stars.—It could not be expected that a category as heterogeneous as “all 
stars’’ could give accurate results by this method. Nevertheless the total star 
density (0-33 stars/pc*) and the luminous emission (0-0013 (M, =o stars)/pc*) 
found by this method are not much greater than values from the luminosity 
function (about 0-18 and 0-0006 respectively). The small value of y found in this 
case is probably due to some flattening of the local system of stars and may not 
represent stars as a whole. 

White dwarfs.—When m, is plotted against N for the brightest known white 
dwarfs (16) one obtains g = 0-35 or less, whereas objects as close as this should have 
q=0'60. It may then be noticed that the three brightest stars « CMa B, 40 EriB 
and « CMiB are all faint companions which would not be detectable at great 
distances. The analysis has therefore been made separately for faint-companion 
white dwarfs and for free ones. For faint companions Fig. 4 shows the slope g 
to be much less than the theoretical 0-6. The theoretical slope is used in this case 
and we must conclude that, unless « CMaB and 40EriB are chance exceptions, 
the faint-companion white dwarfs are very numerous. For the free white dwarfs 
the analysis is straightforward. 

Globular clusters.—Data are from Shapley (17), Graff (18) and Lohmann (19). 
Since the distribution of globular clusters in galactic longitude is asymmetric 
there should be an analysis of /, 7 to give « before one could assign a clear 
meaning to Dp, D, or Dy. 
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Open clusters.—Data are from Graff (18) and Collinder (20). The first four 
clusters in the list were Sco-Cen, UMa, Hyades and Pleiades. The smallness 
of q for the initial plot that includes all these suggests that there is a systematic 
change in the manner of selection with increasing cluster total magnitude. The 
adopted curve in Fig. 4 is obtained by omitting the first two clusters. The curve 
then becomes more normal and is probably more representative of open clusters. 

Extra-galactic nebulae.—Data are from Holmberg (21) and Stebbins and 
Whitford (22) and have been corrected to the new distance scale. Since our 
own galaxy is a member of the local group (23) the nearest nebulae will not be 
representative of all space. In that case a fairly safe procedure is to start the 
list at the point where non-local-group galaxies appear. ‘The first five nebulae 
on the list are Large Mag.Cl., Small Mag.Cl., Andromeda Neb., M33 and 
NGC 5128, of which NGC 5128 is the first non-member. ‘The adopted objects 
for Fig. 4 start one before this at M33. This has approximately the effect ot 
giving unit weight to the local group. 

Normal novae.—Supernovae and recurrent novae have been omitted from 
the list (24). Of the 19 novae listed only 4 occurred outside of the period 
1900-1950, hence 50 years might be regarded as the effective observing period. 
The magnitudes used refer to the intensity maxima. ‘The graph defined by 
the points in Fig. 4 is more curved and the slope smaller than expected. If this 
is interpreted as indicating that many m,>3 novae have been missed we may 
proceed to obtain the results in Table IV. ‘The density obtained of about 10~? 
novae pe-* century ~' is surprisingly high and amounts to more than one nova 
per star in 10° years. ‘This high evaluation of nova density is due to the assump- 
tion that there are many unrecorded novae with m,,,,>3, which assumption is 
however more likely than that there are more close novae than average. It is 
possible that the limitations of the galaxy in the direction of the galactic plane 
have affected this result. 

A stars.—The high value of o,,=2-1 was influenced by the presence of the 
supergiant « Cyg in the list. It was reduced to o=1-6 for the calculations. 

8. Discussion.—The whole-sky statistical analysis is general enough and elastic 
enough to apply to almost any category of celestial object that can be recognized. 
Its weaknesses are inaccuracies caused by (a) the smallness of the number of 
sample objects used, and (4) the simplicity of the assumptions made regarding 
the distribution of objects, their magnitudes and the absorbing medium. 

If 20 objects are used the standard error in Dy, D,, D, from statistical fluctua- 
tions will be 20-1? or about 25 per cent, which might be regarded as satisfactory. 
However, this figure is based on an assumed random distribution, and any 
clustering of the objects that exists will tend to increase this error. Some 
indications of the clustering are generally available from large-scale studies of the 
objects and, furthermore, such effects can be often detected by unevenness in the 
distribution of 6 or / and irregularity in the plot of log,, N against m,. If these 
data have the predicted symmetry or smoothness it may be assumed that the 
results are not adversely affected by clustering. On the other hand, if clustering 
exists one can usually make some allowance for it (e.g. extra-galactic nebulae in 
Section 7). 

The effect of the Gaussian distribution (1) assumed for the luminosity function 
has already been mentioned, in Sections 2 and 6, where it is shown that results 
expressed in the representation M, and D, are not very sensitive to departures 
from the Gaussian form. 

27* 
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The effect of the simplified assumptions on the analysis may be checked 
eventually by comparison of categories studied by other methods, but very little 
direct comparison can be made at present. From Section 5 we find that the 
most recognizable effect of varying a, and f is that for a given y, 5 increases with 
a, or for a measured 4, y increases with decreasing a. The derived y’s for the 
relevant objects (Wolf-Rayet stars 25 pc; O stars 35 pc; open clusters 70 pc; 
novae go pc) have some tendency to be lower than expected and may indicate 
that the adopted value of a, is too large. The adoption of a lower value of 
8 (B=100 pc has been used by Parenago (2)) does not improve the situation. 
More analyses should give clearer indications on these points. 

‘The space densities D, obtained for globular clusters, open clusters, all stars 
and O stars are from 1 to 4 times greater than the corresponding values given by 
Parenago (25). The difference is probably due mainly to the assumption made 
here that the nearer and brighter objects are representative of all space. 

A factor that has considerable effect on the statistics is the fall in percentage 
recognition of objects with increasing magnitude. In the present method of 
analysis this effect can often be appreciated by comparing the observed q with 
q(calc). If the observed g is too small one can assume fainter objects are lost 
from the list by non-recognition. An allowance for the effect can be made but this 
may not be entirely impersonal. If, on the other hand, the observed q is larger 
than q(calc) it may mean that y has been underestimated, since it can be seen from 
Fig. 3 that increase of y (other parameters remaining constant) increases 
d logy) S/dm, and increases g. In principle y could be determined from the 
measured gq, but this is less reliable than by using 5 as in Section 6. 

The effect on the descriptive parameters introduced by changing any observ- 
able parameters may be readily assessed by following through the procedure of 
Section 6. Small changes do not as a rule produce appreciable variation in 
(calc). 

Any systematic comparison of absolute magnitudes and number densities 
of different categories will require a decision as to which of the three systems 
My, Do; M,, Dy; or Mz, Dz should be adopted as standard. ‘The disadvantages 
of M,, Dy have been explained in Section 2, and moreover the usual approximate 
estimates of density are not generally on this system. For example, the density 
of white-dwarf stars in the solar neighbourhood quoted as 7 x 10°/kpc* by Parenago 
(26) does not allow for dispersion in M and might be regarded as an estimate of 
D,. It agrees with D, in Table IV. ‘The discussion in Section 6 suggests the 
M,, D, system to be the best compromise. 


University of London Observatory, 
Mill Hill Park, 
London, N.W. 7: 
1954 February. 
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THE MAXIMUM EFFECT OF CONVECTION IN STELLAR 
ATMOSPHERES ON THE OBSERVED PROPERTIES OF 
STELLAR SPECTRA. I] 
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Summary 

Calculations of the equivalent widths of the calcium lines K and g have 
been carried out for model stellar atmospheres computed in the first part 
of this investigation (M.N., 113, 683, 1953). 

Diagrams connecting the equivalent widths of calcium lines with the 
colour gradients are constructed for both radiative and adiabatic models and 
compared with the same relations established from observational data. This 
comparison shows that actual stellar atmospheres probably are not far from 
radiative equilibrium. 





t. Introduction.—In the first part of the investigations on model stellar 
atmospheres (M.N., 113, 683, 1953), which subsequently will be quoted as 
Paper I, monochromatic fluxes and colour gradients were computed for sixteen 
model stellar atmospheres. Eight of these models were in radiative equilibrium 
(radiative models) and the other eight were computed under the assumption 
that the temperature gradient in the convective zone is adiabatic (adiabatic 
models). ‘The results of these investigations, arranged in the form of diagrams 
connecting the colour gradients ¢ with the reciprocal boundary temperature 
#4 for both radiative and adiabatic models, would settle the question as to the 
amount of energy carried by convection in stellar atmospheres if the boundary 
temperature were accessible to observation: a simple comparison of observed 
colour gradients with the computed ones would then show whether the gradients 
of actual stellar atmospheres approach more to the gradients of the one or the 
other series of models. 

Since, however, the boundary temperature of stars is not known, its replace- 
ment by another quantity easily accessible to observation and _ theoretical 
computations and, on the other hand, closely related to the boundary tempera- 
ture must be chosen as an aiternative way to deal with the problem of convec- 
tion. The equivalent width of the calcium lines K and g seems to be most 
suitable for that purpose. Absorption lines are formed substantially in low 
optical depth, where the temperature does not differ greatly from its boundary 
value; they are therefore closely related to the reciprocal boundary temperature 
#). On the other hand, the equivalent width of the calcium lines K and g shows 
a large range of variation for 9, varying between 0-6 and 1-o which enables us 
to use the equivalent width instead of the reciprocal temperature in spite of the 
low precision of many equivalent-width determinations. 

In the present paper (Section 2) the equivalent widths of the K and g lines 
are computed for the radiative and adiabatic models published in Paper I. 
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The logarithms of the equivalent widths are then plotted against the gradients 
of the corresponding models (Section 3), giving diagrams connecting these 
(computed) quantities for both radiative and adiabatic models; these diagrams 
are shown in Fig. 1. A similar diagram connecting the observed gradients 
and equivalent widths and based on the investigations of 130 stars (Fig. 2) is 
then compared with the theoretical relations shown in Fig. 1. This comparison 
(Section 5) shows that a satisfactory agreement can be achieved between the 
observed diagram and the theoretical diagram for radiative models. This 
seems to favour the conclusion that real stellar atmospheres are not far from 
radiative equilibrium. Finally the possible sources of error are discussed 
(Sections 6 and 7). 

In the way described above the reciprocal boundary temperature 6), used 
as a basic quantity in the computation of stellar models, of their colour gradients 
and of the equivalent widths of their calcium lines, is eliminated in the final 
stage of the investigations. 

2. Equivalent width of calcium lines.—The calculation of line profiles 
followed closely the method of Strémgren (1) as described and illustrated with 
examples by R. v. d. R. Woolley and D. W. N. Stibbs (2). The only change 
concerns the introduction of partition function in the formula for ionization on 
account of the fact that the temperatures of some models are higher than the 
temperature of the model considered by Woolley and Stibbs. This was, how- 
ever, not possible in the computation of the distribution of calcium atoms in 
the singly and doubly ionized states (in the absence of L-S coupling for twice 
ionized atoms); it was assumed here that the ratio of the partition functions is 
equal to the ratio of the statistical weights in the ground states. An error of 
about 20 per cent in the computed equivalent widths is therefore possible for 
hotter models; this is, however, still a tolerable error, changing the logarithm 
of the equivalent width by only about o-1. 

The assumed abundance of calcium in stellar atmospheres is 0-02 of all 
metals. The theoretical values of the oscillator strength (2 for neutral and § 
for ionized calcium) were used in the computations. Since the equivalent 
width is proportional to the square root of the abundance and of the oscillator 
strength, an error in these quantities has the same effect on the logarithm of the 
equivalent width of all models and means only as much as an additive constant. 
This additive constant, however, can be different for the calcium lines K and 
g, since the errors in the oscillator strength may be different for the two lines. 

Apart from radiation damping, only collision damping by neutral hydrogen 
was taken into account. The effect of other sources of line broadening is much 
smaller and probably does not exceed a few per cent of the broadening due to 
hydrogen and natural radiation damping. 

Generally speaking, it may be expected that the logarithms of the equivalent 
widths of calcium lines can be computed with an error of about 0-10-0-20 for 
hot models. For cooler models this error is certainly smaller—disregarding of 
course an additive correction for an error in the assumed abundance of calcium. 

The equivalent width of the K line was computed for all models and that of 
the g line for six models. The equivalent width of the g line for the two hottest 
models, with #,=0-6, is so small that measurements are no longer reliable and 
therefore no comparison between computations and observation is possible. 

The results of the computations are given in Tables I and II. 
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The equivalent width of the g line is not affected seriously by convection 
even in the hottest models for which it was computed; but there is an appreciable 
luminosity effect—for dwarfs the equivalent width of this line is roughly 1-5 
times as great as for giants. 

The K line is more sensitive to the possible effect of convection, but only 
in hot models. In cooler models the effect of convection is small, since the 
absorption lines are formed substantially above the convective zone. In hot 
models the instability begins at a small optical depth and therefore the effect of 
convection is greater; for instance, for 0,=0-6 the equivalent width of the 
K line is approximately 2-5 times greater for adiabatic than for radiative models. 
On the other hand, there is little difference between giants and dwarfs. 












































TABLE | 
Equivalent width of the g line in A 
Dwarfs (logy) g=4°'5) Giants (logis) g=3°5) 
6 Radiative | Adiabatic 6 Radiative | Adiabatic 
’ model model bs model model 
°°7 0073 0-080 
o'8 0°24 0:26 o'8 0°13 O'15 
°'9 0°65 0°67 
1'o 1°39 1°40 are) 089 °'90 
Tase II 
Equivalent width of the K line in A 
Dwarfs (log,y)5 g=4°5) Giants (logy>) g=3°5) 
Radiative | Adiabatic Radiative | Adiabatic 
4, model model 9 model model 
06 0:26 0°67 06 0°22 0°52 
°'7 1°89 2°58 
o'8 6-11 6°51 o'8 4°91 5°89 
o'9 12°6 12°9 
are) 18:2 18-4 1‘o 18-9 19°2 


























3. Relations between computed colour gradient and equivalent width.—P\lotting 
the logarithmic equivalent widths from Tables I and II against the colour 
gradients known from Paper I one gets a set of curves shown in Fig. 1. Thus 
the boundary temperature 6, is eliminated from further investigations. 

It should be noted that for both the K and g lines there is a fairly large 
separation between the curves drawn for radiative and adiabatic models, 
especially for higher temperature. This very fortunate result is due to the fact 
that convection has a considerable effect on the colour gradient but only a minor 
effect on the equivalent widths, since the absorption lines are formed substanti- 
ally above the convective zone. 

On the other hand, there is no great difference between giants and dwarfs, 
the g line for the adiabatic models being the only exception (of minor import- 
ance, however, as further investigations show). This is again very fortunate, 
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since it allows us to disregard any differences in the luminosities of the stars in 
question. 

4. Relations between observed colour gradient and equivalent width.—Plotting 
the logarithms of the observed equivalent widths against the observed colour 
gradients for a number of stars one gets the observed relations between these 
two quantities. These relations are shown in Fig. 2 for both calcium lines. 














Sak Weer 2 
Absolute gradient 





10 . 2-9 
Fig. 1.—Theoretical gradient—equivalent-width relation, 


Radiative models : ————— dwarfs ———— giants. 
Adiabatic models: —-—-— dwarfs +-—-:—-+:— giants. 














1-0 Relative gradient 
Fic. 2.—Observed gradtent—equivalent-width relation. 
we - Weight2 © Weight: @ Weight 4. + Peculiar stars. 
Upper curve : K line. 


Lower curve : g line. 


‘The equivalent widths were taken from publications listed at the end of 
this paper (7-21). The mean value was used where the equivalent width was 
determined by many observers; less weight (4) was given to early determinations 
made with spectrograms of low dispersion. Double stars were included only 
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if their components were of the same spectral type or if the difference of their 
luminosities was at least 2 for components of different spectral types. ‘These 
conditions were not satisfied by two double stars. The K-line determinations 
for six supergiants were rejected, since the gradients of these stars are too high 
(space reddening). Further, eleven early determinations of the g line made by 
E. G. Williams (19) and C. H. Payne (15) were rejected, since they seem to be 
too high (by a factor of about 3) in tomparison with recent determinations by 
other observers. 

The colour gradients of most of the stars are available in the lists published 
by the Royal Observatory, Greenwich (3), by S. C. B. Gascoigne (4) or by 
R. v. d. R. Woolley, K. Gottlieb and A. Przybylski (5). Since all these lists 
form a uniform system of gradients, no corrections were applied. Where the 
gradients were available in two lists, its mean value was used. If no gradient 
was available, the value corresponding to the spectral type of the star was used. 

The observed gradient-equivalent width relation for the K line is well 
defined for gradients over 0-20. The mean deviation of one single star from the 
mean curve is there 0-08. For lower values of the gradient the scattering of 
the stars is larger and probably is due mostly to errors in the determination of 
the gradients. 

The curve for the g line cannot be drawn with great accuracy. Reliable 
determinations of the equivalent width are available in a sufficient number 
only for stars with gradients exceeding 0-80. For early-type stars only three 
scattered points are available and there is a big gap for gradients between 0-10 
and o-80. A straight line drawn in the diagram must be regarded as a rough 
approximation to the relation between the gradient and the equivalent width of 
the g line. 

5. Comparison of the observed and computed relations between gradient and 
equivalent width.—Since the observed gradients are relative, a correction must 
be applied in order to compare them with the computed absolute gradients. 
This correction has been determined by W. M. H. Greaves (6) as 1-1 +0-1 for 
the mean wave-length 5000 A. A further correction must be applied to the 
computed equivalent width, if the assumed abundance of calcium differs from 
the actual mean abundance of this element in stellar atmospheres. On the 
diagram connecting the logarithm of the equivalent width with the colour 
gradient the correction for the abundance of calcium is an additive constant, as 
already stated in Section 2. 

The curve derived from observations must lie between the theoretical 
curves for radiative and adiabatic models, which represent two limiting cases 
of the structure of stellar atmospheres. Using the observed gradient—equivalent 
width relation for the A line drawn on transparent paper and trying to shift it 
on Fig. 1 so as to satisfy the above condition for higher value of the gradient 
(4 >2-0), one readily finds that this is possible only if the observed curve follows 
closely the theoretical curve for radiative models over a wide range. For values 
of ¢ over 1-4 it is possible to achieve this agreement even if the zero-point 
correction for the observed gradients is allowed to take values over a fairly 
large range (between 0-8 and 1-2); for hotter stars, however, this agreement is 
only possible if the gradient correction does not differ greatly from 1-00 (cor- 
rections smaller than 1-00 are, however, impossible, since then the observed 
curve for hot stars lies outside the allowed area). Actually a very good agree- 
ment between the observed relation between gradient and equivalent width 
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and the theoretical curve for radiative models is possible, if a correction + 1:03 
is applied to the observed gradients and a correction —o0-20 to the logarithm 
of the computed equivalent width. The agreement of the relations for the g 
line is then also fairly good; the observed curve for the g line then lies slightly 
below the theoretical one, which is perhaps due to the fact that the ionization 
of calcium in stellar atmospheres may be higher than that computed from 
Saha’s equation. 

On the other hand, the best possible agreement between the observed 
relations for the K and g lines and the theoretical curves for adiabatic models 
needs a correction of 1-40 to the observed gradients, which is certainly too 
high. Even then the agreement between the observed and computed curves is 
not satisfactory, since there are systematic differences at both ends of the 
diagram. The observed relation between colour gradient and equivalent width 
for both K and g lines would then lie partly below the theoretical curves for 
adiabatic models; this is physically impossible. 

The correction + 1-03 for the observed colour gradient, giving a very good 
agreement between the observed and the theoretical curves for radiative models, 
lies within the correction range found by Greaves. This is a strong argument in 
favour of the conclusion that the actual stellar atmospheres are in radiative 
equilibrium. As far as cooler stars with absolute colour gradient over 1-4 are 
concerned this conclusion is valid even if the correction for the observed 
gradient varies within the whole probable range found by Greaves for that 
correction; for hotter stars, however, a considerable deviation from radiative 
equilibrium is possible if the correction for observed gradients differs from 
+1-03. A good determination of this correction would therefore be necessary 
in order to settle the question as to the amount of energy carried by convection 
in stellar atmospheres in the way intended in the present investigations. 

6. Possible sources of error.—The definitive estimation of a possible deviation 
of stellar atmospheres from radiative equilibrium is at present not possible, 
since the quantities used in the present investigations are not free from errors. 
These possible sources of error are briefly mentioned below. 

(i) Models are only models and not real stellar atmospheres. ‘The main source 

of error in the computation of models lies in the mean value of the absorp- 
tion coefficient, which is unknown as long as the definitive structure of 
the model is unknown. ‘The effect of a change in the mean absorption 
coefficient is discussed in Section 7. This difficulty can be solved by com- 
puting model stellar atmospheres in higher approximations. 
The equivalent widths, both computed and observed, are known with 
very low precision, especially for hotter stars. Since, however, the whole 
range of the variation of the equivalent width is large, even an error of 
about 50 per cent is here still tolerable. 

(iii) Observed gradients of early-type stars may be too high on account of 
space reddening. The precision of the gradient determination is not very 
high (mean error about 0-08). 

(iv) Calcium abundance in stellar atmospheres is not well known. An additive 
constant applied to the logarithms of the computed equivalent widths 
removes this difficulty; this correction can be found fairly easily. 

(v) Greaves’s correction for the observed relative gradients is not well 
known. Actually a correction + 1-03 for the observed gradients gives good 
agreement between the observed and computed curves for radiative models, 
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but the final conclusion as to whether stellar atmospheres are in radiative 
equilibrium or not will be possible only when the models are constructed 
in higher approximation and when the correction for the observed gradients 
is known with higher precision. 

In spite of many sources of error it seems to be unlikely that adiabatic 
models can be better than radiative ones. 

7. Blanketing effect—The absorption of radiation in the lines of the Fraun- 
hofer spectrum is an additional source of opacity, which modifies the continuum 
between the lines. Actually the presence of absorption lines in the spectrum 
increases the mean value of the absorption coefficient by an amount known as 


Tase III 


Model stellar atmosphere : 0,=1°:0, logy g=4'5; with blanketing effect 


T 6 logiep log10 Pe 
oo 1000 
o'ol 0'996 3°88 0°09 
0°02 °°993 4°04 0°24 
0°03 0-989 4°14 0°32 
0°04 0-986 4°21 0°38 
0°05 o'982 4°26 0°43 
006 °°979 4°30 0°45 
0°07 0°975 4°34 0°52 
0-08 0°972 4°37 0°55 
0'09 0-969 4°39 0°58 
o'10 0966 4°42 0‘60 
o'12 ©'959 4°46 0°65 
O'14 0°953 4°50 0°69 
0°16 0-948 4°53 0°73 
o'18 0°942 4°56 0°76 
0°20 0°937 4°58 0:80 
0°25 0°923 4°63 0°87 
0°30 o'9gll 4°67 0°93 
0°35 0-900 4°71 0°99 
0°40 o'889 4°74 1°04 
0°45 0°879 4°70 I'lo 
o°50 0869 4°79 1°16 
‘60 0°852 4°82 1°27 
0°70 0836 4°85 1°36 
o'8o0 o'821 4°87 1°46 
©'90 0808 4°89 1°55 
1° 0°795 4°91 1°64 
1:2 0°773 4°93 1°79 
I'4 0°754 4°95 1°93 
1°6 0°736 4°96 2°06 
1°8 0'721 4°97 at? 
2°0 0°707 4°98 2:28 
2°5 0°677 5°00 2°50 
3°0 0653 5 ‘oO! 2°69 
3°5 0632 5°02 2°85 
4°0 0-615 5°02 2°98 
4°5 0°599 5°02 3°10 
5°0 0586 5°03 3°22 
6-0 0562 5°03 3°39 
7° 0°543 5°03 3°54 
3:0 0°527 5°04 3°67 
g'0 O°512 5°04 3°76 
10°'0 0*500 5°04 386 
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the blanketing effect. This cannot be disregarded in cool stars; for instance, 
in the solar atmosphere its contribution to the mean absorption coefficient is 
estimated to be about 20 per cent of the total absorption;. in hotter stars it is 
less important, since the absorption lines are less numerous in their spectra. 

In order to estimate the blanketing effect, a radiative dwarf model, with 
reciprocal surface temperature 6,=1-0 and mean absorption coefficient 1-2 
times the value used formerly, was computed. This model is given in Table III. 
It has a colour gradient of ¢=1-g1 and the equivalent width of the K line is 
17°3A. The same model without blanketing has a colour gradient 2-00 and 
the equivalent width is 18-2 A. The change in the equivalent width is neglig- 
ible, since the lines are formed substantially in the outer layers, whose structure 
is not greatly influenced by the change in the mean absorption coefficient. 
The deeper layers, however, which still contribute considerably to the mono- 
chromatic fluxes and thus have an effect on the colour gradients, are more 
sensitive to changes in the mean absorption coefficient. An increase of the 
mean absorption coefficient substantially means therefore a shift of the model 
towards lower values of the gradients in Fig. 1, and conversely a decrease means 
a shift in the opposite direction. This conclusion is of course also valid for a 
change of the mean absorption coefficient due to other reasons, for instance when 
higher approximations computed according to the formula 


“ Se &, 
k= R | kRdv 
FJo 

are used in computations. Since the slope of the theoretical relation between 
gradient and equivalent width for cool stars is small, inclusion of the 
blanketing effect has little effect on the curve shown in Fig. 1 and little effect 
on the general argument of the paper. On the other hand, a change of mean 
absorption coefficient in the hot stars would require an alteration in our argu- 
ment about the constant connecting Greaves’s relative gradients with absolute 
gradients. 
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Note added in proof.—Recent unpublished determinations of the K line made at the 
Commonwealth Observatory agree fairly well with the curve drawn in Fig. 2. The same 
can be said of the determinations of D. Barbier, D. Chalonge and Nina Morguleff (Ann. 
d’Ap., 4, 137, 1941), converted tentatively into angstrom units; they were, however, 
disregarded in this paper, since such conversion, as the authors themselves remarked, is 
not yet possible with the necessary precision. On the other hand, the determinations of 
J. A. Hynek, Ap. #., 82, 338, 1935 (this reference having escaped attention while this 
paper was being written), disagree considerably with the curves drawn in Fig. 2: the 
logarithms of his equivalent widths of the K line are systematically lower by about 0-3¢, 
whilst those of the g line are higher by approximately the same amount. 


Commonwealth Observatory, 
Mount Stromlo, 
Canberra, A.C.T., 
Australia: 


1954 April 15. 
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Summary 


An analysis of azimuth errors determined from stars “‘ doubled ’”’ with 
the Airy Transit Circle confirms that there is a systematic error in the 
collimation errors adopted for use with this instrument. The existence of 
such an error had been suspected for many years because of the discrepancy 
between clock errors determined with the Transit Circle and with Small 
Transit ‘“‘B’’. Estimates of the magnitude of the error obtained by four 
different methods are in fair agreement. 





1. Introduction.—The existence of a systematic error in the collimation 
errors adopted for the Airy Transit Circle at the Royal Observatory, Greenwich, 
has been suspected for many years. In October 1926 a small reversible transit 
instrument, Transit ‘‘B”’, was brought into use at Greenwich in connection 
with the international programme of longitude determination. A comparison 
of the clock errors obtained with the two instruments revealed a discrepancy of 
about one-tenth of a second, after allowance had been made for the difference of 
longitude of the two instruments and for the different methods used in reducing 
the observations. Comparison with the time signals received from other 
observatories suggested that it was the non-reversible Transit Circle which was 
at fault; the most probable source of error appeared to be in the determination 
of collimation error. The clock corrections derived from the Transit Circle 
observations were not as smooth as those derived from the small transit observa- 
tions. But the R.A.s of stars determined with the Transit Circle showed better 
internal agreement when reduced with the determined clock corrections than 
when reduced with smoothed clock corrections. This is to be expected if the 
collimation is at fault. ‘There was some evidence that the results with the two 
instruments would have been more accordant if a constant collimation error had 
been adopted for the Transit Circle each year, suggesting that the observed 
changes in the collimation error were not real. 

In view of this uncertainty, the Time Service after 1927 July was based on 
observations made with Transit ‘“‘B”. The difference between the two instru- 
ments, referred to in this paper as TC-B, has persisted at about the same value, 
with variations that can be attributed partly to the scatter of the observations 
and partly to the existence of a small annual term. The difference of clock 
errors, TC-B, has been predominantly positive; clock errors are considered 
positive when the clock is slow. If the discrepancy is assumed to be due to an 
error: in the determination of collimation error, it can be shown that the true 
collimation error must be greater than the adopted collimation error by about 
o”-3; -collimation error is considered to be positive when a positive correction 
is required to the times of transit of stars at upper culmination. 
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The hypothesis that the Transit Circle clock errors are erroneous, and that 
the error is introduced in the measurement of collimation, is examined and 
confirmed in this paper. The analysis is based on a method which does not 
require any comparison between the Transit Circle and another instrument, but 
only between ‘‘ fundamental ” azimuths as determined from different stars. 

The collimation error is measured with the aid of two fixed horizontal colli- 
mators placed north and south of the instrument; the process is described 
in detail in Section 7. The level error is determined by finding the micrometer 
reading for coincidence of the movable R.A. wire with its own image formed in 
a pool of mercury. ‘The micrometer reading for coincidence is subtracted from 
the reading for line of collimation and the resulting difference is the level error 
expressed in revolutions of the R.A. screw. ‘This can be converted to arc by 
using the known value in arc of one revolution of the screw. The azimuth 
error of the instrument can be determined only by observations of stars, but no 
azimuth marks are available to control possible diurnal changes of azimuth. 
Circumpolar stars are used for the determination of azimuth, since it is important 
to be able to observe azimuth stars at transit both above and below the pole. 
An azimuth star is said to have been “‘ doubled ” when it has been observed at 
two culminations twelve hours apart, and the azimuth error determined from 
such a pair of observations is here referred to as a “‘ double azimuth”. It has 
the important advantage of being independent of any error in the star’s adopted 
place. 

2. Method of analysis—Suppose that at the first transit the tabular right 
ascension of the star was 7, and that the observed time of transit was O,. “Then 
teu T, —O,=a,A,+6,B,+¢,C,+ Ey (1) 
where a,, 6, and c, are the true azimuth, level and collimation errors of the instru- 
ment, A,, B, and C, are trigonometrical factors depending on the declination 5 
of the star and the latitude ¢ of the observatory, and £, is the true clock error. 
a,, 6, and c, are assumed to be expressed in seconds of arc and (7, —O,) and E, 
in seconds of time. ‘This equation is strictly true only if the observations on 
which it depends are perfect. The effect of diurnal aberration is assumed to 
be included in c, in the usual way. It is assumed that pivot errors are negligible ; 
this point is discussed in more detail in Section 3. 

Similarly, at the second transit of the same star, we have 

T,— Og = a, Ag+ b,B, + gC. + E. (2) 
If we now assume that the azimuth error has not changed in the twelve-hour 
period between the observations, i.e. that a, =a,=a, the azimuth error can be 
determined by subtracting (2) from (1). After re-arrangement we have 
(7 — T2)— (O1 — Og) — (6B) — b2B3) — (4¢Cy — €2C2) — (FE, — E) 

(4,— Ay) | 





‘The quantity (7, — T,) is known with a high degree of accuracy and is independent 
of any small error in the assumed R.A. of the star. (£,—£,) is also known with 
sufficient accuracy, since it depends only on the rate of the clock in a twelve-hour 
period. ‘The quantities O, and O,, A, and A, are also known, so that the azimuth 
error a can be determined. 

We shall now suppose that the adopted collimation errors are both subject to 
the same systematic error Ac, in such a sense that 
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Ac= True collimation error — Adopted collimation error. 
‘The adopted level errors will also be directly affected by Ac; in fact 
Ac= True level error — Adopted level error. 
Equations (1) and (2) now become 


T', — O, = @Ag + (bg + Ac) By + (Cg + Ac)C, + Es, (4) 
where a is still the ¢rwe azimuth error, but 6,, bg, cy, cg now denote the adopted 
level and collimation errors. Combining equations (3) and (4), we obtain 

~~ (7, — T2)— (O, — Og) — (6, B, — bpBe) — (eC, — ¢2C 2) — (FE, — Ey) 
(A, - Ay) 
(B,— By) +(C.— <1] 
+ Ac | 
(4,— Ay) 


(B,— By) +(C,— Cy) 
(A, — 2) ] 








‘True azimuth = Deduced azimuth + Ac 





The formulae for A,, B, and C, are as follows: 


A,=} sec 5 sin (6—8), B,=# sec 5 cos (¢—5), C,=% sec 5. 


These three expressions refer to a transit at upper culmination; the correspond- 
ing formulae for lower culmination may be obtained by substituting (180° —8) 
for 6. 

The quantity (B, — B,)/(A,— A.) is a constant independent of the declination 
of the star, and is equal to tan ¢, or 1-256 at Greenwich; indeed, this result 
holds even if two stars at different declinations have been observed. It follows 
that any systematic error in the determination of level affects all determinations 
of azimuth equally; they have, as will be seen, no effect on the determination 
of Ac by the present method. 

The quantity (C,—(C,)/(A,—A,) is equal to sec 4 cosec 5, or 1-606 cosec 5 
at Greenwich. If now several stars are doubled within approximately the same 
twelve-hour period, the true azimuth will be common to all equations, but the 
deduced azimuth will vary, apart from the accidental errors of observation, with 
the star used, through the factor 1-606 cosec 5. The variation of the deduced 
azimuths with 6 can therefore be used to give an estimate of Ac. 

Suppose that two stars have been doubled and that the deduced azimuths 
are a’ and a” respectively. ‘Then, since the true azimuth is the same for each 


star, we have (a’ —a")=1-606 Ac (cosec 5” — cosec 5’) 


aie (a’—a’") 
= = 1-606 (cosec 5” —cosec 5’)" 
For the azimuth stars used at Greenwich the denominator cannot exceed 0-66 even 
in the most favourable case, and is more likely to be about 0-2. The probable 
errors of azimuth determinations made with the Airy Transit Circle are about 
o”-2 or 0”-3, so that, if « is the probable error of Ac from a single pair of stars, 


e> [(0”-2)? + (0”-2)?]}"2/0-66 
> 0"-43. 
In practice, « is probably at least 1”, so that some hundreds of pairs of azimuth 


stars must be analysed if the existence of a systematic error Ac of the order of 
0”-3 is to be demonstrated conclusively. 





28 
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The impersonal R.A. micrometer was not brought into use until 1915; 
observations prior to this date are not considered suitable for analysis. From 
1915 to 1922 the determinations of azimuth error were based on a list of stars all 
within 34° of the pole. The range of the factor 1-606 cosec 5 among these stars 
is only 0-002, so that this period also must be discarded. A fundamental catalogue 
was finished in 1938 April, and after that date only a very few double azimuths 
were obtained. After 1940, regular observations of double azimuths were 
discontinued and have not been resumed, apart from two short periods in 1953, 
when a special series of azimuth determinations was made for another purpose. 
Thus we are left with only the period 1922-1938; however, in this period about 
a thousand double azimuths are suitable for analysis, so that the data should be 
adequate. 

Between 1922 and 1938 azimuths were determined also by observations of 
two of the stars in the azimuth star list, one observed above pole and the other 
below pole; the observations were not usually more than a few hours apart. 
Azimuths so determined are relatively free from the effects of possible diurnal 
changes in collimation, but are affected by any errors in the azimuth star places, 
and they have not been used in the present discussion for this reason. 

The method of least squares was used to obtain the best estimate of Ac. An 
example of the way in which the calculation was set out is given below. Column 1 
gives the date, column 2 the name of the star doubled and column 3 the deduced 
azimuth. ‘This information has been extracted from the annual volumes of 
Greenwich Astronomical Results. The fourth column, x, is the deduced azi- 
muth minus the mean of the azimuths in column 3. Column 5 is the factor 
f=-— 1-606 cosec5. Column 6, y, is the factor minus the mean factor. Columns 
7, 8 and g give the products x*, xy and y® respectively, each multiplied by 
100 for convenience. ‘The estimate of Ac from these entries alone is 


Lxy/Ly? = + 4°270/3°537= +121. 
The column showing 10ox is needed for the determination of probable errors. 






































I 2 3 5 6 7 8 9 
Date Name of Dodused Factor y 100x" 1ooxy 1ooy" 

star azimuth : 

1922 4 ‘ 

Feb. 13 B Cas +§5°52 | +-25 | —1°878| —-o45 | + 6-25 | —1°125 | +0-202 
a Cas +4:96 | —-31 | —1°933 | —-100 | + 9°61 | +3:-100| +1-000 
y Cas +548 | +-21 | —1°847| —-o1g | + 4°41 | —0-294] +0-020 
5 Cas +5°56 | +-29 | —1.856 | —-023 | + 8-41 | —0-667 | +0:053 
e€ Cas +4°39 | —-88 | —1°797| + :-036 | +77°44] —3°168 | +:0-130 

50 Cas +5°71 + -44 | —1°687| +:-146 | +19°36| +6-424 | +2-132 
Means + §°27 — 1833 





The columns of roox*, rooxy and 1ooy? have been summed for each month 
throughout the period under consideration. The results are summarized in 
Tables I, If and III respectively. Table IV shows the number of double 
azimuths m, on which each entry in Tables I, II and III depends. Each table 
consists of two sections, the first for the period 1922-1930 and the second for 


1931-1938. The rows and columns in each section of every table have been 
summed and the results are shown, 
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3. Possible sources of error.—In this section we examine the assumptions 
inherent in the double azimuth method and consider whether the approximations 
involved could have any significant effect on Ac. The following possible sources 
of error have been considered. 

(a) Pivot errors.—The effect of errors in the shape of the Transit Circle pivots 
is to make the telescope trace out in the sky a curve which departs from a true 
great circle. The effect of a collimation error is to make the telescope describe 
a small circle instead of a great circle. It follows that the effect of pivot errors 
can be regarded as a collimation error which varies according to the direction in 
which the telescope is pointed. Before we can treat Ac as a constant quantity, 
it is therefore necessary to consider the magnitude of the pivot errors. 

The last determination of pivot errors was made in 1924, after new bearings 
(Y’s) had been installed. The Introduction to the 1924 volume of Greenwich 
Astronomical Results gave the results of the investigation in the form of x and y, 
the horizontal and vertical components respectively of the motion of the western 
end of the axis relative to the eastern. For the present investigation, smoothed 
values of x and y have been resolved into a displacement of the telescope in the 
plane containing its optical and mechanical axes, equivalent to a variable collima- 
tion error, and a component which is equivalent to a small rotation of the telescope 
about its optical axis, which can be ignored. The first component is plotted in 
Fig. 1, in which a smooth curve has been drawn through points at every 10° of 
N.P.D. (North Polar Distance). Positive displacements of the curve from the 
zero-line are equivalent to a positive collimation error. 


Jae 120" 60" 0” 300" 


————————————_r- 
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Fic. 1.—Pivot errors of the Airy Transit Circle. 


From a study of this curve it was found that: (i) if pivot error corrections 
had been applied to the double azimuths used in the analysis for Ac, the value 
of Ac would have been increased by +0”-11; (ii) if pivot error corrections 
had been applied to the original observations for collimation error, the adopted 
collimations would have been decreased by 0”-06; (iii) if pivot error corrections 
had been applied to all the original observations of collimation, level and times 
of transit of clock and azimuth stars, the clock errors derived with the Transit 
Circle would have been o*-018 greater on average, so that the TC-B difference 
would have been increased by o*-018. 

Pivot error corrections have not been applied in the past to observations 
made with the Transit Circle, because of their small size and the difficulties 
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involved in their determination. They have not been applied to the results 
given in this paper for the same reasons, but the above corrections can easily be 
applied to results relating to periods after 1924 if desired. 

(b) Variation of azimuth in a period of twelve hours.—It was assumed in 
Section 2 that the azimuth error, a, of the telescope was unchanged after a period 
of twelve hours. Suppose, however, that the true azimuth changes by Aag, i.e. 
a,=a,+Aa. The deduced azimuth is increased by A,Aa/(A,—A,). One 
would expect this correction to be sometimes positive and sometimes negative, 
since it depends both on the sign of Aa and on the order in which the transits 
were observed. The effect of ignoring the correction would then be to increase 
the scatter of the double azimuths without affecting Ac. 

The assumption that Aa can be ignored is not necessarily correct. The 
great majority of double azimuths are obtained by observing the above-pole 
transit in the morning, usually at about 10" U.T., and the below-pole transit in 
the evening, at about 22" U.T. Below-pole transits are much less commonly 
observed in daylight because of the difficulty of seeing stars in daylight at low 
altitudes. Hence, if there is an appreciable systematic difference between 
yo, the true azimuth error at about 10" U.T., and ayo, the error at about 22" U.T., 
the estimates of Ac obtained from double azimuths will be in error. 

An estimate of the size of (@j) — ag) can be obtained by comparing an azimuth 
(not a “‘ double ” azimuth) deduced from two stars observed at about 10? U.T. 
with an azimuth from two stars observed at about 22" U.T. on the same day, or 
the preceding day. Unfortunately, the difference due to Aa cannot, in general, 
be easily separated from that due to Ac, although it is theoretically possible. In 
the period 1916 to 1921, however, the two effects are easy to separate, and the 
azimuths for this period have accordingly been analysed. Azimuths in this 
period were obtained by combining an azimuth star within 34° of the pole with 
a clock star observed at about the same time. All azimuth errors deduced in 
this way are affected by Ac to the same extent, so that any difference between the 
deduced azimuths obtained at 10" and 22" must be due to a real change of azimuth. 

The result of the analysis was that 


(449 — @g2) = — 0"-04 + 0”-06. 


The effect of allowing for this term would be to increase Ac by +0”-04. For 
the purpose of the present investigation, it seemed reasonable to suppose that 
(449 — 2) does not differ significantly from zero and that the effect on estimates 
of Ac is not much greater than 0”-04. The data for the years 1916 to 1921 provide 
some evidence that the sign of (a) — a.) may reverse during the course of a year ; 
if this is true, it may provide an explanation of the large annual term in Ac found 
in Section 5. 

(c) Neglect of short-period terms of nutation.—The azimuth errors tabulated 
in the volumes of Greenwich Astronomical Results were calculated without including 
the short-period terms of nutation in the tabular places 7, and 7,. Their 
omission introduces a small scatter into the deduced azimuths which is without 
any significant effect on Ac. 

(d) Personality.—With the use of the impersonal micrometer, it is unlikely 
that the personality differences between different observers, and between observa- 
tions made by the same observer above and below pole, are sufficiently large to 
have any appreciable effect on Ac. A following error, in which an observer sets 
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the wire on a moving image either slightly in front of or slightly behind the true 
centre, has no effect on double azimuths obtained by that observer. A bisection 
error, in which the observer sets to the right or left of the true image centre, 
irrespective of the direction of motion, would affect the determinations of azimuth. 
It is known that the differences of personality among observers are small with the 
impersonal micrometer, and it must be assumed that the absolute personalities 
are also small. Nevertheless, it was considered advisable to omit from the 
analysis all stars which transit south of the zenith at upper culmination. This was 
because there may be a discontinuity in observations at the zenith introduced by 
changing from the “ feet North ” position usually adopted in observing azimuth 
stars to the ‘feet South” position. For near-zenith stars the observer’s 
position is not always recorded. In practice, those stars with factors f numeri- 
cally greater than — 2-000 have not been used. A star with a factor of — 2-000 
transits nearly 2° north of the zenith, and all the stars used in the analysis would 
be observed “‘ feet North ” in the great majority of cases. 

4. The mean value of Ac; comparison with clock errors.—If we divide the 
grand total of Table II by the grand total of Table III, we obtain the mean value 
of Ac for the period 1922~—1938. The result is 


Ac= +0"-124 + 0"-083. 


This value is of the right sign, but is hardly significant in view of the probable 
error. 

Values of the clock error difference, TC-B, are available from 1927 July 
onwards. A comparison with Ac from double azimuths is therefore possible 
for the period 1927 July to 1938 April. The mean value of Ac for this period is 
(again from Tables II and IIT) 


Ac = +0":240 + 0-096. 
If the collimation errors adopted for the Transit Circle were increased by Ac 


seconds of arc, the TC clock errors would be reduced by 0-25Ac seconds of time. 
Assuming that the TC-B difference is due entirely to Ac, we have 

TC-B= +0-25Ac. 
This result has been quoted in many of the annual volumes. 

The mean value of Ac= +0"-24 means that the average value of TC-B 
over the period considered should have been about + 0°*-060. A direct comparison 
with the mean difference of clock errors would be misleading, since the value 
obtained for Ac is necessarily a weighted mean, depending more on observations 
in the winter months than on the relatively few double azimuth observations in 
the summer. However, the monthly means of the difference T'C-B, which 
are given in Table V, can be weighted in the same way, by multiplying each 
monthly mean by the appropriate monthly value of 1ooy? given in Table III and 
dividing the sum by L100y?. The numbers of determinations of clock error, 
with both instruments, and their accuracy, are such that the weight depends only 
on the doubled azimuths. The result is TC-B= +0*-064, equivalent to a 
collimation error Ac= +0"-26. 

The agreement between the values of Ac obtained from double azimuths 
and from the TC-B difference is therefore very good. If the observations are 
divided into two groups, 1927-1932 and 1933-1938, for example, it may be shown 
that the close agreement for the whole period is to some extent fortuitous, but 
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the agreement is still sufficiently good to confirm that the TC-B discrepancy is 
primarily due to a systematic error in collimation. 

The yearly mean values of Ac, as deduced from double azimuths, are tabulated 
in Table VI. The corresponding weighted estimates of Ac from the clock error 
difference, TC-B, are also given. The probable errors in column 2 were 
deduced on the assumption that the probable error of a single determination 
of azimuth remained at its mean value of about +0”:29 throughout the period 
considered; they are inversely proportional to \/(Zy?) for each year. 


TABLE VI 
Yearly means of Ac 





Year Ac from double azimuths | Ac from TC-B 





1927 (part) Fo'50+o-4I +0°35 
1928 +107 +0°29 +036 
1929 -0:08 + 0°30 +009 
1930 —0'22+0°28 +0°'14 
1931 —1'17+.0°34 +0°14 
1932 +o'60+ 0°44 +0°30 
1933 | Foor +0°35 + 0°34 
1934 -0'22+0°40 +028 
1935 Fo'75+ 0°22 +0°33 
1936 +0°:22+0°34 +0°31 
1937 -1'27+0°56 +0°05 
1938 (part) +1:26+ 1-20 +0°30 




















It will be seen that the values of Ac given in column 2 show variations from 
year to year which are much larger than the variations in column 3. In view of 
the probable errors, however, there is no significant disagreement between the 
results derived by the two methods. Moreover, the values of Ac in column 2 
show a significant correlation coefficient of about +0-83 with the values in 
column 3. 

The correlation between columns 2 and 3 would appear to imply that the 
variations from year to year are real; it is difficult to understand why the varia- 
tions in column 2 should be so much larger than the variations in column 3. 

5. The annual term in Ac.—In Tables II and III the columns have been 
summed to give totals of rooxy and 1ooy? respectively for each month. These 
values have been analysed for an annual variation in the value of Ac. The data 
were divided into two groups, 1922 January—1930 December and 1931 January— 
1938 April. The inclusion of the data for the period 1922 January—1927 June, 
which cannot be compared with TC-B differences, appeared to be desirable in 
order to obtain as much information as possible. ‘The two periods considered 
coincide with those covered by the Second Greenwich Catalogue for 1925-0 and 
the First Greenwich Catalogue for 1950-0 respectively. 

The annual term in Ac was expressed as 


Ac=uy+u, cos 0+, sin 0+u, cos 20+, sin 26, 


where @ is an arbitrary phase angle equal to 30° in January, 60° in February, 90° 
in Marchand soon. A least-squares solution for the coefficients was carried out, 
taking account of weights. The results for the two periods were as follows: 
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1922 January—1930 December 

Ac = —0":035 +0"-714 cos 0+0"-522 sin 6+0":251 cos 20—0”-ogo sin 20, 
1931 January—1938 April 

Ac= +0":389 +0"-532 cos 6+0"-150 sin 6—0”:214 cos 20+0":275 sin 20. 


A comparison of these two results suggests that there is a real variation in Ac 
with a period of twelve months, but the existence of a six-monthly term is doubtful. 
The analysis was therefore repeated, with the terms in 20 omitted. The results 
were 


1922-1930, Ac= +0”:037+0"-689 cos 4+ 0"-509 sin 9, 
1931-1938, Ac= +0":266+0"-579 cos 6 +0":133 sin 0. 


A harmonic analysis of the clock error differences given in Table V yielded 
the series ‘ 

1927 July—1938 April, Ac= +0":245 +0”:049 cos 0 +0”-089 sin 0. 
This result does not confirm the existence of an annual term of the amplitude 
suggested by the double azimuth method. 

6. Other methods of estimating Ac.—In view of the uncertainty in the deter- 
mination of Ac it is of interest to consider other methods of obtaining an estimate 
of its value. Two more methods which have been investigated are discussed 
below. 

(a) From the R-D discordance in R.A.—Regular observations of stars by 
light reflected in a mercury pool were made with the Transit Circle for many 
years. Observations in both coordinates, R.A. and Declination, were taken 
until about 1909, when the measurements in R.A. appear to have been discon- 
tinued. The results obtained by reflected light, R, showed systematic differences 
from those made by direct light, D. The R—D discordance in R.A. can be 
explained as the effect of a systematic error in the adopted level error, since the 
effect of a level error changes sign at the horizon, whereas the effects of collimation 
and azimuth errors do not. 

An error Ac in the collimation will not give rise to an R-D discordance 
directly, but, since the line of collimation is used in deducing the level error 
from nadir readings, it is clear that the R-D difference can be equally well 
explained in terms of a collimation error or a level error. Since there is no 
apparent reason why the nadir readings for level error should be systematically 
wrong, it is reasonable to suppose that the collimation error, Ac, is responsible 
for the R-D discordance. 

In the annual volumes of Greenwich Astronomical Results the value of the 
systematic error, Ad, in level error which could account for the R-D discordance 
has been tabulated for each star observed by both methods; weighted yearly 
means for stars transiting north and south of the zenith have been formed. 
Table VII, giving the values of Ac for the five years 1900-1904, is based on the 
assumption that Ac= Ad. 

The north stars give a mean value of Ac= +0"-376 for the five years and the 
south stars give Ac= +0":397. The combined mean is Ac= +.0”-390. An analysis 
of the north stars, which are all azimuth stars, for an annual term gave 


Ac= +0"-438+0"-111 cos 6+0”"-132 sin 0. 
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Taare VII 
Yearly means of Ac from R-D discordance 





North stars South stars 





Ac Weight Ac Weight 
103 158 
+0°44 92 224 
49 125 
73 131 
46 96 


























The observations could only be grouped into hours of R.A. so that it was necessary 
to assume, in deriving this series, that certain hours of R.A. would always be 
observed in certain months, e.g. 1" o™ to 2" 59™ in November, 3" o™ to 4" 5g™ 
in December and so on. ‘This assumption is not entirely satisfactory, since the 
observations of a particular star are likely to be spread over several months. 
The coefficients of cos @ and sin @ in the above expression are probably too small 
for this reason. 

(b) From an analysis of Greenwich corrections to the FK3.—The mean effect 
of Ac on the determined places of clock stars is very small. Nevertheless, there 
is a differential effect varying with declination. This effect allows us to estimate 
Ac, provided we can compare the places derived with the Transit Circle with a 
system of star places which is free from error. Assuming the right ascensions of 
of the FK3 catalogue to be correct, and the Greenwich observations to be in error 
only as a result of collimation error, we obtain 


FK3 place=TC place + Ac(y —¥). 


y is a factor defined by 
y=C+B+ 314 


and ¥ is the mean value of y for all clock stars. The coefficient 3-1 is the average 
value of (1-256 + 1-606 cosec 8) for all stars used in deriving the azimuth error. 

The above equation allows us to calculate Ac from the tabulated differences 
of (TC—FK3). The TC places are based on observations made with the 
Transit Circle in the years 1922-1940. A least-squares solution for Ac yields a 
mean value of +0”-449; when the cos @ and sin @ terms are included, the series 
for Ac becomes 


Ac = +0":446—0"-031 cos 6—0”-000 sin @. 


It was again necessary to identify hours of R.A. with the months in which they 
were most likely to be observed. 

7. Summary and discussion.—F our different methods have been used to obtain 
estimates of Ac. ‘The mean values obtained, and the periods to which they refer, 
are summarized in Table VIII. 

The values of the coefficients wu», u, and v, derived by the various methods are 
shown in Table IX. It must be noted that the coefficients u, and v, given in 
Table IX for the R-D and Greenwich corrections methods are probably smaller 
(numerically) than their true values. 
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Mean values of Ac; summary 





Method Period 





Double azimuths 1922-1938 
Double azimuths 1927-1938 
TC-B (weighted) 1927-1938 
R-D 1900-1904 
Greenwich corrections 1922-1940 

















It is suggested that two main conclusions can be drawn from a study of Tables 
VIII and IX: 

(i) there is an appreciable systematic error, Ac, in the collimation errors 
which have been adopted for use with the Airy Transit Circle; Ac is probably 
between +0”:2 and +0”:5, although there have probably been real changes in 
its value in the present century ; 

(ii) there is probably a real annual variation in Ac. ‘The coefficients u, and 
v, are probably positive, but not so large as suggested by the double azimuth 
method alone. 

TABLE IX 


Values of the coefficients ug, Uy, V; 





Method Period Ug 





Double azimuths 1922-1930 +9037 
Double azimuths 1931-1938 9266 
TC-B 1927-1938 0'245 
R-D I1QO0—-1904 +-0-438 
Greenwich corrections 1922-1940 + 0°446 























The advantage of the double azimuth method, considered in detail in this 
paper, is that it proves that it is the adopted collimation errors which are at fault. 
The analysis depends only on observations of times of transit of stars, of colli- 
mation and level errors and on the adoption of a clock rate. It is most unlikely 
that times of transit or adopted clock rates are systematically wrong. It was 
shown in Section 2 that an error in the adopted level error affects all determina- 
tions of azimuth equally, and has no effect on the determination of Ac. Thus, 
by a process of elimination, we are left with only the collimation error to account 
for the value of Ac. 

The three other methods of estimating Ac provide valuable confirmation of 
its value, but are not in themselves “‘ fundamental”’. The value derived from 
TC-B clock differences rests on two assumptions, firstly that the tabular places 
of stars used with the Transit Circle are free from error and secondly that the 
determinations of clock error with Small Transit “ B” are perfect. The R-D 
discordance gives a determination of (Ac + Ad), where Ad is a possible systematic 
error in the nadir readings. We can only equate the value of (Ac + Ad) with Ac 
on the reasonable assumption that Ad is insignificant. The analysis of Greenwich 
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corrections rests on the assumption that the FK3 places of clock stars are entirely 
free from any error in R.A. which varies with declination. It is to be noted 
that each estimate of Ac obtained in this paper might be altered if it were possible 
to allow for the effects of pivot errors. 

No satisfactory reason for the existence of a systematic error in the measured 
collimation can be suggested. The routine of observation for collimation error 
in the period under review was as follows. 

The north collimator was first set on the south collimator, the observation 
being made through openings in the cube of the Transit Circle. The Transit 
Circle was then set in turn on the south and north collimators; the mean of the 
readings obtained was the deduced line of collimation. The setting of the north 
collimator on the south collimator was then repeated; any discrepancy between 
the original and final setting was attributed to “ motion of the collimators ” 
and an appropriate correction, usually small compared with the scatter of the 
readings, was applied. The Transit Circle, as designed by Airy and constructed, 
did not havea pierced cube. Toset the north collimator on the south collimator 
the instrument had to be lifted out of its bearings. After it had been in use for 
some time, Airy decided to have the cube pierced, so that the setting could be 
made without raising the instrument, by viewing through the cube with the 
telescope in a vertical position. ‘The design of the cube, with stiffening webs, 
was such that it was not possible to cut a circular hole in the cube. It has a 
number of small apertures, which produce diffraction effects, so that the definition 
of the south collimator wires, when viewed from the north collimator through the 
pierced cube, is always extremely poor. 

Once a week, on Mondays, the telescope was raised out of its bearings so that 
an uninterrupted view of the south collimator could be obtained with the north 
collimator. It was found that the setting of the north collimator with the tele- 
scope raised differed systematically from the setting obtained when the south 
collimator was viewed through the cube. A correction to the deduced line of 
collimation was therefore always applied, on the assumption that the setting 
obtained with the telescope raised was the correct one. The values of this 
so-called “‘ Up minus Down” correction applied from 1900 to 1938 are sum- 
marized in Table X. The corrections to the setting of the north collimator 
wire have been converted to corrections to the deduced collimation errors. 


TABLE X 
** Up minus Down ”’ correction 





Correction to deduced 


Period . - 
collimation error 





1897-1905 +O'41 
1906-1914 +0°27 
1915-1921 +0:°28 
1922-1930 +0°31 
1931-1940 +0'23 














It will be seen that the ‘“‘ Up minus Down ” correction is of approximately 
the same magnitude as Ac; if the correction had not been applied, the adopted 
collimation errors would have been smaller and the value of Ac approximately 
doubled. 
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The illumination used in the measurement of line of collimation was diffused 
daylight reflected into the collimators. In the absence of artificial illumination, 
collimation error could net be measured at night when the star observations 
would normally be obtained. In practice, collimation error was determined 
about once a day, usually in the morning. The mean of the observations for 
several days or weeks was taken and was used in reducing all observations in that 
period. The existence of a diurnal variation of collimation could introduce a 
systematic error into observations of stars made at night, since the adopted 
collimation error would differ from the true collimation error at the time of 
observation. However, when artificial illumination of the collimators was 
adopted in 1944, it was found that any diurnal term was very small and could 
certainly not account for a value of Ac of the order of a third of a second. 

The collimation error is the difference between the observed line of collima- 
tion and the line to which times of transits of stars are referred. It seems most 
unlikely that the determinations of the latter could be in error by the equivalent 
of a third of a second of arc. 

The star observations are corrected for the effect of diurnal aberration by 
subtracting o”-19 from the collimation error. 

No significant correlation could be found between monthly values of Ac and 
monthly means of the measured line of collimation. If any annual term in Ac 
exists, it is of the opposite sign to the annual term in the measured collimation 
error itself, and the true collimation error is probably more nearly constant than 
the measured collimation error. ‘This conclusion has also been suggested in 
Greenwich Astronomical Results for 1927 and 1928, in discussions of the TC-—B 
differences. 

The most probable source of error appears to lie in the determination of the 
setting of the north collimator. The poor definition of the image of the south 
wire system and the existence of an appreciable ‘“‘ Up minus Down ”’ correction, 
referred to above, are unsatisfactory features in the determination of collimation 
error, although it is difficult to see how either could introduce a systematic error. 

It is probable that the collimation errors of a reversible transit circle can be 
freed from systematic error more easily than can the collimation errors of a 
non-reversible instrument; it is a matter of regret that the Airy ‘Transit Circle 
cannot be reversed in its bearings. 

Royal Observatory, 


Greenwich : 


1954 March 4. 





SECOND-ORDER ‘TERMS IN THE FIGURE OF SATURN 
Sir Harold Feffreys, F.R.S. 
(Received 1954 April 6) 


Summary 


The theory of the secular perturbations of Saturn’s satellites requires 
several corrections of the order of the square of the ellipticity, and the 
result is a substantial increase of the estimate of the parameter /. The data 
agree fairly well with the model S, of Miles and Ramsey. The ratio of the 
masses of the Sun and Saturn is 3494°8+ 1-1, but the uncertainty may need 
to be increased to 2-0 to allow for that of the micrometer screw. 





1. In my study of the masses of Saturn’s satellites and the figure of Saturn* 
I adopted from H. and G. Struve the expressions for the motions of the apses 
and nodes to the first order in the constants k and /. Professor D. Brouwer 
has called my attention to the fact that terms in k? have effects larger than the 
standard errors of some of the observed values. Further, though for the first- 
order terms the predicted motions of the apses and nodes are equal and opposite, 
this is not true of the second-order terms. 

Brouwert took the potential due to Saturn in the form 


k, , k , i ) 
u="{a4 (1-3 sin? B) + (1-10 sin? 8 +7 sin’ B) 


and derived the expressions 
du ( gk,? ¥ 10%) 


— =n - an 
dt ef 2 a‘ 
dd 3k. 27k, k, 
— = —n(- = —-—_— +10), 
dt (5 2a °° a 
where w is the mean motion of the satellite and a is its mean distance; they are 
connected by 
n®a = u(1 + 3k,/a?). 
In comparison with the notation of the Struves 
ko=tk ; ky=H. 
Then 
dw ‘ ‘ 
= ="\|-— ) 
dt 2 (6) 
dé 
77, (7) 
If it was not for the term in &, in (4), the coefficients of & in these formulae 
would be proportional to m**. Those given by the Struves did not satisfy 
this condition accurately, and I was unable to find out whether they had made 
*H. Jeffreys, M.N., 113, 81-105, 1953. 
+D. Brouwer, A.7., 51, 223-231, 1945. 
29 
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a correction for k in estimating a, or whether they had used observed values 
ofa. Direct examination of the individual values of a determined by G. Struve* 
from observation in Hefte 2 and 3 showed systematic differences from his 
adopted values, and the only safe method was to redetermine the values of a. 
The & term in (4), if not allowed for, would introduce other errors of order k? 
in the calculated values. 

The observed mean distances showed discrepancies definitely larger than 
the apparent uncertainties would account for, and I took simple means. The 
results were as in Table I, in seconds of arc. 


TABLE I 
Adjusted 
26°843 
34°438 
42°633 
54602 
76-255 
176°763 
515125 


- 
° 


Mimas 26°848+ 0-011 
Enceladus 34°441+0-018 
Tethys 42°642+0°013 
Dione 54°608+ 0-014 
Rhea 76:260+ 0013 
Titan 176°716+ 0°064 
Iapetus 515°046+ 0:086 


~ 
MW AW Os! 


The Sun reduces the mean motion of Iapetus by 2-46 x 10~* deg./14, and this 
must be added to the actual mean motion to give that due to Saturn alone. For 
other satellites this correction is negligible. Precession is smaller still. Then 
these data, substituted in (4), give a set of equations for 1 andk. On inspection 
it was clear that they could add nothing of importance to the determination of 
k, and I took k/a,?=0-0243, with a,=8"-75, and estimated p. 
= (28149 + 9)10°(1°/19)°(1"). 

x? was 2:6 on 6 degrees of freedom. ‘The adjusted values of a are given in 
Table I. 

With these values 1 recomputed the coefficients of k/a)” and //a,* in the 
motions of the apses and nodes. ‘These were as follows : 

Mimas Enceladus Tethys Dione Rhea Titan 
k/ag? —-1°48260 X10! 6°1955 X10 2.9342X10° 1°2338X10*% 3°8328x10° 2°0207XI0 
lag 3939 SX 10® = r-00c0 X10* 3:090 XI0% 7-921 XIO 12°62 071238 

The coefficient of k?/a,* is — 0-2 or —0-6 of that of //a)', according as it refers to 
the apse or the node. 

The largest effect of the terms in k® is +1°-40/1¥ for the node of Mimas. 
I revised the previous solution, taking into account the additional equations 
on p. 95 for the mass of Titan, and found : 
Ho = FO1Z L155 wy = +108 +021; py= +13 4155 y= —3°542073; 
k ; l 
72 = 10 *(2'4303 + 0°0038); =] = 10 “(14°56 + 1°47); 
ay” ay 


lk? =2-47+0°25; 


My, =(3°2 + 3°8)10°* | my, =(2°4115 + 0°0182)10°'; 


I My, = 4147 + 32. 
x”=9°6 on 7 degrees of freedom. 





* G. Struve, Veréff. Berlin-Babelsberg, 6, Hefte 1-5, 1933- 
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The observed and computed secular motions are as in Table II. 


TABLE IT 


Calc. Obs. O-C 

Mimas apse 365 86 365 ‘60 —o'26 +0:20 

node 364°93 36523 +0-30 +020 
Enceladus apse 152°51 
Tethys node 72°231 72°227 —o0'004 +0°065 
Dione apse 30°72 30°75 +0°03 +0°42 
Rhea node 10°05 10°20 +o-15§ +0°08 
Titan apse 0°4994 o'5012 -+0°0018 + 0:0070 

node 0°4994 0°492 —0°0074+ 0'029 


©0008 has been included for Titan to allow for the effect of Iapetus. 

2. The mass of Saturn.—The above value of » for Saturn, compared with 
the value of n?a* for Saturn’s orbit, gives a ratio (2-°8650+0-0009)10-*. This 
needs two corrections. Hill* gives the effect of Jupiter on Saturn’s mean 
distance, for given mean motion, as 


5 logy) a’ = 0-00018522. 


Other planets raise this to 0-00018538. This must be eliminated by multiplying 
y. for (Sun+ Saturn) by 1-001284 ; then we have for the mass-ratio 


Sun/ Saturn = 349404 + 1-11. 


This agrees with H. Struve’s value 3495-4. Hillf found 3502-20+0°53 from 
observations of Jupiter, but remarked that the satellites gave a smaller value 
and therefore continued to use Bessel’s 3501-6. lapetus provides the chief 
equation in the estimate of « for Saturn, so that the mass of Titan has effectively 
been added to that of Saturn. For this reason the above ratio should be increased 
by about 0-8. Systematic error in the measurements of satellite distances is 
possible. Newcomb stated that the calibration of a micrometer screw could not 
be trusted within 1 part in 6000; this alone might give an additional uncertainty 
of 1-7 in the mass-ratio, which may therefore be taken as 3494°8+2-0. H. G. 
Hertz} has recently given 3497-64 + 0-27 trom perturbations of Jupiter; only an 
abstract is yet available but it appears that there is now little disagreement 
between results from Jupiter and the satellites. 

3. Structure of Saturn.—In attempting to compare with the models given 
by Miles and Ramsey§ I found an inconsistency. I had taken m=0-1353. 
Their data correspond to m=o0-162, with an uncertainty of about 3 per cent. 
This was traced to H. Struve’s original solution, which took the period of rotation 
as 10"14™; I had taken 1038". ‘These are observed periods for two spots in 
different latitudes, and the difference would affect m by 10 per cent. Struve 
took the equatorial angular diameter at mean distance as 17”-50, whereas Ramsey 
followed Russell, Dugan and Stewart, who gave 17"-25. Further, the definition 
of m is different, since Struve did not allow for the factor 1—e. For the last 
reason alone Struve’s ellipticity 4 derived from the perturbations of the 


95 


satellites should never have been trusted within 10 per cent. But the formula 


* G. W. Hill, Astr. Papers for the American Ephemeris, 4, 20, 1890. 
+ G. W. Hill, Astr. Papers for the American Ephemeris, 7, 17, 1895. 
tH. G. Hertz, A.7., 58, 42, 1953. 

§B. Miles and W. H. Ramsey, M.N., 112, 234-243, 1952. 
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(16) on p. 241 of Miles and Ramsey’s paper depends on de Sitter’s definition of 
m and is inconsistent with the adoption of Struve’s ellipticity as a datum. 
Having regard to the uncertainties I now take 
m= 0-147 + 0-006. 
If a, is the equatorial radius, 
Jaz=k; Dat=*% 1. 


12 
The present determinations of k/a,? and //a,‘ refer to a,=8"-75, whence, with 
a) = 8" -625, 
J =0°0250, D=0-00449 + 0°00045. 
We have also 


J=e—im+e(— be+im), D=ie—- imet+4n, 


whence e=0-101 + 0°003, K =0°00149 + 0:00025. 

Since e and m are fairly well determined, the most convenient way ot 
comparing with the models seems to be to use Miles and Ramsey’s equation (16) 
to determine 7,, which depends very directly on the model. ‘This gives 

= 1657 + 0-042. 
Their model S, gives 7,= 1-67, so that the value of J corresponds to a body at 
least as strongly condensed to the centre as S, (17 per cent of the mass in the 
central particle). ‘They preferred S, (y,= 1-34); for convenience of computation 
I used S, (7,= 1-04), but it is clear that S, is better than either. 

I computed « for S, and S,; the surface values were 0-00185 for S,, 0-00134 
for S,. This favours S, slightly, but the evidence from 7, is stronger. 7, could 
be used to give an estimate of the central condensation and hence of «; but 
since a large part of the uncertainty of the relation between / and « comes from 
that of m, little would be gained by using the calculated value of « to re-estimate /. 

All the models used by Miles and Ramsey, except J, and S,, have a massive 
particle at the centre. It is worth while to notice that for such models direct 
study of Clairaut’s equation makes e behave like r* for r small, so that y—-3._ But 
the Radau approximation would make 7->2°5. 


160 Huntingdon Road, 
Cambridge : 
1954 April 5. 





THE INFLUENCE OF STELLAR RADIATION ON THE RATE OF 
ACCRETION 


L. Mestel 
(Received 1954 March 31) 


Summary 

Bondi’s discussion of spherically symmetric accretion (1) is extended 
10 take account of the heating of the cloud by ultra-violet radiation from the 
star. It is shown that if the accretion rate is to be high, the radius of the 
HII zone must not be more than about one thousand stellar radii. The 
absorbing power of a cloud is enormously increased by the high densities 
near the star, so that a cloud of moderately high density at infinity is able to 
keep the H II zone small enough for rapid accretion to take place. For a 
star of given type, there is a lower limit to the cloud density that will keep 
the H IT zone small; when a star, accreting from a cloud of given density, 
reaches a critical mass, the H II zone increases rapidly to Strémgren’s value, 
and accretion effectively ceases. ‘The conditions under which the state of 
rapid accretion can be set up are investigated : it is shown, for example, that 
this state can be reached with the cloud of density 10° particles/em* and 
temperature 100 deg. K, provided the star entering the cloud is of mass less 
than 1°75 





1. Introduction.—The first discussions of accretion treated the problem purely 
dynamically, it being assumed that radiation losses would keep the thermal energy 
low enough for pressure to be negligible. ‘The rate of accretion was found to be 
limited by the velocity V of the star through the undisturbed cloud (2) : 

dM 
dt 


where M is the mass of the star, ¢ the time, p, the density of the cloud at infinity, 
G the gravitational constant and « a constant of order unity. 

More recently Bondi (1) has discussed the limiting case of spherically symmetric 
accretion, in which a star is considered at rest in an infinite cloud. In this case the 
rate of accretion is controlled by the thermal energy of the gas: 


oF . ape). (2) 
dt 
where c, is the velocity of sound at infinity and A a constant. ‘The maximum 
value of A (A,) depends on the amount of heating by compression suffered by the gas 
asitflowsin. Ifthe gas contains powerful radiators, or if y, the ratio of the specific 
heats at constant pressure and constant volume, is near unity, the flow will be 
nearly isothermal, and A,=1-12. In adiabatic flow of a gas with the greatest 
possible value of y (2), A, = 0°25. 

By comparing (1) and (2), Bondi was led to conjecture for the intermediate case, 

when c,, and V are of comparable magnitude, an accretion rate of order 
dM 

dt 


‘This has been confirmed by Dodd (3) for the case of supersonic isothermal flow. 


= 2m (GM) V-3p,,, (1) 


~27(GM)?(V2+¢,.2)-3%p.. (3) 
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McCrea (4) has applied (1), (2) and (3) to a specific cosmogonical problem: the 
origin of O and B type stars. ‘The lifetime of these stars as determined by their 
rate of hydrogen consumption is much shorter than the galactic lifetime, so they 
must be the result either of late condensation or of growth by accretion of stars of 
moderate mass born earlier. If the explanation by accretion is correct, McCrea 
points out that the mechanism cannot be purely dynamical, as summarized by (1). 
In order that accretion should be of importance, V cannot be much greater than 
1km/sec; with larger velocities, the increase of stellar mass from (1) would not 
affect the star’s evolution unless the gas clouds were much denser or much more 
extensive than those observed. But the velocity of sound in ordinary neutral 
interstellar clouds ~1 km/sec(corresponding to a temperature = 100 deg. K (5)), so 
that the accretion rate must be governed by formulae (2) and (3). McCrea shows 
further that a star entering a cool dense cloud of reasonable extent (density 
= 103 particles/cm*, temperature ~ 1oodeg. K, radius ~ 10 parsecs) with a velocity 
ofabout 1 km/sec will be reduced effectively torest by the retarding force of the gas 
before it has penetrated through the cloud, so that ultimately the accretion process 
will be spherically symmetric. ‘This result removes the necessity for clouds of 
great dimensions. Finally, McCrea remarks that the small number of stars with 
velocities as low as 1 km/sec is not necessarily a disadvantage to the theory, as the 
proportion of O and B type stars is in any case very small. 

Since this work, the case of spherically symmetric accretion can no longer be 
regarded as a mathematical problem of little astronomical interest, but may itself 
turn out to be the most important accretion mechanism. It therefore deserves 
closer study. 

Bondi’s results show that as long as the variations in temperature in the gas cloud 
are due solely to compression of the inflowing gas, the effect on the accretion rate is 
not larger than a factor of about 4. However, since the work of Stroémgren (6), 
it is known that a bright star is surrounded by a zone of almost completely ionized 
hydrogen (an “‘ HII zone”’) at a temperature nearly the same as that of the star’s 
surface, and it has been asserted (7) that this heating of the gas near the star will 
inevitably cut down the accretion rate to a negligible value. ‘This possibility is 
examined quantitatively in the present paper. In Section 2 the maximum 
accretion rate is found as a function of the radiusr of the HII zone. ‘The ionization 
of hydrogen ceases once the bulk of the ultra-violet light (of frequency greater than 
the threshold frequency for hydrogen ionization) has been absorbed. ‘Thus the 
value of r depends on the properties of the star (principally its surface temperature) 
and on the densities in the cloud. In Section 3 the absorption of ultra-violet light 
is analysed, and in Section 4 the relation is found between the maximum rate of 
accretion by a star of given temperature and the density at infinity of the cloud in 
which it is embedded. 


2. The variation of the accretion rate with the radius of the H II zone 


2.1. The equations of flow.—The problem is idealized as follows: a star of 
mass M and radius R is at rest in an infinite spherically symmetric cloud, consisting 
mainly of hydrogen, which is at rest at infinity. A sphere of radius ¢ (zone II) is 
almost completely ionized, and is at aconstant temperature of 104 deg. K, while the 
rest of the cloud (zone I) is un-ionized and at 10? deg. K. Once again c,, and p., 
represent the velocity of sound and density at infinity, while p, p, v, all functions 
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of the radial coordinate r only, are the pressure, density and inward velocity. 
A steady state is assumed to have been reached. Further discussion and justi- 
fication of these assumptions is postponed until later. 

The treatment follows closely that in Bondi’s paper (1). The equation of 
continuity for steady spherically symmetric flow is 


A=4rr' pv, (4) 
where the constant A is just the accretion rate. 
For each zone the equation of motion integrates to give Bernoulli’s equation 
(dp GM 
hut + oe x —— = constant, (5) 
Jp 


As each zone is isothermal, we have 
Pp - Kp, | (6) 
p=Kzp, J 
in zones I and II respectively, where K, = 107%, K,=2 x 10'# and & is the gas 
constant. ‘The factor 2 in K, arises from the contribution to p from the free 
electrons in zone II. The velocity of sound is constant in each zone because of the 
assumed constancy of temperature, and equals \/ K, in zone I and \/K, in zone II. 
For zone I, (5) reduces to 
$v? + Kj ln x. 
Pew 
on using the condition of zero velocity at infinity. In zone II we have 


Jot+ Kylnp— 2™ = 45,24 K,Inp,— A (8) 


GM 
=o 


where bars refer to conditions at the interface and suffixes to zones I and II 
respectively. 
At the interface the conditions of fit are 
py 0, =p2Ve (9) 
Py + Py Oy? = py + P20”. (10) 
Equation (9g) expresses the conservation of mass and (10) equates the change of 
momentum per second across one square centimetre of the interface to the 


difference of the pressures acting. 
The following transformations are made: 


” GM = K, , oe Ky II 
T= 28 x, v= Co KR yV, p= K.?= z, ( ) 


where the sufhix 7 takes the values 1, 2 in zones I and II respectively. Then (4) 
reduces to 


and 


xtyz=A, (12) 
where 
A=47\(GM)P ca pos (13) 
identical with Bondi’s equation (2). Equations (7) and (8) become 


iy? +Inz—— =o, (14) 


K a\*_,-—, Ky, - 1)? 
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The interface conditions (9) and (10) reduce to 


121 = Vo, (16) 


‘ “ - K, . 
3,(1 +5)=3,(r4 8), (17) 


which combine to give 





. +9" o. +¥a"{Ke Ky) =20, (18) 
M1 Ve 
say. 

2.2. Conditions at the interface.—In order that y, and y, should exist to satisfy 
(18) for given », the condition »* > K,/K,(>1) must hold. When this is satisfied, 
it appears that there should exist two possible values of y, for given y2, and similarly 
two values of y, for given y,. However, the choice is in fact restricted to values 
giving either subsonic or supersonic flow on both sides of the interface.* In the 
narrow zone in which the temperature rises steeply from 10? deg. K to 10* deg. K 
the velocity of sound increases monotonically. The condition (10) holds over 
this region (neglecting the small variation of gravity) so that 

p+ pv? =p(c? + v*)=const. (19) 
Combined with the equation of mass conservation 


pv =const., (20) 


. «@ 
c(£ +=) <const. (21) 
a Tia 


this gives 


€ 
where c is the velocity of sound. Thus from the graph of w/c against c it 
is seen that in a region of monotonically increasing c the velocity must be either 
subsonic or supersonic over the whole range. In our notation, v;=c, is given 
by y,;=1 and v,=cy by vo= \ K,/Ky. Hence if y,<1, yo<1/K,/Kg, and if 
Yi > 1,22 VKy/ Re. 


From (18) - K.\./) IK 
val fee ~tnt al (e- E)} $4 


yi=7+t V(y?-1). (23) 
Thus either the positive or the negative sign must be chosen on both sides of the 
interface. 
2.3. Conditions in zone I1].—From (12) and (15), 
K K K, 1) 
1,2_ 1 a ae on Bed 
(19 K, iny) K, (2inx+ K, oJ +a, (24) 
where « is a constant. 
The expression (2In«+(A,/K,)1/x) has a minimum at x,= $(K,/K,), and if 
x >x,, this minimum will be attained. Also, ($y? —(K,/K,)Iny) has a minimum 
for y?=K,/K,. Hence, if x >x,, the constant « cannot be less than 


BS, Ky\_ %, =z) 
oR (1+in x)- a (2+2insg ° (25) 
For K,/K,=2 x 10%, this gives 

(Jy*—5 x 10-*Iny)—(10-#In + *5 10") 6:57 x10 (26) 





* The writer is indebted to Professor Cowling for drawing his attention to this limitation. 
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for all y in the range x(star)<x<x. In particular, at x equation (26) provides 
a new restriction on the allowed values of y,: either 

, K 
V22 ) > lz 


Ve < Y, ~ Jz . 


The different possibilities are summarized in Fig. 1 








It is seen that for y,= Y, or yy= Y,, the solutions pass through the point C 
for which the gas velocity equals the velocity of sound. The initially supersonic 
solution from A may either pass smoothly into subsonic flow or, by means of 
a finite jump in the acceleration at C, remain supersonic up to the surface of 
the star; similarly, the initially subsonic solution from B may either become 
supersonic or remain subsonic after passing through C. All other solutions, 
such as those starting from F and L, remain either supersonic or subsonic for 
all x in the range. 

If x<x,, there is no limitation on the range of y,; (2In x +(Ky/Ky )1/x) 
attains its minimum at x, and the smallest possible value of « is given on 
substituting x for x and \/(K,/K,) for y in (24). Solutions that are either super- 
sonic or subsonic at x remain so up to the surface of the star; for x=x,, 
there are two solutions branching from the point ~,= \/(K,/K,), one into the 
subsonic and the other into the supersonic region. 

2.4. Conditions in zone I.—The solution in zone I must satisfy the condition 
of zero velocity at infinity. From (12) and (14) we have 


Ind=(Z24 Inx)~(Jy*—Iny) (29) 


for all x from x to infinity. Suppose first that x<o-5. Then at x=0-5, 
(1/x+2I1nx) attains its minimum; ($y?—Iny) has a minimum value for y=r. 
If for «=0-5, y<1, then y remains less than 1 for all x in zone I, and similarly, 
it y(o-5) >1, then y>1 throughout the range; thus only the subsonic sets satisfy 
the condition y=o at infinity. If y(o-5)=1, then again there is a branching of 
solutions, and the flow that is subsonic from x= 0 to x=0-5 may, by a discon- 
tinuous jump in acceleration, remain subsonic, or may pass smoothly inte 
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supersonic flow for x<x<o-5. For this critical case Ind attains its maximum 
value 0-1138, while the other solutions correspond to smaller values of A. The 
yun of the solutions is shown in Fig. 2. 


¥7 

















Fic. 2. 


When x 20:5, (1/x+2Inx) increases monotonically outwards, and so there 
is no possibility of a flow with a supersonic velocity at x becoming subsonic and 
so vanishing at infinity. Thus if x>0-5, y, is necessarily <1. 

To sum up: for x >0-5, y, may take any value up to unity. For x<o:5, it 
is allowed just one supersonic value, namely the value at P, the point where 
the smooth solution through M meets x= x (see Fig. 2). The subsonic values 
allowed must all be less than the value at Q, where the other branch of the maximal 


solution cuts x=x. The value of Ind, and hence the accretion rate, is fixed by 
(29) when y, is known. 

2.5. Solutions for the complete range.—The results of the last three sections 
must now be combined so as to determine what solutions satisfy all the conditions 
at the interface, in zones I and II, and at infinity. The following cases arise: 


(a) x 205 


By (27) and (28) either 


The supersonic cases may be rejected at once, as the corresponding values 
of y, must also be supersonic, by the condition of fit, and so cannot, as shown 
in 2.3, lie on a solution satisfying the condition at infinity. 

Any solution with y,<Y,, satisfies all the conditions. The corresponding 
value of y, is found from (18) to be ~ yg, as the velocities are so much below 
the velocity of sound that (10) reduces effectively to equality of pressure at the 
interface. ‘The solutions behave as in Fig. 3. The largest value of A, A,, corresponds 
to either of the two solutions oBCD and oBCE, while smaller values 
characterize solutions such as oFGH. A set of values is given in Table I, assum- 
ing K,/K,=2x 107. These should be compared with InA,=0-1138, Bondi’s 
maximum. In the limiting case x= «©, when the whole cloud is at a temperature 
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104 deg. K, A, is reduced by (K,/K,)°”, due to (K,/K,)*” from the increase 
in c,%, and (K,/K,) from the decrease in density. It is readily verified that (2) 
gives precisely the accretion rate in Table I, when c,,* and p,, have these new 
values substituted. 
TABLE | 
In A 
II‘I5 
~I2°1§ 


13°14 











(b) 0075 <x <0°5 





For this range the solutions behave similarly to the set (a). Again there are 
solutions in zone II for y, < Ys, with corresponding values for InA. For y.= Yz, 
there are two solutions, «0 BCD and «BCE (cf. Fig. 4), the first being supersonic, 
the second subsonic near the star; for y.<Y, the solutions are everywhere 
subsonic. The condition of fit at x again reduces effectively to continuity of 
pressure, so that y, yp. It is found that within this range of x, the point P 
(corresponding to a supersonic velocity) gives values of » not satisfying the 
condition »? > K,/K, of Section 2.2, so that only solutions subsonic in zone I 
exist. As zone I now includes the point x=0-5, the solutions in zone I all have 
a maximum. The lower limit to this range of x is fixed by the coincidence of 
the points Q and B, giving for Ind, the value 0-1138, Bondi’s maximum. With 
our chosen value of K,/K,=2 x 10”, this limit is 0-075, as in the heading. For 
different values of K,/K,, the limit varies slightly; e.g. for K,/K,= 107, 4 x 107, 
it is 0-087, 0067 respectively. Typical values of Ind, are as follows: 


TABLE II 
In A, 
~II'I5 
3°125 
o'r135 


(c) 2-45 x 10 3<x <0-075 





The diagram for this range is similar to Fig. 4, except that Q falls below B. 
Thus all the solutions remain subsonic right up to the star. Again there are no 
solutions passing through P. The maximum allowed value of InA is again 
0-1138, corresponding to the curve «oMQ(star). 
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(d) x<2:45x 10% 


The subsonic solutions are similar to those in case (c), but in addition there 
is now a solution ooMPR(star), as the value of 7 corresponding to P satisfies the 
condition 7? > K,/K, (cf. Fig. 5). Again, theupper limit x = 2:45 x 107% is the value 
corresponding to K,/K,=2 x 107; for all K,/K,, this upper limit is slightly less 
than x, = 4K,/K,, so that all the zone II solutions are monotonic (cf. Section 2.3). 
For oMPR(star) and oMQ((star), Ind is again 0-1138, with correspondingly 
smaller values for all the other subsonic solutions. 



































3. The theory of radiation heating 

3.1. Physical discussion.—It was Eddington who first drew attention to the 
heating effect of stellar radiation on the interstellar gas (8). Ultra-violet radiation 
of frequency greater than the photoelectric threshold for hydrogen will liberate 
electrons from the hydrogen atoms in the cloud, the mean kinetic energy of the 
expelled electrons corresponding to a kinetic temperature about two-thirds the 
temperature of the exciting radiation (assumed to be given by a black-body 
distribution). The electrons suffer many Coulomb collisions with each other 
and with the heavy particles in a time short compared with the time for recom- 
binations to occur, so that the particles take up a Maxwellian velocity distribution 
with a definite kinetic temperature at each epoch. The slower electrons are 
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much more easily re-captured than the faster, so that an initially cool cloud is 
gradually heated up by the net transfer of ultra-violet light from the central star 
into kinetic energy. 

Neglecting all reactions but the photoelectric ionization of hydrogen and the 
recombination of electrons and protons, the ultimate steady state will be 
determined by conditions of statistical but non-thermodynamic equilibrium. 
The thermal radiation from the star liberates electrons at a rate depending 
primarily on the temperature of the radiation and its dilution; recombinations 
occur to both the ground and excited states of the hydrogen atom, removing 
electrons from the gas and transforming some of the kinetic energy into photons. 
The temperature of the gas and the degree of ionization are determined by a 
balance between gains and losses of electrons and kinetic energy. As the 
radiation field is so dilute, and the lifetimes of the excited states of hydrogen so 
short, the population of the excited states is low enough for photoelectric 
ionizations from these states to be negligible; also, stimulated emission is 
negligible compared with spontaneous emission. 

If all the proton-electron recombinations were to the ground state, there 
would be no reason, in the absence of other cooling mechanisms, for the sphere 
of hot hydrogen not-to extend continually outwards. The dilution of the stellar 
radiation reduces the density of ionizing photons, and so also the proportion of 
free electrons in a steady state, but does not change the high mean energy with 
which the electrons are expelled; the process outlined above might therefore 
be expected to continue indefinitely, the extent of the hot region being limited 
simply by the longer time that diluted radiation takes to heat up the gas, and so 
ultimately by the lifetime of the star. However, in the presence of diluted 
thermal radiation the hydrogen atom fluoresces, emitting fewer high-energy 
photons than it absorbs, the balance being made up by emission in the Balmer, 
Paschen, etc. series (9). ‘Thus there is a net extinction of ionizing photons. 
Strémgren has shown (6) that this leads to a sharp division between a sphere 
of almost completely ionized hydrogen (HII zone) and the outer HI zone of 
almost completely un-ionized hydrogen. Spitzer and Savedoff (5) have performed 
detailed estimates of the temperatures of static HI and HII regions, taking into 
account all the heating and cooling mechanisms occurring; they give 10? deg. K 
as a good average for the HI zones, while HII zones take up temperatures of 
the same order as the exciting stars. 

3.2. Mathematical discussion.—The energy density of radiation of frequency v 
at distance r from the star is given by 


82 hv dv 

c (ehikT 7) 
where W is the distance dilution factor R?/4r? and e~*” is the extinction factor 
in frequency v, defined later. Let Ny, N,, N, be the number densities of 
hydrogen atoms, protons and electrons respectively, and let N,, No, ... N, be 
the densities of H atoms in the ground, 2nd, ... sth quantum state, so that 
Ny==N,. 

From Gaunt’s work on continuous absorption by hydrogen-like atoms (10), 

the absorption cross-section from the sth quantum state is given by 





p,dv= We-*, (30) 


25xime® 1 1 ay 


(31) 
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with the Gaunt factor put equal to unity; m and e are the electronic mass and 
charge respectively, c the velocity of light and A Planck’s constant. The 
application of detailed balancing (9g) leads to the formula 


2h?y? = 2h? 2h"dg i I (32) 
——— CT. = 
mc" mc ov sy 3 


B,’= 


for the recombination probability to state s for unit current of electrons with 
velocity v. The factor s* arises from the degeneracy of the sth energy level. 

The binding energy of the sth level is y,/s®, where x, =2-15 x 10" erg. The 
frequency of the radiation emitted or absorbed is related to the electron velocity 
by 

hv = hmv" + x,. (33) 


From (32), and an assumed Maxwell distribution for the electron gas 


3 
(27mkT,)*” v2 dv e —mv*2kT', wis 


(7, being the electron temperature), we have for the rate per cm* of recombination 
to the sth state: 


n,v)dv= 











_ 8rh®og N.N, 1 I I u.(” et dx ‘ 
R,= 2 (2mm)? y, (RTo) = =(we | , ), (35) 
where 
eS 
° at, 


The total loss of electrons from the gas is then just ©R,. ‘Thus must be balanced, 
in statistical equilibrium, by the rate per unit volume of photo-ionization from the 
ground state: 

p,o,'dv 





Q=Nuc | 


v=zilh hv 


' (36) 


where p, and o,' are given by (30) and (31), and Ny has been written for N,, 
owing to the low occupation of excited states (cf. Section 3.3). Expression (36) 


reduces to 
8a ty ety 
O=Nyz%W | —agemmnre GO, 
c J ylkr * (e* — 1) 


The rate of photoelectric absorption of energy per cm? is given by multiplying 
the integrand in (36) by hv before integrating : 


(37) 





NaS a 2 70WRT | 25° (38) 
Hence the mean kinetic energy of the — electrons is 
ary 
ulkT wn5 fen? Ha 
[2 e-*»dx —Xu (39) 
ulkT x(e* ig 1) 


If the extinction factor e~*” is put equal to unity, this expression reduces te kT, 
Eddington’s result. 
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Similarly, from (32) and (34), the rate of emission of radiant energy by the 
electron gas by recombinations to na and excited states is given by 


8rrh8o, oy 
(2am)8? c2 (RT, aa, NV ( 2 5) (40) 
and the rate of loss of kinetic energy by the gas is given by subtracting from (40) 
the expression 


. 82rh8o, r? e-* dx 

> 9 

3 xRe= Gempta ETN > (Fame os ): (41) 
Hence the mean loss of kinetic energy per initiates is 


(rue JS") 
mem | a pa 
I r® e- * dx . (42) 


x 











If a steady state has been reached, the electron gas will have warmed up the 
heavy particles to its own temperature by Coulomb collisions, so that on balance 
no energy is transferred from the electron gas by collisions. In the absence of 
any other heating or cooling mechanisms, the condition of thermal equilibrium 
is satisfied by equating (39) to (42); anticipating that e~’» differs sensibly from 
unity only at the edge of the HII zone, it is then found, on using tables of the 


integral [ e~” dy/y (x1), that 7, = 7. In fact, there are in general other cooling 


mechanisms active (§); in particular, it is shown in Section 3.3 that collisional 
ionization will prevent the temperature of the gas from rising much above 
10‘ deg. K. Thus the precise temperature of the gas is determined by a much 
more complicated equation of thermal balance, but the results of Spitzer and 
Savedoff (5), and of Section 3.3, show that the assumption of a standard value 
of 10 deg. K for 7, is a good approximation. 
From tables (11), it is found that 
io a rO er 
p> ~ ues | iat 
s=1 § 
where 7, has been put equal to 10 deg. K. As the term for s=1 has the value 
unity very nearly, it follows that three-fifths of the photons absorbed are 
re-emitted in low, non-ionizing frequencies. The degree of ionization y, defined 
by 


= 2°5, (43) 


Us 


N, ” N, = oN, ~ Nn) - YN, (44) 


is given as a function of r by equating the recombination rate = R, to the rate 


&s=1 


of ionization (37), and making use of (43); the result is found to be 
yN _ 2 (2amkT)? (3)" x1 {" e~*» dx 





G9 FS T) RT J nr%(e*—1)’ (45) 


where N is now the density at r of the heavy particles (protons and H-atoms). 
The extinction factor e~”” is given by the equation 


dr, = N(1 —y)(a(v)o,1)(3-08 x 101%) dr, (46) 
where (3:08 x 10'%) has been inserted so that r is measured in parsecs. The 
factor a(v) is the proportion of the absorbed photons of frequency v that are 
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re-emitted in low frequencies. ‘The emission from the gas is isotropic, and not 
just in the outward directions; the use of «(v) in (46) will therefore tend to 
under-rate the extinction slightly, and the correct factor should be between «(v) 
and 1. However, it is easily seen that the error involved is not large, and the 
value ? found from (43) for the integrated extinction factor will be used. 7, as 
defined is the analogue of the optical depth, but the distinction is clearly drawn 
between absorption and re-emission in ionizing frequencies, and re-emission in 
lower frequencies. 

From (37) it is seen that the bulk of the photoelectric absorption occurs 
very near the absorption edge; also, the high-frequency photons emitted by the 
gas are from (35) and (41) nearly all of energy near the ionization energy of 
hydrogen. It is therefore a legitimate approximation (6) to give e~* in (45) its 
value w at the absorption edge, so that (45) becomes 

2N = 2 (2amkT)8? ( T, \¥? kT 
ii a). wn 


We now introduce the variables x and z of Section 2, i.e. 


r=(5-19x 10° °M/0)x, | 


8 
N/N,, =(§x 10-3)z. (48) 
Then (46) and (47) reduce to 
din 
pais BN ..*°(y?x?2?), (49) 
1-y) _., [{2x* 
C -cn. (=), (50) 
where 
logig C = 2 logy) M’ — 2 logy, R’ — logy) T’ + 12°34 + 13°53/7", (51) 
and ‘e 
logip B= 4°595 + logig C + logy, M’. (52) 


The value 7,=10* has been used here. Quantities with primes have been 
expressed in solar units. Equations (49) and (50) differ from Stro6mgren’s through 
the presence of z, measuring the density variation. 

Integrating (49), we have 


w—1= —BN,? | y*x?s* dx, (53) 
~ ro 


where use has been made of the boundary condition that 7, vanishes when x 
has the value x, corresponding to the surface of the star. Hence ionizing radiation 
is completely extinct at a distance x given by 


1= BN,,? ; yx" 3? dx, (54) 
J x, 
where y is given by (50). As Strémgren remarks (6), the value y=1 may be 
used in (54) as y differs sensibly from unity only within a narrow region near x: 
as soon as the proportion of neutral hydrogen becomes important, the absorption 
of photons increases rapidly, and the extinction of the ionizing frequencies is 
soon complete. 
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Substituting z=A/x*y from (12) into (54), and putting y=1, we have an 
equation for (AN,,), proportional to the rate of accretion, in terms of x: 


I = B(AN..)* |” dx/x2y2, (55) 
From Section 2 the allowed values of ) as a function of x are known; (55) therefore 
fixes the density at infinity so that all the ionizing radiation may be extinguished 
within the sphere x. 

3.3. Discussion of the approximations made.—The following assumptions are 
implicit in the theory so far: 

(a) In computing the degree of ionization, the contribution from absorption 
of energy by excited states of hydrogen has been neglected. 

(b) No modification has been made to take account of the motion of the gas. 
If the time taken to reach ionization equilibrium were comparable with the time 
of flow of matter through the HII zone, then the results found from (55) would 
seriously overestimate the required value of N,, for extinction, as the proportion 
of neutral hydrogen at any point would be significantly greater than that given 
by (50). 

(c) The flow is taken as isothermal in both zones; i.e. it is assumed that the 
heat generated by compression is rapidly absorbed and radiated away. In the 
HII zone, collisional ionization becomes active as a cooling mechanism once 
the temperature rises much above 10‘ deg. K. It has to be shown that the 
consequent drop in the proportion of neutral hydrogen does not diminish seriously 
the absorbing power of the gas. 

(d) The diminution of gravity due to the pressure of absorbed radiation 
has been neglected. 

These effects are now discussed in turn. 

(a) The equation of balance for transitions to and from the sth energy level 
for s>1 is 


x 


ea) s 8-1 
R,-Ne| "22 dv=N, = Py + © N,Pry (56) 
r=] 1 


4 Aalh hv 





r=s8 


where P.,, is the spontaneous transition probability from state 7 to state k. (The 
dilution of the radiation justifies the neglect of both line absorption and 
stimulated emission.) From (30) and (31), 





[" p,o,* 820,Ww [" dx 
y= —_——. 
hv css) er ¥(€7—1) 


(57) 

~ dalh 
It is known that the lifetimes of the excited states of hydrogen are of order 10~* sec, 
and so the probabilities P;, are of order 10%/sec. On substituting numerical 
values, it follows that for W appreciably below unity the value of (57) <10%. 
Thus equation (56) reduces effectively to one of balance between photoelectric 
recombination and spontaneous transitions, and in the generalized equation of 
balance for the electron gas 


oe | 4 
x Nie ——dv= XR, 8 
s=1 oe Xglh hv ’ s=1 : (5 ) 


only the first term on the left-hand side makes a significant contribution. 


30 
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(6) Consider an element of matter with its degree of ionization y considerably 
lower than the equilibrium value as given by (50). Then by (37), the rate of 
decrease of (1 — y) is given by 

d 87 , kT 
ges aaro—o Ni — 9). (59) 
From (47) the space-variation of (1—y) in matter in ionization equilibrium 
satisfies 
d(i-y) 2dr 
~ Bas * (60) 
-—" 
Thus matter ionized to a degree y(r) given by (50) will attain the value y(r— dr) 
in a time of order 
a % x1 ar 
pi walk AL =< 
S2o,Ww° kT?’ (61) 
while the time for matter to flow through dr is dr/v, where v is the local value 
of the velocity. Insertion of numerical values shows that for all r 
oa — e-er RT 4 
Xi 2 
so long as 7 is greater than about 10° deg. K. As the present theory is irrelevant 
to such low temperatures, it follows that in all cases of interest deviations from 
ionization equilibrium are negligible. 

Similarly, it can be shown that matter flowing into the HII zone is heated 
up to about 1o* deg. K and ionized rapidly enough for the approximation of 
assuming a sharp discontinuity between the two zones to be adequate. 

(c) Consider first the HII zone. The work done in compressing unit mass 
of the inflowing gas is 


(62) 


> __t.®% 
p?>Dt ip? dr’ (63) 


On substituting (11) of Section 2.1, and using the standard values 10 and 10? 
for the temperatures of the HII and HI zones respectively, (63) reduces to 


1-88 x 108 / y (dz “)) 
a)! mM’ ( re : (64) 


If there are no other cooling mechanisms active, this energy must be absorbed 
by collisional ionization, if the zone is to remain isothermal. The probability per 
second of collisional ionization of a hydrogen atom by an electron is taken (12) as 

_  i2meb kT, 
“12 (2amkT,)'? x, 
where 7) is the temperature of the gas and 4 a factor of order unity. At each 
ionization the kinetic energy of the electron gas is decreased by x, erg; the rate 
of loss of energy per second per gram of gas is therefore 
Nmy 


Some of this energy will be returned to the gas by the reverse process of collisional 
recombination, while the rest will be transformed into both high and low frequency 
radiation by photoelectric recombination. As photons of energy less than x, 
are lost to the gas, collisional ionization is a powerful cooling mechanism. In 





~ 


—, (65) 


=7 x 10%57,!? exp (— 1°57 x 10°/Ty)y(1 —y)N. (66) 








No. 4, 1954 The influence of stellar radiation on the rate of accretion 451 


a gas at rest, the temperature is automatically cut down toa value at which collisions 
are negligible. In the present problem the temperature of the gas at any point 
is given by including in the equation of energy both (64) and (66); as, however, 
the other terms in the energy equation (cf. (39) and (42)) are insensitive to changes 
in Ty, it is adequate to equate (64) to (66). 

The value of y is given by the generalized equation of ionization equilibrium 


ER, +P NS =O+7(1—Y)N*Sp (67) 
where Q and R, are defined in (36) and (35) respectively and S,, is the probability 
per second of collisional recombination. 

Provisionally, both S,. and S,, in (67) are neglected, so that the equation 
for y reduces to (50). The heat of compression is absorbed by collisions when 
1-88 x 108 / y(dz/dx)\ — Sixx, N. Na 
M’ z Nm 
=7 x 10767,1% exp (— 1°57 x 10°/Ty)y(1 —y)N. (68) 
From (68), (35) and (43), y(1—y)N?S,,. / & R* is found to be 
1 


8= 


5°6 x 104 (- =) /wN.. 


which has a maximum value 3:75/M'N.,.. Thus for any value of N,, large enough 
for the accretion rate to be significant, neglect of the S,, term in (67) is a good 
approximation. The S,, term is negligible a fortiori. From the principle of 
detailed balancing (12), S., is given by 


I- 
So = ( “aie ), Sia (69) 


where the suffix 7 refers to conditions of thermodynamical equilibrium. In the 
presence of diluted thermal radiation, the value of y*V/(1 —y) is from (47) very 
nearly W(y?N/(1—y))7, and so collisional recombination is much less probable 
than collisional ionization, which is itself negligible. Hence y is given accurately 
by (47), and values of N,, calculated from (55) will be trustworthy. From (68) 
it is seen that the temperature of the gas is not allowed to rise above about 
1-5 x 104 deg. K, so that the temperature variation through the HII zone is very 
slight. To sum up: collisional ionization is able to preserve nearly isothermal 
conditions without seriously increasing the degree of ionization over the value 
given by (47). 

If radiative processes are not strong enough to prevent some compressional 
heating of the HI zone, the maximum allowed accretion rate will be reduced, 
and the values of N,, as computed from (55) will be increased somewhat. The 
results of this paper therefore refer to the best possible conditions for accretion. 

(d) As emission from the gas is isotropic (stimulated emission being negligible) 
the rate of absorption of momentum per cm* of the gas from the radiative stream 
is from (38) 








810, wre 
~—— WRT we" N(1 —y), (70) 
or from (47) 
Saragk hs «By 5 N23, (71) 
c  (aamk)*? T,1? * 


30% 
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Hence the ratio to gravitation of the force due to radiation absorption is 
Sarayk h* I [ee 

2 

r 





. > aria & N® 
ce} (2mmk)3? T, 12 * 
where z and x are given by (11) of Section 2.1. 

We confine attention to the cases of maximal accretion (sequence (4) solutions 
of Section 4.1). The product zx*=(1-12)/y, and this attains its maximum 


value at x=X; hence for a given HII sphere, neglect of radiation pressure is 
justified for all x <x if 


= 2:14 x 10°°M’'N, 2x?, (72) 


B=2'39 x 10° °M’'N,,/y.<1. (73) 


From (24) and (55) it follows that 8 decreases with increasing N,, until x reaches 
the value ex,; for larger N,, and correspondingly smaller x, 8 increases to infinity 
with JN... 

From Tables III and V, we have the following results: with x= 2-45 x 10°, 
B has the values 0-33, 1-68, 8-98 for Bo, O8, O5 type stars respectively, while 
with x=10-°, the corresponding values are 0:04, 0-25, 1°34. Thus for stars 
earlier than type Bo radiation pressure from photoelectric absorption is not small 
even for moderate values of N,,. The reduction in gravity will lead near the star 
to smaller velocities and hence higher densities than those given by the present 
theory. The consequent increase in the absorptive power of the gas reduces 
the radius of the HII zone for a given value of N.,. For bright enough stars the 
radiation pressure will exceed gravity; we must therefore impose an extra 
condition that the material entering the HII zone has enough momentum to 
reach the star in spite of the deceleration within the zone. It is hoped to extend 


the theory so as to take account of both this effect and also of radiation pressure 
on dust particles in the cloud. 


4. The combination of accretion theory with the theory of radiation heating 


4.1. Numerical results.—On applying the ionization theory to the accretion 
problem it is at once seen that the results obtained depend critically on the 
behaviour of the flow near the star. In solutions which are subsonic near the 
star the velocity decreases to a small value as the surface is approached, and the 
finite accretion rate is maintained by the density becoming large. From (12) 
and (24), y and z behave near the star like 


respectively. Substituting this expression for y in (55), the integral becomes 
so large that only very small values of N,, are required for extinction of ionizing 
radiation to be complete just outside the star’s surface. Hence for “‘ subsonic” 
type solutions the edge of the HII zone will in all cases of interest lie very close 
to the star, so that Bondi’s theory effectively governs the accretion rate. But 
the densities near the star given by (75) are so enormous as to be quite unplausible 
—with x= 10 °, 2 turns out to be 10”, corresponding to a mass-density of order 
107° g. While it is difficult to rule out rigorously these “‘subsonic’”’ solutions, 
it is unlikely that the potential energy of an element of gas in the field of the star 
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should all be ultimately converted by compression into thermal energy (subse- 
quently radiated away) and none into kinetic energy. In the non-spherical 
problem of accretion by a star embedded in a cloud possessing angular momentum, 
it is conceivable that a very high density could be built up near the star by the 
joint action of gravitation and centrifugal force; but this effect could not be 
combined with a large accretion, as the centrifugal force would prevent flow 
into the star. It seems likely that the “‘ subsonic”’ solutions arise from the artificial 
assumption of a steady state, and a rigorous discussion of the time-dependent 
equations of flow would eliminate them by showing that the natural tendency 
of a spherically symmetric system is to approach a state in which matter falls 
freely on to the surface. The “subsonic” type solutions can thus be eliminated 
by an extra boundary condition demanding that the star’s gravitational field 
should lead to rapid inflow, and not build up an enormous pressure gradient. 
In the rest of the paper only “supersonic”’ type solutions will be discussed. 
The following stellar types are considered, with the Sun for comparison 
(6, 15): 
Tasce III 
Spectral type ] R/R T o/T 
Os 30 0063 
O8 12 0'126 
Bo 14 8 o'2 
B2 10 6°: 0°245 
Ao 2°34 ‘ 0°470 
Sun I I 


oe 


\ Xo logi9 B 
19 14°965 
‘765 13°42 
"503 12°57 
‘85 12°82 
‘5 9°61 
38 2°46 


NNN O 


-+ WN 


The problem is then the determination of \,, from (55) for the two sequences 
of “supersonic” type solutions discussed in Section 2. 

(a) Sequence with x>7-5x10-*.—In order to estimate numerically the 
integrals in (55), a curve of 1/x*y® against x is plotted, y being given by (26) with 
the inequality sign replaced by equality. The following integrals are estimated 
numerically : 

‘75X10 0 


dx/x*y? = 4:65 x 107; 
| J 


“ Lo 


‘1 
dx/x*y? = 1:03 x 10°. (76) 


To 


The slight variation of x) with spectral type is of no consequence. If x>o-1, 
then 


dx/x*y? = 8-67 x 10!(x3 — 10-8) (77) 
O1 
to a good approximation; use has here been made of 

—Iny=13:14+Inx?, (78) 
to which (26) reduces for large x. The results are as shown in ‘Table IV with 
computed values of N,, and AN,, for a typical case, a Bz star. 

A remarkable feature of this sequence is the peak in the computed values 
of N,,atx=%. As«xincreases beyond 0-075, the accretion rate drops off sharply ; 
the increase in the volume of the HII sphere is therefore at first more than 
off-set by the drop in the densities near the star, so that the value of N,, necessary 
for the ultra-violet extinction to be complete within x increases with x. After 


x=%, InA is so near its asymptotic value — 13-14 that with increasing x the 
extra volume of the absorbing sphere more than balances the slight drop in 
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density, and so the computed values of N,, decrease with increasing x. For 
large x, the density through the bulk of the H II zone is nearly constant, and the 
theory reduces effectively to Strémgren’s. 


TABLE IV 
7 


r 
| Lo xty® 7 Ne ssa 


4°65 x10° 01138 = 1309 X 107 

o'l 10°3 10° _— Sr 3a8 a2 X16" 
o'2 it Mate 8:17 1°33 X10 
o'5 1°08 x10!” II‘l5 6-69 10 
0667 a°57 x x0" —11°64 7°09 X10 
I 8-68 x 10! 12°15 6:48 x10~ 

118 1°425 x 1017 13°14 4g x<20°° 


2 
8 


Neo 


wn 
~ 


57 
< 107 38 
x 104 14 
< 10° 3°77 
5 x 105 2°4 
x 10° 1°3 
< 107 1°03 x 10~* 


oe 


Huns Ax] 


MV nNNM 
wm 


The accretion rate, proportional to AN,, decreases monotonically with 
increasing x. Even in the most favourable case, x=0-075, AN, is too small 
for accretion to be an important factor in the star’s evolution. It is further shown 
in Section 4.3 that the solutions between x=0-075 and x = % are unstable. 

(b) Sequence of solutions with x < 2-45 x 10-*.—All solutions in this sequence 
have the same parameter InA=0-1138. The solutions for y and z in the HII 
zone vary as x varies, though they all have similar asymptotic behaviour for 
x0. Two values of x have been chosen: the limiting case x= 2-45 x 10% 
and, for comparison, the case x=10-°. The results are summarized in Table V. 


TABLE V 
o-4n we 3 1078 
Spectral (7#*! 2dx NoiW/B N« . A? dx NoVB Nw 


me ao a 7 “Ac x 3 
type Fa xy? (x=2°45 X 107°) 
Ao 2°65 x 10° 1°35 X 107% 21 2°64 x 104 6-16 x 107% 96 
B2 2°65 x 10° 1:9395x<107° 5:25 x 10° 3°15 x 104 "16 X 107% 2 x 10° 
I 
I 


22 = & 
= «xy (s=10~°) 


O8 2°65 x 10° % 3:23 x 108 ‘16 X 10~* 104 


5x10 2:6 X10 
O5 2°65 x 10° ‘35X10 44 X10° 3°82 x 104 waxie* 57x 96" 


The results in the fourth column are the most important, as they give the 
minimum density necessary for a star of given spectral type to be abie to accrete 
at a significant rate. Although somewhat high, the estimates are certainly not 
unreasonable for cool dense clouds (13). For the accretion rate to be high, the 
cloud must be dense, and the HII zone must be limited to a small sphere 
surrounding the star; densities of order 10° particles/cm* perform both functions. 

It is to be noted that the satisfactory results of Table V are due directly to the 
shielding effect of the high densities near the star, caused by the star’s gravitational 
field. This is shown strikingly by computing N,, from (54) with z put equal 
to unity, as assumed by Strémgren (6): the minimum values of N,, that will 
keep x <2-45 x 10°-* are then 10’ times those given in the fourth column of ‘Table V, 
while densities of order 10° particles/cm* give values for x = 100, with correspon- 
dingly negligible accretion rates. Even for a star with the same temperature 
as the Sun, neglect of the density variations near the star leads to a minimum 
value for N,, as high as 8 x 10‘ particles/cm® in order that the HII zone be kept 
small enough for rapid accretion to take place. The density gradient is therefore 
crucial for the theory, and the conditions under which we may expect the 
mechanism to work will require careful discussion. ‘This is postponed until 
Section 4.3. 
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4.2. The heating effect of the accreted material.—The effect of accretion on the 
atmosphere of a star is studied in the Hoyle—Lyttleton—-Bondi theory of the 
solar corona (14). ‘The incoming material is stopped by collisions with the heavy 
particles in the solar atmosphere, and the kinetic energy of inflow is converted 
into heat. A system of convection currents is set up, transporting downwards 
the energy accreted. In a star of surface temperature below about 10‘ deg. K, 
all the ionized hydrogen will become neutral by the time it reaches the photo- 
sphere, the bulk of the emission taking place at temperatures near 10‘ deg. K. 
In a star of surface hotter than 10* deg. K, the photospheric layers consist of 
largely ionized hydrogen, but Hoyle’s results show that there will still be an 
appreciable ultra-violet emission as the convection currents approach the 
photosphere. An upper limit to the effect of this emission on the accretion theory 
can be obtained by assuming in all cases that all the hydrogen becomes neutral 
before finally being swallowed by the star. 

If A is the accretion rate, the number of ultra-violet quanta emitted outwards 
per second per cm? of the star’s surface is 


1/2A 
bal oes 2 
2 ( 5 my )/47R ‘ (79) 


the factor # being roughly the proportion of recombinations to the ground 
state. 
At distance r from the star, the number density of ionizing quanta is 


Fs. 


5 my 47R® cc’ (80) 


where W and w are defined in Section 3.2. The rate of photo-ionization in the 
cloud by the quanta (80) is therefore 

3 
Fe ue (81) 
47R x1 
where the value of the absorption coefficient at the absorption edge has again 
been used. As most of the quanta will be emitted when the matter is near the 
photospheric temperature 7’, the temperature of the liberated electrons will be 
of order 7°. 

This ‘‘ accretion heating’ becomes important if (80) is comparable with the 
number of ionizing quanta in the star’s natural radiation, i.e. if 


Ww(1—y)N 


4 4 
5 my 


2 
A m= 4om*R? EX kT eit”, (82) 


From (13), the maximum accretion rate is given by 
A=5:25 x 105M"N.,. (83) 


‘Then (82) and (83) give, for a star of temperature 7’, the density at which accretion 
heating becomes comparable with natural heating. The results found must be 
compared with those of Table V, which give as a function of JT the minimum 
values of N,, which allow maximal accretion. 

For an Ao star, N,, is found from (82) to be about 10°, which from Table V 
is larger by a factor 50 than the lower limit for maximal accretion. For brighter 
stars the accretion heating is, for the case x = 2-45 x 107%, an even smaller fraction 
of natural heating, so the lower limits given in ‘lable V are trustworthy. If the 
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star is accreting rapidly from a cloud of much greater density than the critical 
value, the calculation of the actual radius of the HII zone from (55) would 
need to be modified by the inclusion of the accretion heating; thus the radius 
of the HII zone would be rather bigger than that given by Table V, but still 
less than the critical value 2-45 x 10-8, so that maximal accretion can occur. 
The important result is that the lower limits from Table V do not require 
modification. 

For stars a little less massive than those of type Ao, accretion heating necessarily 
predominates—the density given by (82) is much less than the lower limit found 
from (55) for maximal accretion. A newequation of ionization equilibrium can be 
calculated from (35) and (81) to replace (45), and the result substituted in (46). 
The temperature of the HII zone will again be nearly the same as the surface 
temperature of the star, so that K,/K, of Section 2.1 will be 27/107. On sub- 
stitution of numerical values the analogue of (55) is 


T'52 


N,, appears linearly, as the accretion rate A is itself proportional to V,,. The value 
of the integral is taken from Tables IV and V (the modification of the theory due to 
the temperature of the H II zone being lower than 10‘ deg. K is negligible). The 
minimum values of N,, are given by (84) with x = 2-45 x 10~, and are all found to 
be quite small; e.g. for the Sun, N,,=1, as compared with o-o11 given by (55). 
Thus for values of N,, high enough for accretion to be significant, (84) determines 
the radius of the HII zone, and the values found are necessarily less than the 
critical radius 2-45 x 107%. 

4.3. Interpretation of results.—Consider a star embedded in a cloud of density 
greater than the lower limit for maximal accretion. ‘The star will then increase 
its mass substantially in a time short compared with its time of exhaustion of energy 
sources; but the increased mass leads to an increase in luminosity and surface 
temperature, and so to a greater emission of ultra-violet light. Hence the radius 
of the H II zone must gradually move outwards. Provided the supply of matter in 
the cloud does not run out, sooner or later the critical radius x = 2-45 x 10-* must 
be reached. Further increase in ultra-violet emission will expand the HII 
zone beyond this limit; but, as shown in Section 2.5, the pressure is then unable 
to drive cool material into the hot zone, so that a steady state cannot be set up. 
Thus the ultra-violet light is not extinguished, but heats up a still larger zone. 

At first, it might appear plausible that a new steady state is reached, with the 
radius of the HII zone x~o-1; for it is known from Table IV that with this value 
of x there exists a steady-state solution belonging to Sequence (a), with the density 
at infinity precisely that corresponding to the limiting solution of Sequence (4). 
The value of A for this solution is about 4 of the maximal value, so that the 
accretion rate would be considerably lower than in the previous state, but we should 
still have the H II zone limited—because of the high densities near the star—to a 
much smaller sphere than that predicted by Stromgren. However, it is easy to 
show that such a state is essentially unstable. A slight increase in ultra-violet 
emission, due either to a perturbation or to a secular increase in mass, immediately 
expands the H II sphere ; but, in striking contrast to the behaviour of Sequence (4) 
solutions, and those of Sequence (a) beyond the peak found for « = , increase in 
x lessens the absorptive power of the HII zone, as shown by the increase of the 
computed values of N,, with increasing x (cf. Table IV). Thus a slight excess of 


I =3°5 x orton No y 2*x2dx. (84) 
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ultra-violet light over the value for equilibrium immediately leads to a state in 
which the gas is still more unable to absorb it, so that the expansion must continue. 
If the cloud is of moderate mass, the large increase in the thermal energy of a large 
fraction of the mass may destroy the gravitational stability of the cloud and scatter 
it, with interesting dynamical consequences. A more massive cloud will be able 
to preserve its stability; the ultimate steady state reached will be the one in 
Sequence (a) with x~118, since a sphere of this radius corresponds to the same 
density at infinity as does the limiting Sequence (4) solution. ‘The accretion 
rate is quite negligible, and the system is effectively as described by Strémgren. 
The next problem concerns the way in which the initial state of accretion is 
reached. It must again be emphasized that the moderate densities given in Table V 
are due to the shielding effect of the star’s gravitational field. A star condensing 
from a dense cloud, and continuing to accrete further matter from its birth, satisfies 
the shielding condition from the start, so that the use of the density function z in 
the integral in (55) is justified. On the other hand, a star which moves into a cloud 
from outside will initially ionize a zone of roughly uniform density. Ifthe cloud is 
massive enough to withstand the heating of part of itself up to 10% deg. K, then an 
initial state of equilibrium will be approached as expansion of the cloud and 
consequent decrease in density lead to equality of pressure. The radius of the 
HII zone will be that given by Strémgren’s theory, with 4x* substituted for 


9222 dx in (55). However, once the expansion has effectively ceased, infall of 


~ Xe 

matter on to the star must commence, and the consequent increase of densities 
near the star must lead to greater ultra-violet absorption and a shrinking of the 
H II zone. 

The ultimate state reached depends critically on whether the quantity NV, \/B 
is less than or greater than the maximum value 0-71 found in Table IV for x=. 
If N,, is so small that N,,\/B<o-71, then the contracting H II zone must reach a 
new stable equilibrium radius x >%, with a resulting accretion rate that is quite 
negligible. ‘The fact that the lower limit to the densities for maximal accretion is 
less than 0-71// B (cf. Tables IV and V) is of no avail, as the system cannot reach 
a state with x <2-45 x 10-3. 

If N,.\/B>o-71, the situation is quite different. The initial value of x, found 
from Strémgren’s equation 

x? 


BN?, eh (85) 


may be greater than or less than 3, depending on the value of N,,; but even if it 
is greater than §, the subsequent contraction of the H II zone cannot be halted at 
asmaller radius, x > 3, as by hypothesis the cloud density is too great. Neglecting 
for a moment the increasing influence of accretion heating, the contraction must 
continue until a Sequence (4) state is reached, so that rapid accretion can take place. 
As an example, consider a cloud of density N,,=10%. ‘Then the critical 

stellar parameter B is, from 
N,V B=071, (86) 


given by log ,»B=7-7. Assuming 7” « M’, L' « M’*?, (51) and (52) give M’ =1-75. 
From (48) and (85) the initial radius of the HII zone is about 0-016 parsec. 
From (86), a B2 star will need to move into a cloud of density greater than 2-75 x 10° 
if it is to be able subsequently to accrete; the corresponding value for an Ao star 
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is N,,=1-1x 104. For stellar types earlier than Bz, it is seen from (73) that 
radiation pressure at the edge of the H II zone defined by (85) is no longer negligible, 
but in any case the densities required to satisfy (86) for such stars are larger than 
any observed. Thus radiation pressure may become important only when the 
system in any case cannot reach a state of rapid accretion. 

It remains to be shown that accretion heating is not important in this context. 
Consider a system in which only accretion heating is operative; then the equation 
analogous to (55) is 


M * dx 


39x 10° Fesg (No) | sae =F (87) 


This reduces to (84) when A is put equal to 1-12. Then for varying x, this equation 
will behave similarly to (55); in particular, the computed values of NV, will possess 
a maximum, which from (77) and (78) is found to be near x=}. Then (87) 
reduces to 

No a =2. (88) 
In order that, with N ,, given by (86), accretion heating should not prevent contrac- 
tion to a Sequence (b) state, by providing a steady state solution with x>4, the 
value of N,, from (86) should be greater than that from (88). By inspection this is 
shown to be the case for the three values of B just discussed. For very much 
smaller stars, V,, from (88) will exceed that computed from (86), but the density 
given by (88) is then far too small for accretion to be significant: any cloud of 
reasonable density must automatically be denser than both the critical values 
from (86) and (88). We may therefore conclude that in all cases of interest, if a 
star of parameter B entering a cloud of density N.,, is to reach a state of maximal 
accretion, the condition N,,.\/B>o0-71 is sufficient (as far as the considerations of 
this paper are concerned). Accretion heating is never able to prevent the H II 
zone from contracting to the required size, although as shown in Section 4.2, it 
may very well dominate over natural heating in the ultimate steady state. 

It must of course be remarked that this paper deals with only one aspect of the 
accretion problem, and the results found must necessarily be combined with those 
found by other workers in the field. (Thus McCrea’s conclusions on the velocity 
with which a star enters the cloud (4) are of crucial importance, as higher-velocity 
stars will never be reduced to rest, and what accretion does take place will be 
described by Hoyle and Lyttleton’s formula (1). Departures from symmetry due 
to the presence of a magnetic field, or angular momentum in the cloud, are likely to 
play an important part in the theory, but discussion of these problems is beyond 
the scope of this paper. 
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Note added in proof.—Since this work was completed, the author has learned that 
Dr E. Schatzman has arrived independently at similar conclusions, to be published in 
the Proceedings of the Second Symposium on Gas Dynamics of Interstellar Clouds, 
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STAR-STREAMING FROM THE STARS OF SPECTRAL TYPE B 


T. R. Tannahill 
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Summary 


‘The proper motions of 1536 stars of spectral types O and B in the Boss 
General Catalogue have been analysed, according to the two-streams hypo- 
thesis, by Eddington’s method. It is found that the phenomenon of star- 
streaming is clearly shown. The galactic coordinates of the vertex of 
star-streaming are (343°°3, —1°°6), and the equatorial coordinates of the 
solar apex (260°-7, +37°°3). The stars in the regions containing the 
Scorpio—Centaurus group exhibit anomalous behaviour. 





1. ‘The B-type stars have always been regarded as an anomalous group as 
regards the phenomenon of star-streaming. Eddington (1), in a review of the 
existing knowledge of star-streaming in 1914, summed up the position as follows: 
‘‘ A more minute investigation of this phenomenon shows that it is complicated 
by differences in the behaviour of stars according to their spectral type. An 
analysis which treats the heterogeneous mass of the stars as a whole without 
any separation of the different types will fail to give a complete insight into 
the phenomenon. But, until a great deal more material is accumulated, this 
interrelation of stream motions and spectral type cannot be worked out very 
satisfactorily. ‘The outstanding feature, however, is that the stars of the Orion 
Type (‘Type B) seem not to share to any appreciable extent in the star-streaming 
tendency. ‘Their individual motions, which are always very small, are nearly 
haphazard, though the apparent motions are, of course, affected by the solar 
motion. ‘They thus form a third system, having the motions of neither of the two 
great streams, but nearly at rest relatively to the mean of the stars.” Later, the 
position as regards the B-type stars is again unequivocally stated: ‘“... the 
stream-motion appears with remarkable suddenness. ‘Throughout Type B, even 
up to B8 and Bog, the two star-streams are unperceived ; but in the next type A, 
the phenomenon is seen in its clearest and most pronounced form.” 

Opinion has not changed since then. ‘There are many reasons why investi- 
gations of the B-type stars from this particular point of view tend to be dis- 
couraged. ‘They are not very numerous, and investigations on star-streaming 
generally require a fair number of stars grouped in limited areas of the sky. ‘Their 
proper motions are generally small; Eddington, indeed, removed them (2) from 
the material of the Boss Preliminary General Catalogue before investigating star- 
streaming with the data of this catalogue, for this very reason. Their radial 
velocities are affected by two principal uncertainties: the high proportion of 
spectroscopic binaries and a large K-term. ‘Their distribution in the sky is far 
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from uniform; apart from their concentration near the galactic equator (which 
is not in itself a disadvantage) the distribution within this zone is very irregular, 
and a considerable number of the stars appear to form small groups or clusters. 

2. The present paper is an attempt to exploit the accurate proper motions of 
the Boss General Catalogue (3) in order to investigate the problem anew. There 
are in this catalogue 1625 O- and B-type stars, of which seven belonging to the 
Pleiades and all those outside the zone + 30° galactic latitude have been omitted 
in this investigation. There remain 1536 for discussion—a number which is 
just sufficient to populate ten galactic regions adequately. Consideration of the 
treatment to be given to the “ clusters’ delayed the start of this investigation 
for some time. In the end, the treatment accorded was the simplest possible— 
the clustering tendency was ignored. 

The material was grouped into ten regions, each consisting of a “‘ square”’ 
of north galactic latitude combined with the diametrically opposite region of 
south galactic latitude. Details of these regions are given in Table I. The 
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coordinates adopted for the centre of a region have been chosen after considera- 
tion of the distribution of stars in the region. ‘The centre is a median position 
with respect to galactic latitude and longitude. Systematic corrections to the 
catalogue proper motions were applied as in a previous paper on the later-type 
stars (4). The position angles of the proper motions were then obtained. It 
was deemed advisable to convert these to galactic position angles, and this was 
done using galactic parallactic angles obtained from charts which were constructed 
from the Lund Tables (5). This procedure avoided any need for applying 
corrections for convergency of meridians. ‘The system of weighting adopted 
was the same as in the previous paper (4), namely, a star is assigned half-weight 
if the probable errors of the centennial proper motion in right ascension and in 
declination both exceed 0”-075, or if either exceeds 0”-085. Counts of the 
numbers of stars, in each region, moving in galactic position angles between 
6—5° and 6+5° (@ being a multiple of 10°) were then made, the figures being 
smoothed by taking the running means of three adjacent sectors. These observed 
frequency curves are indicated by the full lines in Figs. 1 and 2, which also show 
the “‘ theoretical” drift-curves for each region which follow from the analysis 
by Eddington’s “ trial-and-error’ method (2). The constants hV, 6 and N— 
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the projection of the space-velocity of the drift on the tangent-plane at the centre 
of the region, the position angle of the direction to the drift apex, and the weight 
of stars belonging to the drift—are given, for each drift, in Table I. 

3. It is clear from an examination of Table I and the drift-curves of Figs. 1 
and 2 that the proper motions of the B-type stars are not sufficiently represented 
by a single drift. In each region it has been necessary, in order to secure a satis- 
factory fit, to assign a considerable proportion of the stars to a second drift, the 
number in some cases exceeding that allotted to Drift I. In two regions (8 and 9) 
there is a great preponderance of Drift I stars, accompanied by values of AV, 
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Fic. 1.—Frequency curves with galactic position angle, regions 1-5. 


——— observed values, running means of three adjacent sectors. 
--—-- theoretical drift-curves from Eddington’s analysis. 
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which are exceptionally high. The divergences in these regions are so consider- 
able that the regions concerned have been omitted in arriving at the final solutions ; 
but it should be noted that the divergences relate to Drift I, the existence of 
which is not in question. It is probable that in these regions an alien population 
exists whose characteristics approximate in some respects to those of Drift I— 
this would account for the large number of stars assigned to Drift I, and for the 
high value of hV,. Divergences in this region of the sky were not indeed un- 
expected: this is the area containing the troublesome Scorpio- Centaurus group, 
which cannot, evidently, be ignored. 
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Fic. 2.—Frequency curves with galactic position angle, regions 6-10. 
observed values, running means of three adjacent sectors. 
---- theoretical drift-curves from Eddington’s analysis. 
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4. From the analysis given in Table I (omitting regions 8 and g) two equations 
of condition are obtained for each region: 


xsinl — y cos 1 = —AVD sin 0. 
x cos / sin b+y sin / cos b—2z cos b= —AV cos 6. 


Here x, y, x are the linear components of the space-velocity of a drift relative to 
the Sun, referred to the usual axes associated with the galactic equator (5), / and 
b the galactic longitude and latitude of the centre of the region as given in Table I. 
In forming the normal equations equal weight was given to each of the eight 
regions used. ‘The solution of the normal equations gives x, y, z from which the 
galactic coordinates L and B of each drift apex may be found from the equations 


tan L=y/x, tan B= 2/(x? + y?)'2, 
and the space-velocities hW of the drifts from 


hW = (x? + y? +2?)2. 


TABLE I] 
The drift constants 
Drift I Drift II 

—1'165 + 0°126 +0028 + 0:058 
—O'132. o'r16 —0'490 0°054 
—0°375 0:'087 —0°'409 0°04! 

186°°5 5°°6 ajs°°3, 6's 
—aee £9 — 398 4:2 

1°231 0°123 0639 = 0049 

521 526 

+0063 + 0°044 +0'134 + 0-021 
+ 1°194 120 -0'005 0°055 


87°: 


° ‘I 3578 24°°3 
13°°8 


—_ 77°°9 1°°8 


° 
°: 

—0'294  0°143 —o'625 0:067 
2 
o 


The results are given in Table II, with the probable errors. ‘The corresponding 
equatorial coordinates are also given to facilitate comparison with results from 
other types of stars. Table III gives a comparison of the present results with 
those obtained from the galactic zone of type A5—M stars of the previous 
paper (4). 
Tas_e III 
Drift I Drift II 
As-M B As-M B 
A 04 °°1 87°-o 305°°2 357 °8 
D —10°'6 13°°8 —70°°9 99° °9 


Bearing in mind the probable errors involved, the results indicate satisfactory 
accordance. ‘The right ascension of the Drift II apex for the B- -type stars is 
somewhat high—the angular distance between the two Drift II apices is 15°. 

‘The important feature of the solution, however, is not the actual position of 
the Drift II apex, but the fact that it is determinable. ‘The analysis for the 
separate regions, as exhibited in Table I and in Figs. 1 and 2, shows that a single- 
drift hypothesis is insufficient to explain the observed distribution, and gives 
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the constants describing a second drift, in each separate region, necessary to 
produce accordance. The solution given in Table II indicates that the second 
drifts determined from the separate regions can be interpreted as the regional 
symptoms of a drift which is common to all the regions. In other words, the 
B-type stars show similar drift-characteristics to the other types of stars. 

5. From the values of L, B and AW as given in Table II it is possible to calcu- 
late the theoretical values of AV, and @, in each region. ‘These are shown in 
Table IV, together with the residuals hV.—AV and @,.—@. Considering the 
paucity of the data, the residuals are not unduly large, except for region 5, where 
6,—9, is 49°. Reference to Fig. 1 shows the difficulty of determining 6, in this 
region. As expected, large residuals occur in hV, in regions 8 and g (which were 
excluded from the solution). ‘The directions @, in these regions, however, are 
normal—-an indication that, as suggested in Section 3, the abnormal population 
has some characteristics similar to those of Drift I. ‘The Drift II residuals are 
quite satisfactory. It is to be noted that this is true even in regions 8 and 9, 
where the difficulties of determination of #, and AV, are much greater on account 
of the apparent preponderance of Drift I stars. 

The convergence of the directions of the drifts, as determined from each 
region, towards the appropriate apices, is best indicated by calculating the 
angular distance 8 of the apex from the great circles defined by the directions. 
Values of 8 are shown in Table IV, together with the distances A of each region 
from the drift apices. 

TaBLe IV 
Drift I 
Region ‘ Bp, a hV.—AV, 
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1090 0°05 

127°2 o'5I 0:06 
y 143°2 0°38 0°22 
8-0 143°6 0°38 0°02 
12°0 124°8 0°52 o'1s 
9°7 93°2 0°64 0°16 
149° 20° 16°9 56°0 0°53 0°13 
181° II° 8-7 47°8 0°47 -0'03 
222 II° 9°5 57°6 0°54 -0°36 
228-4 - 16 1°6 88-4 0°64 -or1g 
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156° 
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6. ‘The linear components x,, y,, 2, of the motion of Drift I relative to Drift 
I] are shown in Table V, which also contains the galactic longitude /. and the 
galactic latitude 6, of the vertex of star-streaming and the relative space-velocity 
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© of the drifts. 'The most noticeable feature of these results is the unexpectedly 
close agreement of the vertex position obtained here with the values which are 
normally found in such investigations, namely, (343°, —2°).. Considering the 
probable errors involved, the closeness of the agreement must be regarded as 
fortuitous. But the result removes any shadow of doubt that the B-type stars 
exhibit the phenomenon of star-streaming. 

Table V also gives the characteristics of the solar motion relative to the 
B-type stars. The weights of stars in the two drifts (omitting regions 8 and 9) 
are 521 and 526. ‘The ratio of these is so close to unity that equality was assumed 
in the solution for the solar motion. The results are not notably different from 
those obtained in other investigations. 





TABLE V 
The vertex of star-streaming and the solar motion 
Vertex Solar motion 
iy +1°193 + 0°139 ly 28°°7 
Vy —0o'358 0-127 bo +31°°1 
zy —0°034  0°096 hU, 0°757 
ly 343°°3 5°°9 Ay 260° *7 
by —1°6 = 4°"4 Do +37°'3 
Q 1'247 0138 n/N 0-99 
University Observatory, 
Glasgow, W.2 : 
1954 March 29. 
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VELOCITY ELLIPSOID FROM THE CAPE PHOTOGRAPHIC 
PROPER MOTIONS OF 19 839 STARS 
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Summary 


The proper motions of the Cape Catalogue of 20 554 Faint Stars, 1939, 
have been analysed on the two-stream and the ellipsoidal hypotheses. 
A total of 19 839 stars of spectral classes later than B8 has been used in 
the discussion. 

The results agree well with those from the analysis of the proper motions 
of the Cape Zone Catalogue, 1900, and confirm to some extent the variation 
with spectral class determined in that investigation. 

There is a slight decrease in the longitude of the major axis of the velocity 
ellipsoid for the faint stars compared with the brighter stars of the Zone 
Catalogue, 1900. A vertex deviation of some 15° still remains. 





1. The proper motions which have been examined in the present work 
are those of the Catalogue of 20 554 Faint Stars, 1939, published by the Royal 
Observatory, Cape of Good Hope (1). The stars in this catalogue are all in 
the Cape Astrographic Zone and the catalogue is a sequel to an earlier volume (2), 
published in 1936, which contained the proper motions of the stars in the Zone 
Catalogue of 20 843 Stars, 1900 (3). The 1936 and 1939 volumes are hereafter 
referred to as the first and second volumes respectively. 

The photographs used in the determination of the proper motions and the 
star positions were taken in three series—18g2 to 1896, 1897 to 1905 (with a few 
exceptions) and 1923 to 1928, the last series being taken through the glass. The 
first volume contains all the stars brighter than magnitude 9-00 on the C.P.D. 
scale, including 8560 standard stars observed with the 8-inch transit circle for 
the reduction of the plates. ‘The second volume contains all the stars of magni- 
tudes g:1 to g-5 on the C.P.D. scale and also lists 2 432 stars fainter than magnitude 
9°5—1 120 not being in the C.P.D.—which showed total relative proper motions 
exceeding 5”/century. ‘The two catalogues thus contain over 40000 of the 
127 877 C.P.D. stars in the zone, and are complete to at least magnitude 10-2 
on the Cape Photographic Magnitude Scale. 

The reductions from relative to absolute proper motions were made by 
applying three series of corrections. They were :— 

(a) a curvature correction to reduce from rectangular coordinates to R.A. 
and Dec., 

(6) a correction due to the position of the star image on the plate and its 
magnitude, and 

(c) a correction for the systematic motion of the reference stars due to solar 
motion and galactic rotation. For the latter correction Van Rhijn and Bok’s (4) 
results were used, the mean magnitude of the reference stars being 11:2. The 
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measurements and reductions were performed initially under the direction of 
Sir Harold Spencer Jones and later under Dr J. Jackson. 

2. The proper motions of the first volume have been analysed on the two-. 
stream theory by Professor W. M. Smart and Dr T. R. Tannahill (5, 6) and on 
the ellipsoidal theory by Dr J. Jackson (7). ‘Their results clearly demonstrated 
the presence of preferential motion and were in substantial agreement with the 
results of previous investigations of other catalogues. ‘The value determined 
for the declination of the solar apex was, however, rather lower than that normally 
found. 

In this investigation the proper motions of 19 839 stars in the second volume 
have been analysed both on the two-stream theory—by Eddington’s trial and 
error method—and on the ellipsoidal hypothesis by Schwarzschild’s automatic 
method. Stars of spectral type earlier than B8 have been omitted by reason of 
their restricted numbers and extreme galactic concentration, and binary stars 
have been treated as single stars. Variable stars have also been omitted. 
Included in the 19 839 stars are 1735 with unclassified spectra, these stars being 
mainly among the stars fainter than C.P.D. magnitude 9-5. 

3. For the purpose of this investigation the zone was divided into 16 regions 
each of one and a half hour’s width in R.A., and scatter diagrams were prepared. 
Details of the centres of the regions and the number of stars in each region are 
given in Table IV. From the scatter diagrams counts of the numbers of stars 
in successive 10° sectors of position angle were made. In preparing these 
counts the position angles of stars in the two outer half-hours of each region 
were corrected for convergency of the meridians. The counts thus obtained 
were then smoothed by taking running means of three successive sectors, hence 
yielding the final distribution counts. ‘These counts are the subjects of the 
analyses. ‘The methods of analysis are described in references (§) and (7). 

In both methods of analysis the closeness of the agreement between the 
observed and calculated distributions of the stars in position angle depends on 
the smoothness of the observed distribution counts. In the Cape Astrographic 
Zone the smoothness of the corresponding curves depends mainly upon the 
number of stars involved, and consequently upon the galactic latitude. ‘The 
effects of irregularities are reflected in the uncertain or anomalous values 
obtained for the various parameters. In the two-stream analysis the effect of the 
irregularities is thrown almost entirely into Drift I] as an increased uncertainty 
in the determination of its local velocity and direction. It should be noted here 
that as much of the zone is near the Drift II apex, the determination of Drift II 
in any region is always much more uncertain than that of Drift I. 

In the analysis by Schwarzschild’s automatic method, the effect of irregularities 
in the distribution counts is mainly reflected in the determinations of the position 
angles of the major axis of the velocity ellipsoid and of the solar motion, and, to 
a lesser extent, in the derived values of the axis-ratio and projected solar velocity. 
If anything, the automatic method is more sensitive to irregularities than the 
drift analysis, as the five constants determined for each region in the latter method 
are not so interdependent. The irregularities may, however, be so distributed 
as to have no effect at all in the determinations of the parameters by either method. 
There does not appear to be any evidence in the present case to decide in favour 
of one or other of the methods of analysis, as, although in the two-stream analysis 
the investigator is frequently required to decide which, of two or more calculated 
curves, is the best representation of the observed distribution, the evaluation of 
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the ellipsoidal parameters can often be quite largely affected by a single irregularity 
in the observed distribution or a small error in the derived values of the two position 
angles 6, and 6). 

4. In this investigation five groupings of the material have been employed 
and analysed by each method. ‘The first solution was for all 19 839 stars. The 
other groupings were by spectral classes. As it was found in the analysis of the 
first volume that the behaviour of the A-type stars (spectral types B8—-A5) was 
anomalous, they were not treated separately. An analysis of this group of stars 
would also be hampered by their extreme galactic concentration. The three 
groupings used were (1) As5s—Fs, (2) F8-Gs5 and(3) Ko-M. ‘The final grouping 
was of spectral types A5s—M. 

The complete results of the analysis are given in the appendix, ‘Tables V to 
XVII, in which the symbols A@ and AhV refer to the differences between the 
calculated and observed values of the respective position angles and drift velocities 
(ie. AO =6,.,.-A,.3 AAV =AV.,.—AV,,,.). The other symbols are standard, 
e.g. Smart (8). 

Tables V to IX list the details of the drift analyses for each region and 
group. ‘The convergency of the great circles defined by the drift position angles 
is best shown by the angle 8 between the apex and the great circle defined by the 
position angle obtained for the region. It will be noticed that, although the 
values are satisfactorily small for Drift I, there is a greater degree of scatter for 
Drift I]. ‘The degree of convergency is less in each case than was the case for 
the first volume. ‘This is a consequence of the curves in this investigation being 
less smooth than in Smart and Tannahill’s investigation. ‘This may have been 
caused by the effects of accidental errors in the small proper motions. 

Tables X to XIV give the results of the analyses by Schwarzschild’s automatic 
method. In this case the divergences between the calculated and observed 
position angles are generally much smaller and there is little doubt about the 
degree of convergency. 

In the case of the A5—F’5 stars it was found necessary to omit 4 regions, their 
distribution counts being so irregular that no weight could be attached to their 
analyses. 

Tables XV and XVI give the final results of the analyses by the two-stream 
method, and the derived positions of the vertex of star-streaming and the solar 
motion. ‘The results of Smart and Tannahill’s investigation of the first volume 
are given in brackets. ‘Table XVII gives the values of the ellipsoidal constants 
and solar motion, Jackson’s results being given in brackets for comparison. 

Finally in ‘Table XVIII the results of the two methods of analysis are com- 
pared for the vertex and the solar motion. 

5. It can be seen from these tables that the positions of the drift apices 
and velocities of each group differ little from the values obtained from the first 
volume. When the probable errors of the determinations are considered then 
it can be seen that the differences are so small that we may take the sets of results 
as being identical for Drift | and almost so for Drift Il. The increase of hW, 
with advancing spectral type noted by Smart and Tannahill (6) is not so clearly 
evident in this case, but may still be present when the larger probable errors are 
considered. 

In Table I below the drift constants for each volume are compared with 
those obtained by ‘T'annahill (g) from his analysis of the Boss G.C. It will be 
observed that the results are again in good agreement, the main difference being 
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in the higher values of hW, given by the G.C. This, however, was to be antici- 
pated as the G.C. covers the whole sky and not a zone near the Drift II apex. 

In Table II the values for the vertex of star-streaming and the solar motion 
are compared. The deviation of the vertex is clearly shown, the positions of 
the vertex from each analysis being in close agreement. As a result of the higher 
values of hW, given by the G.C., the G.C. values of the relative speeds of the 
drifts (k) are higher than those obtained from the Cape data. 

For the solar apex, a slight decrease in the right ascension and declination of 
the apex with advancing spectral class is noted from each set of results. For the 
right ascension the three sets of data for each group agree well. The low value 
of the declination of the solar apex noted by Smart and Tannahill and by Jackson 
is still evident. It must, however, be recognized that the Cape declinations 
will have a probable error of about 5° (as shown by the ellipsoidal results). 

An outstanding feature of the table is the high values of N,/N, obtained in 
this investigation. The smallest value obtained in this investigation is 1-4 for 
the Ko-M group, whereas the other investigations average 1-05. ‘This is probably 
a result of the uncertainties in the determinations of Drift II. 

In Table III the constants of the velocity ellipsoid are compared with Jackson’s 
results for the first volume and with Delhaye’s (10) analysis of the G.C. It should 
be noted, however, that the grouping used by Delhaye is rather different from 
that employed in this paper. His groups are B8—As5, Fo-F9, Ko—K2, and were 
further divided into regions of high or low galactic latitude. The figures given 
above are thus approximate means of his results. 

There is some slight indication of a decrease in the longitude of the major axis 
for the faint stars of the second volume compared with the brighter stars of the 
first volume. This is not certain as it is not fully evident in the results from the 
drift analysis. There is no doubt that the deviation of the vertex from the 
centre of the galaxy remains. 

6. In conclusion, the following points are to be noted: 

(1) The results for the analysis of the proper motions of the second volume 
are in complete agreement with those obtained by Smart and Tannahill and 
Jackson from the first volume. The uncertainties are somewhat greater, possibly 
as a result of the effects of accidental errors on the small proper motions. 

(2) The increase in hW, with advancing spectral type noted by Smart and 
Tannahill from the first volume is not fully confirmed. It may well have been 
masked by the larger probable errors involved. 

(3) A decrease in the right ascension and declination of the solar apex with 
advancing spectral type is indicated. 

(4) Abnormally high values of N,/N, are obtained. ‘These may be spurious, 
being probably due to uncertainties in the determination of Drift IT. 

(5) Aslight decrease in the longitude of the major axis of the velocity ellipsoid 
for the faint stars compared to the bright stars of the first volume is suggested by 
the results. 


Finally I wish to tender my thanks to the Department of Scientific and 
Industrial Research for the award of a maintenance grant during the tenure of 
which the investigation was performed. 

University Observatory, 

Glasgow, W.2: 

1954 August 14. 
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APPENDIX 


TABLE IV 
Centres of Regions and Numbers of Stars 
Number of Stars 

Region R.A. All As5-F5 F8-Gs5 Kc-M As5-M 
m 

© 45 22 79 291 
15 394 38 84 244 
45 525 3 128 361 
15 630 164 523 
45 1175 : 293 841 
15 2137 35 q 345 1082 
9 45 2437 5 431 1380 
II 15 1604 381 1174 
I2 45 1526 407 1166 
14 15 1519 E 375 1116 
15 45 2308 q 434 428 1291 
17 15 2419 c 441 464 1258 
18 45 955 315 265 731 
20 15 627 264 196 549 
21 45 600 264 207 532 
23 15 561 5 231 158 439 


Totals 19 839 5168 4405 12978 
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TABLE V 
Analysis of the Drift Curves—All Stars 
Drift I Drift II 
A, Aé, AhV , AV, 9, Ne A» 


fa 


II 0°65 210 216 
‘32 0°90 215 #1294 
12 060 175 2714 
‘04 0:25 180 3764 
04 O30 230 5§:748 
0:20 210 8 782% 

19 045 210 8355 
14 0°50 170 §29 
0°35 160 5194 

‘04 06 86150 395 
‘19 O°5 170 = 752% 
O29 405 «6155 = 7034 
—o-09 ©3. 200 357 
0°31 03 180 276% 
0°03 06 215 293 
° 05 200 314 
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125° 
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119° 
107° 
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TaBLe VI 
Analysis of the Drift Curves—Sp. Type As—F5 
Drift I Drift I 
A, AO, AAV, 6, " A; AG, AhVs, 


> 
_ 


- 23°2 —o'25 
—0'03 
0°28 


36°3 179 0°23 
36°3 7% —Os7 
45°3 4°8 0°02 
59°0 0°25 
74° —o'18 
go o-'18 
104 0°03 
115 -0°08 
122 O'57 
122° —0'07 
115 0°42 
0°33 
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TABLe VII 
F8-G5 
Drift I Drift IT 
ny A, ‘ 


ad 
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Co 


895 75° 
73;  66- 
943 46° 
38° 
3871 
46'9 
66°8 
75°8 
go"! 
105°8 
117°5 
124°4 
124°4 
1754 117°5 
1524 105°8 
129 go'l 
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Tas_e VIII 
Ko-M 
Drift I Drift II 
A, Aé, AhV, A> AhV . 


~ 
> 
_— 
= 
” ts 


gr 


76-9 
61°5 
47°7 
38°6 
37°9 
46:2 
59°5 
748 
90°4 
1050 
117°0 
124'2 
124°8 
118-9 
106°8 
92°4 


0°45 
0°32 
‘12 
24 
05 

‘10 
31 
04 
‘21 
“16 
27 

— 0°35 
—0'36 
o'41 
O'13 
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85 
574 
40 
20 
335 
310 
290 
2724 
2573 
2424 
221} 
195 
175 
155 
115 
1024 


fad 


eSCxnaukwWN 


Drift I 


ny 


1184 
157} 
224) 
263% 
467 
649 
8464 
772 
777% 
709 
796 
773 
3774 
338% 
280 
2314 


kih 
0°570 
0581 
0°654 
0-708 
0683 
0°703 
0°599 
O'551 
0656 
0667 
0668 
0°807 
0816 
°°799 
0649 
0598 


ki/h 


°°745 
0°700 
o'751 
0°532 
0631 
0624 
0667 
°°555 
0882 
0-909 
0816 
0°944 


A, 


74°9 
59°5 
45°9 
37°2 
37°5 
46°4 
60-2 
756 
gor! 
105 °4 
117° 
123°7 
123°5 
116°6 
104°8 
900 


TABLE IX 
As-M 


AhV, 


0°20 
0°35 
0°34 
0°02 
o'l2 
—Or15 
-0'24 
—o'09 
— 0°05 
—o'1o 
—0°30 
—0'29 
—0°30 
0°20 
o'll 
0°25 


TABLE X 


hV, 


O'4 
o'8 
o'7 
o'4 
o'4 
O's 
o'4 
o°5 
O'5 
o'5 
o%3 
o'5 
03 
03 
06 
o'5 


UP 


200 
220 
190 
185 
220 
210 
190 
160 
180 
175 
180 
185 
205 
240 
220 
220 


Ellipsoid Analysis—All Stars 


9% 


67°1 
48-6 
22°1 
12'0 
3446 
324°4 
306 °3 
290°1 
276°0 
265:0 
237°! 
204°4 
170°8 
132°2 
g2°2 
83°9 


9 


Ig'l 

O's 
3189 
306-2 
291°3 
275°0 
281°3 
240°6 
197°2 
177°4 
135°3 
105°5 


Sh 
0536 
0°523 
0382 
0°314 
0°449 
0805 
1-022 
1:087 
1-072 
1-116 
1090 
1°154 
0°957 
0660 
0-699 
0680 


TABLE XI 
Ellipsoid Analysis—A5—F'5 


Sh 
0°385 
0°490 
0°740 
1:098 
1°15 
0-956 
I‘O15 
1°094 
1-465 
0-909 
0°595 
0°843 


9, 


292°7 
263°6 
257°4 
208°1 
140°4 
118-4 
100°8 
86-0 
71-9 
54°4 
35°° 
12°8 
357°7 
343°0 
321°3 
304°7 


9, 


178°9 
132°'0 
I19g'2 
98°8 
82:8 
69°7 
Lg 
29°8 
10°3 
349°9 
329°3 
306°7 


Vol. 114 
Drift II 

No A. Aé, 
172% 46°73 °'9 

86% 51°71 —23°9 
136% 548 1°5 
2594 57:2 - OF 
374 578 —40°4 
433. 569 —36°7 
533% 543 —22°2 
402 504 3°2 
3884 45°55 —I9'I 
407 40°9 —Irer 
495 352 —II'S 
485 31°79 —I10°7 
353% 306 —24°2 
210 3271 —49°5 
252 36:0 22°9 
207% 41°0 —20°3 
AG, Aé, 

2°7 —I11°'0 

46 2°5 

I1°9 —11'4 
— 04 2°5 

o°7 15°6 

o'8 — 08 
—— —~ 47 
— 04 — 58 
— 30 — 63 
— 77 - 9 
— 19 — 26 
— 33 — 06 
—13°4 ~ 7-0 
— 8-3 —12°9 

9°8 — g2 

1°4 — 84 
AO, AO, 

2°5 18°4 
<9 8-9 

14°4 — o9 

6°3 — 12 

3°7 — 14 

471 -~ $3 
— 183 — 54 

29 20 
+ 18-2 — 18 
— 22 — 4°6 

2:0 — 52 

61 — 1'0 





AhV 


0°05 
~0°32 
-O'1g 
o'12 
O13 
0°02 
ovll 
—o-'02 
—o'06 
—o-09 
0:06 
—O°17 
0:02 
0:03 
—0°23 
—o-09 
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TaBLe XII 
Ellipsoid Analysis—F8-G5 
8 Sh 6, AO, 


fae 


65:0 0°529 290°6 6°8 
50-4 0°442 256:2 4°6 
19°7 0-292 244°5 
8-9 0°250 220°0 4°2 
346°1 0°388 140°6 2° 
325°0 0-518 105°0 °6 
312°1 0852 95'0 6-1 
294°8 I 000 80-9 — §°5 
277°1 1084 69°4 3°4 
264°1 0874 50°4 7°3 
234°1 1°307 36-0 °o'9 
204°1 1317 iy 18 
161°3 1:046 08 - 09 
150°0 0°745 341°3 —23°2 
88-8 0°731 319°7 
88-2 0671 310°8 — 08 


TABLE XIII 
Ellipsoid Analysis—Kc-—M 
A Sh 6, AO, 


oCOn auf wWN & 


67°5 0°43I 320°'0 O'5 
53°9 0°444 248°2 — 2°2 
13°5 0°444 261-0 
177 0°322 213°2 — 66 
334°2 0483 150°0 
327°5 °°949 122°9 -~ a9 
302°5 0-891 103°3 6:0 
289-2 0932 88-4 2‘1 
276°7 1086 748 °'7 
261°8 0962 61-4 — 24 
238-2 1:064 390 o'8 
222°5 1°116 17°7 
176°8 0°736 3°9 
124°0 0641 348-3 
104°0 0°605 326°4 
86-7 0°493 3081 
TABLE XIV 

Ellipsoid Analysis—A5-—M 

A Sh A, 


I 
3 
4 
5 
6 


onl 


69°5 0°443 294°3 
49°'0 O'511 242°0 
22°4 0°317 250°8 
13°6 324 205°0 
347°4 "504 137°5 
323°6 "727 I115°0 
307°6 ‘968 99°3 
292°8 036 85-0 
2768 ‘062 710 
266°5 ‘O41 51°1 
237°7 ‘073 35°0 
204°0 1165 13°4 
167°7 0888 359°2 
138°0 0°693 341°2 
93°8 0689 320°4 
86-7 o'591 308-1 
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All Stars 
XxX, —0°013+0°053 
(0°052 0°035) 
y; 1°486 0°053 
(1°446 0°032) 
Z —O'IgO (0°054 
(—0'217 0°044) 
A, go"5 20 
(B8°-o = =1°*4) 
dD, — one: a9 
( 8°-5 o°°6) 
AW, 1°498 0°053 
(1°463 0°035) 
N, 12 513 
(10 745) 
All Stars 
X, 0°037 +0°067 
(0-093 07029) 
Y. —o'108 0:067 
(—o'140 0°026) 
Z: —o'603 0068 
(—0-°628 0:036) 
Ay 2890 33°°3 
(303 °6 9°°5) 
Dz 79°O0 7"! 
(— 7g°O 3°°$) 
hw, 0614 0068 
(0o°650 0°035) 
Ve 7326 
(7578) 
Soin. All Stars 
A, 271° G42 °3 
(a71°-5 = 1° 6) 
D, “a Ses 
(—14°°5 2°:0) 
G 343 "4 
(343 °2) 
g, ol 
(+0°:5) 
k 1648+ 0041 
(1-638 0-043) 
A, 269° -6 
(265 -1) 
D, 25°°4 
(+ 26°-1) 
hU, 0800 
(0879) 
Ni/N¢ 1-708 
(1°4T0) 
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TaBLe XV 
Drift Constants 
Drift I 
As5-F5 F8-Gs5 Ko-M 
—0'000+ 0°072 0001 +0°074 —0'036+0°'078 
(—0-028 0°054) (—o'oI9 0°042) (—0'029 0°052) 
1°548 0°074 1°449 0°074 1°369 0°078 
(1-677 0-041) (17533 07038) (17442 0°047) 
—0'296 0-092 —0'226 0076 —0o'212 0-080 
(—0:'220 0:067) (—0°136 0°053) (—o-099 0°064) 
goo 2°6 go""o  3°°0 ms 2s 
(g1°-o 1°°8) (90°-7_ _1°°6) (ors 2°-1) 
~10°8 4°5 — 89 2°4 — 8°-8 3°°3 
(— 7s 29 (— 5°11 2°) —- yee 
1°576 0°074 1°615 0°068 1°385 0°078 
(1-692 0-042) (1°539 0°038) (17446 0046) 
2052} 30784 25684 
(1760) (2015) (2191) 
Drift I 
A5-F5 F8—G5 Ko-M 
0°081 + 0°083 0°056-+-0°070 —0°'047 -|-0°079 
(0-182 0°033) (0°046 0028) (0040 0°037) 
0'064 0086 —o'189 0°070 —o'188 0:°079 
(—o'141 0°030) (—0°139 0°025) (—0'086 0°033) 
—0'583 0°107 —0°666 0071 —o'585 o-o81 
(—0°619 0:040) (—0°670 0°034) (—0°587 0046) 
321°°7 40 286°°5 20°°5 256°:0 23°°4 
(ga8°:2 «-9°"9) (288°-4 10°9) (295°°2 21°°2) 
80°-0 19°:2 —- 74°" §°°9 — 71°6 7°4 
(— 6°43 20 (—23F7 2s). (— Bs . 3°'1) 
0592 0°106 0°694 0°023 0616 0-084 
(0°661 0°039) (0686 0°034) (0595 0°046) 
11525 20893 18363 
(1503) (1941) (2030) 
Smart and Tannahill’s results are given in brackets. 
TaBLe XVI 
The Vertex of Star-Streaming and the Solar Motion 
As-F5 F8-G5 Ko-M 
272°'9+3°°9 271°-9+ 3°°6 269°-6+ 4°71 
(270°-6 2°-0) aya: 1°45) (272°°6 2°-4) 
— ws 2 —rg° 0 3°°s ct ie ae 
(—12°°3 2°°5) (—18°'2 2°) (—17°-7 -3°°0) 
347 “8 343 °° 343 2 
(344°°5) (340 -2) (340°°6) 
< ol 2°°6 
(+2°°8) (—1°-2) (—1°-4) 
1°639+ 0°044 1-697 + 0°075 1‘601 + 0084 
(1:873 0°052) (1°756 0-049) (1-606 0-060) 
268° -3 268° -6 273 ‘2 
(265 -3) (268 -9) (269 -5) 
22°'4 25 °5 27°83 
(+25°°6) (-+-29°-2) (+25°°3) 
1047 0'972 0808 
(0°934) (0o°817) (0-780) 
1°781 1°473 1°398 
(1-173) (1-038) (1-080) 


s results are given in brackets. 


Vol. 114 
As-M 
0-013 + 0°068 
(—0'022 0°045) 
1°433 0068 
(1°538 0-041) 
—0o'246 0-069 
(—o0'144 0°056) 
So°°*s 2°97 
(o's 2°97) 
ee Fe 
[3:3 
1°454 0°067 
(1°545 07041) 
7781 
(6002) 
As5-M 
0°024--0°058 
(o°092 0°024) 
—o'145 0°058 
(—0O'I3I 0°022) 
—o'604 0058 
(—0'649 0°030) 
279 "5 22 "4 
(30479 8-4) 
ws -§°3 
6s 2°36) 
0°622 0058 
(0°668 0-024) 
5197 
(5446) 
A5-M 
270 ‘4+3°°2 
(a73°°S 1°8) 
—12°8 3°-4 
(—16°:8 2°-2) 
344 °2 
(342° :o) 
ss 
(—2° 3) 
1°618+0°100 
(1-747 0°048) 
268° -8 
(267°-5) 
25 9 
(+ 27°-3) 
o'SgI 
(0-838) 
1°497 
(1-102) 
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TaBLe XVII 
Ellipsoidal Constants and Solar Motion 
As-F5 F8-—Gs5 
278°-7+ 4°°3 270°°7+ 3°°6 
(276° -6) (272° 6) 
15°44 5°°2 ~15°°5$3°°5 
(—23°9) (— 16°-4) 
345 °°8 342°°0 
(348°:1) (342°°1) 
—6°0 +o°6 
0676 + 0'072 0°554+0°075 
(0°552) (0°56) 
266° -9+ 3°°5 271°:2+4°"4 
(266°-o) (269°-9) 
19°84 4°1 30° 9+ 3°'8 
(+ 23°:2) (+ 31° 6) 
1'056+0°064 0'950+0°063 
(105) (0°89) 


Ko-M As5s-M 
269° 6+ 2°-6 271° -6+2°°6 
(268°-5) 
—19°-7+2°-6 16°-2+2°°5 
(-— 16°-2) 
3378 341°°8 
(340°°3) 
~0°-4 —0°3 
0-606 + 0:046 0-601 + 0°046 
(0°65) 
275 “4$5°'5 269°:5+4°°6 
(272°*3) 
24°24 5°°2 29°°s+4°°2 
(+31°°8) 
0°831+0°074 o'819+ 0:065 
(0°82) 


Jackson’s results are given in brackets. 


TaBLe XVIII 


Comparison of Ellipsotd and Two-Stream Positions of Vertex and Solar Motion 


Soln. All Stars 
A 271°-6+2°:0 
(270°:2) 
D —17 ‘2+2 ‘0 
(—13°-4) 
G 340°°'0 
(343°°6) 
g +o°'8 
(+2°-1) 
K/H 0602 + 0:036 
(0-625) 
Ay 271°°5+4°:2 
(266° -9) 
D, +25°-0+3°'9 
(+ 28°-3) 
AU, 0°948 + 0064 
(0°94) 
G 
Sol. TS. Ell. 
All Stars 343°°4 340°°0 
P.E.~+2 
As5-F5 3478 = 345°°8 
ae 
F8-G5 343°°0 342°-°0 
+4° 
Ko-M 343°°2  337°°8 
1 y 
As-M 344°:2 341°'8 
ee 


g 
T.S Ell. 
o°'l oe 
rs an? « 
—o°'r °° 
2°6 —o° 
2°°0 —o”: 


Nw 


Ay 
T.S. Ell. 

8 269°6 271 
+4 
o 268°:3 266 
+3 
6 268°-6 271° 
+4 

‘4 273°°2 275 
+5 

3  260°8 269 
+4" 


Auuet PHUS 


Dy hU, 
T.S. Ell T.S. Ell. 
+25°°4 25°°0 0800 0-948 
+ 3°°9 +o-c6 
22°°4 19 ‘8 1:°047 1-056 
> -<3 +006 
25 °5 +3°°9 0°972 0°950 
+ 3°6 + 0:06 
27°°'0 24°'2 0808 o-S31 
+52 + 0°07 
25 °5 29°°5 o'SgrI o'sig 
+ 42 + 0°07 








REMARKS ON THE REVISION OF FK3 AND ITS RELATION TO N3o 
A. Kopff 


(Received 1954 September 29) 


Summary 


In the first part of the paper the principles adopted for the revision of FK3 
are explained. About 57 catalogues with epochs between 1918 and 1950 
are to be used, including both those with independent instrumental systems 
and those tied by differential observations to the system of a fundamental 
catalogue. For the improvement of the position of each individual star in 
the system of FK3 all the catalogues are used; the method of correcting 
position and proper motion is that given by Boss in the P.G.C. For the 
correction of the system FK3 itself only the independent catalogues can be 
used. 

The second part gives a comparison between FK3 and N30, shown by 
the tables at the end of the paper. Using the system of the G.C. as a basis 
of comparison, it is found that there are some systematic differences in the 
N3o system between those stars which are in the FK3 catalogue and those 
which are not, particularly in two zones of declination. The differences are 
shown in Tables I-IV. A satisfactory explanation of this difference cannot 
be given. 





The revision of FK3 

1. The many new star catalogues which have been published since 1930 
have made it desirable to revise the catalogue FK3. The mean epoch of the 
Auwers stars in the FK3 is about 1900; that of the additional stars is somewhat 
later. The need to undertake this revision was confirmed in 1952 by the I.A.U. 
in Rome; moreover, for various reasons, an increase in the number of stars in 
the fundamental system is required. Recommendation No. 18 of Commission 
4 is concerned with this matter. In the meantime the desire for an increase in 
the number of fundamental stars—which was particularly emphasized by 
geodesists—has been met by the compilation of a supplementary catalogue 
(FK3 Supp) for which positions and proper motions have been given (1). 

2. In response to a circular letter of 1952 December (2), various observa- 
tories have forwarded observations not so far published or have promised to 
forward them within a very short time. About 57 newly observed catalogues 
are available for the revision of FK3. The revision was commenced in the 
middle of 1953. 

3. The first survey of the expected correction of FK3 results from a com- 
parison with the catalogue of Standard Stars by Morgan, called N30 (3), 
about which details are given later. For the part of the sky observed by the 
observatories in the northern hemisphere it is noticed that individual corrections 
of the single stars within the system of FK3 are small. But in the southern 
hemisphere the number of catalogues available for the compilation of FK3 was 
rather small, so that now the individual differences are appreciable. 

4. The main task is, however, the systematic correction of FK3. The 
N30 catalogue shows some large systematic differences from FK3 ; but it is 
not possible to obtain valid information about the necessary correction of FK3 








No. 4, 1954 Remarks on the revision of FK3 and its relation to N30 479 


by comparing FK3 with N30. Though the catalogue N30 is to a great part 
free from systematic errors at its mean epoch 1930, this does not apply to the 
catalogue positions for the epoch 1950 (compare p. 480). These positions are 
derived from the 1930 positions by the use of proper motions which are based 
only on the systematically improved positions of the G.C. at the time of its mean 
epoch. The proper motions of the G.C. were completely left out of account ; 
this has already been mentioned previously (4). 

The systematic revision of the G.C. at its mean epoch, however, has not been 
carried through adequately. The corrections in R.A. of the G.C. were obtained 
without regard to the declination and are given in the form of a Fourier series 
(Cat. N30, p. xvi). According to a letter from Morgan, corrections of the 
right ascensions, depending on Dec., were taken into consideration only south 
of —30°. For the declinations there are small deviations from the Fourier 
series north of +30° and south of —30°. In contrast, when compiling FK3 
the relations of the catalogues were determined by graphical methods, owing 
to the very complicated and irregular nature of the catalogue differences. Thus 
the change of Ax, with Dec. was completely taken into account after subtracting 
the derived Aa; (5). In the same way the values Ad, also showed a dependence 
on Dec., which had to be applied (6). This procedure is being used again in 
the present revision. 

Moreover, when FK3 was compiled, none of the catalogues which were 
observed before 1845 were used, because they contain many systematic errors. 
On the other hand, the proper motions of G.C. were considerably influenced 
by the catalogues of these epochs before 1845. 

5. For the revision of FK3, the catalogue N30 cannot be used at all, as 
strictly speaking it is not a fundamental catalogue on account of its special 
composition and character. This was previously pointed out in 1950 at the 
Conference on Astronomical Constants in Paris (7). 

There are particularly two parts of the sky where it is necessary to carry out 
systematic corrections of FK3. One is in the extreme northern sky. When 
compiling FK3 the discussions of Kahrstedt (8) in particular showed that in 
the north polar region there were strong marked discordances between the 
individual instrumental systems ; in this region the system of FK3 was largely 
based on the instrumental system of Pulkovo. Secondly the southern part of 
the FK3 system (south of — 30°) needs a more extensive systematic revision. 
The number of newly observed catalogues which can be used is, however, very 
small, and is even smaller than those used when FK3 itself was compiled. It 
is therefore most important to increase the number of meridian circles in the 
southern hemisphere. 

6. The revision of FK3 will be carried out, as was the case with the com- 
pilation of FK3, in two steps: it begins with the determination of the individual 
corrections within the system FK3 ; then the correction of the system itself will 
be deduced separately. The discussion will, in the first place, be concerned 
with the corrections depending on the declinations ; the corrected residuals will 
then be discussed for the dependence on R.A. Only those catalogues which 
have been observed on an independent system can be used for the systematic 
correction of FK3. The other catalogues are tied in a differential way to an 
existing fundamental system and they can only be used for the individual cor- 
rection ; it would not be logical to relate such observed catalogues to a common 
preliminary corrected system and to combine them with the absolute catalogues. 
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‘The method to be used for the revision of FK3 differs from the procedure 
which was used when FK3 was compiled. According to a method which was 
already proposed by Boss (g), it is intended to correct the positions and the 
proper motions of the fundamental catalogue within its own system. ‘The 
mean errors given in FK3 for the position at the mean epoch and for the proper 
motion of each star provide the basis of the method. ‘This has already been 
pointed out at the colloquium in Paris (7). This procedure was also used by 
Nowacki to correct the position and proper motion of 61 Cygni (10). The 
correction of the system itself has then to follow. 

For uniformity the equinox of G.C. and N30 is adopted as the equinox of 
the revised FK3, though the old value could be used in every respect (compare 
FK3, part I, Appendix, p. 106). 

The revision of FK3, commenced in the middle of 1953, should be completed 
in 1958. Immediately afterwards the final formation of FK3 Supp will follow. 

The further necessary extension of FK3 is mentioned in the last part of 
Recommendation 18 ; the catalogue FK3 is to be supplemented by a funda- 
mental catalogue of faint stars. It should thus be possible on the one hand to 
fix in a precise manner the system of the fainter stars of the old catalogues ; 
and on the other hand to provide a reliable basis for the establishment of a true 
fundamental system with the help of observations of faint extragalactic nebulae. 


The relation between FK3 and N30 


7. In the following a comparison between FK3 and N30, carried through by 
Nowacki, is given in the same form as for the reduction of G.C. to FK3 (11). 
A comparison between the positions and proper motions of FK3 and the standard 
catalogue N30 has already been made by Morgan (3, p. xx1) and in more detail 
by Nowacki (4). ‘These two comparisons differ considerably in the method of 
their derivation. Morgan gives the difference between the catalogue positions 
and proper motions according to R.A. and Dec. in one table. Nowacki 
first determined the systematic differences between the catalogues depending 
on Dec. (Aas, Aus, Ads, Av’;), following the method formerly used by the 
Astronomisches Rechen-Institut. After correcting for these differences, the 
residuals were analysed according to R.A. (Aw,, Ayu,, Ad,, Au’,). The depen- 
dence of these variations on the declination become obvious by this method. 
‘The results of the comparisons are given at the end of this paper. 

8. The other investigations, carried out by Nowacki concerning the relations 
between FK3 and N30, are only mentioned briefly ; a more elaborate publication 
will appear in the Verdéffentlichungen des Astronomischen Rechen-Instituts =u 
Heidelberg. In these discussions it was necessary to include the G.C. in its 
unrevised form. It should be noted that the Standard Catalogue N3o is given 
for the equinox of 1950, and is therefore influenced to a very great extent by the 
adopted proper motions which are based on the G.C. positions. 

The relation (N30—G.C.) was derived anew as a difference between the 
relation (FK3—N 30) obtained in the investigation (4) and the (FK3—G.C.) 
formerly derived (12). ‘These values (N30o—G.C.), derived as a control, are 
based only on the stars of FK3. They are compared with (N30—G.C.),, 
derived from the tables of Morgan (3, p. xIx and 13), which were based on all 
N30 stars listedin G.C. These two comparisons should agree. ‘The differences 
in R.A. and Dec., arranged according to Dec., i.e. 


A=[(FK3—G.C.) —(FK3—N30)] minus (N30—G.C.),, (1) 
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are given in Fig. 1 only for the proper motions. It can be seen from it that on 
the whole the agreement is satisfactory. 

g. There are, however, both in R.A. and Dec. regions with large differences. 
For the R.A. there are the regions at + 50° to +70° and — 40° to —60°. Inthese 
parts the differences formed for the FK3 stars alone diverge appreciably from 
the differences formed for all N30 stars. This was the reason for a more special 
investigation of these areas. The following data are confined to the positions 
and proper motions in R.A. within the two areas mentioned above. 











jrcaa tare 2 











4 4. n walls £ shin 4 iu a L 4 L 1 
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+60 +60 +40 +20° ‘o -20 -60 -80 
Fic. 1.—Differences Su, and Au’, (compare formula (1), p. 480). 





From the individual values Ax and Ay, listed for each star in the catalogue 
N3o, the mean values for the two areas have been formed anew first for the 
FK3 stars and then for the other stars of catalogue N30. The means were 
calculated in two ways: at first for individual Dec. zones of 3° breadth for all 
R.A.’s for both areas, and secondly in R.A. for two hours each extending over 
all declinations of these areas. In both cases overlapping means of the 
consecutive values were formed. Finally the differences of these mean values 
for the FK3 stars and the non-FK3 stars were formed. ‘The results are shown 
in TablesI-IV. ‘The differences in R.A. and the proper motions in R.A. between 
the two groups of stars, which can be seen quite clearly from Fig. 1, are given in 
these four tables. 

The tables give the overlapping means Ax, Ayu (centennial) and their differ- 
ences ; the mean epochs and the number x of stars of these means are also given. 
For the last column of each table, values 


Ap = 


were calculated at first for FK3 stars and non-FK3 stars separately ; these 
values are not given here separately, but the last column Ax.) in Tables I-TV 
contains the differences of the values. 
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Taste I 


Means of (N30—G.C.) for %s9 and yw (centennial) in the area +50° to +70° with Dec. as argument 
(compare for 4x39) formula (2), p. 481) 
Unit for 4a and Ap: 08-001 





FK3 stars non-FK3 stars FK3 stars minus 
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TABLE II 


Means of (N30—G.C.) for %9 and y (centennial) in the area — 40° to —60° with Dec. as argument 
(compare for 4a,39) formula (2), p. 481) 
Unit for 4« and Sy: 08-001 





non-FK3 stars FK3 stars minus 
non-FK3 stars 





mean epoch 
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TaB_e III 
Means of (N30—G.C.) for «59 and j (centennial) in the area +-50° to +-70° with R.A. as argument 
(compare for Aa39) formula (2), p. 481) 
Unit for 4x and A: 08-001 
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10. At present it is not possible to give sufficient reason for the fact that, in 
both areas mentioned above, the FK3 stars systematically differ from the other 
stars of catalogue N30. It must, however, be mentioned that the FK3 stars 
are to a large extent identical with the standard stars of Lewis Boss (14). So it is 
possible that in the uncorrected G.C. the standard stars differ systematically 
from the other stars. This conjecture is supported by the following circum- 
stance. For the A« the differences in Tables I-IV for the epoch 1950 are 
much greater than for epoch 1930. The figures for 1950 have been influenced 
by proper motions, but those for 1930 have not; the proper motions, however, 
are based on the positions of G.C., so that systematic differences between different 
groups of stars in G.C. will influence the catalogue N30 for the epoch of 1950 
through the proper motions adopted in N30. 

It is also possible that the systematic deviation of FK3 stars from the other 
stars is the result of the additional correction applied to G.C. for 1900, before 
the derivation of the proper motions of N30. The corrections applied to the 
positions of G.C. for 1900 were valid only for the mean epoch of 1g00. In 
deducing the proper motions of the individual stars, the exact difference of 
epochs for each star was used ; but as the epochs of the standard stars, which 
in general correspond to the FK3 stars, deviate to some extent from the other 
epochs, it is quite possible that the difference found in N30 has been caused by 
the insufficient additional systematic correction of G.C. 


32* 
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TaB_e IV 


Means of (N30—G.C.) for %59 and (centennial) in the area —40° to —60° with R.A. as argument 
(compare for 42,39) formula (2), p. 481) 
Unit for 4a and Ap: 08-001 

















] 
FK3 stars non-FK3 stars FK3 stars minus 
non-FK3 stars 
a mean epoch mean epoch 
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h 1900 + |1930+ 1900+ |1930+ 

0-0 +55 t+ 97 | +69] +1°5] 15 +-62 +115 +8-0 | +2°3] 30 — 7] —-18 | — 5 
1-0 +39 + 58 | +5-4] +1°5] 15 +58 +108 +81 2°8| 32 —19 | —50 | —12 
2-0 + 39 + 571 +5°5| +1°8] 13 + OI +116 | +7°9 2-83] 31 —22 |] —59 | —12 
3:0 + 31 + 53 +7°7| +1°8] 14 +51 + 98 | +7°2] +2:-9] 26 —20 | —45 | —13 
4°°}]) +21 37 bee PSST Fo + 36 e Fe Eevee. ce 31 a ia = 
A tT 23 + SOE TIO Tess 4 +20 F-45 } $S-3 1 +3°2} 92 3 + 8 I 
6-0 +2 + 62 +-6°8 | +0°7] 14 +13 + 40 +79 2-8] 30 45 +22 6 
7:0 +15 + 50 +5°8 | +0°4] 13 +19 + 64 +8°7 2-6] 30 — 4 —14 — 3 
8-o +13 + 48 | +66] +0°3] 15 +160 + 60 | +8-6}] +2-6] 31 - 3 —12 — 2 
9°0 r1o r 4! +5°7| +02] 14 II 49 | +7°O| +2°7| 34 =. 3 = ise 
10-0 — 2 + 14 | +66] +0°6] 13 +22 + 72 +7°9 | +2°7] 33 —24 —58 —15 
11-o +10 | + 47 +8-4] +2-1] 12 +34 |] +104 +9°2] +2°7] 31 —24 | —57 -14 
12-0 +2 + 85 | +7°7]| +1°5| 13 +41 +127 | +8-9] +2-9] 30 —17 | —42 —II 
13°0 +26 + 94 | +7°5| +0°3] 14 +45 +139 |] +8-1] +3:°1] 29 —1I9 | —45 = 
14°0 +- 2° + 83 | +7°5| +1°0] 15 + 44 142 | +8-9] +2-°8] 26 —21 | —54 | —13 
15:0 +20 + 77 +7°5| +20] 15 + 37 +119 +3-5]) +2°8] 25 —17 —42 —tII 
16-0 +42 | +123 +7°7| +1°2] 14 + 38 +115 | +8-6] +3:-1] 30 + 4 + 8 I 
17°0 +42 +112 +O°5 0-0} 13 +50 tI4t | +9°4] +3°0] 34 - 8 —29 -— 6 
18-0 || +41 +102 70 o'r} 13 |} +58 | +149 | +88] +2°5] 27 17 | —47 | r1 
190 t 54 r 123 7°3 1°7 | 13 £58 | +133 | +8°5] +3°1 | 27 ~*~ ese 2 
20°0 +50 t+ 96 -6°3 | +2-0] 14 +56 +1ig | +8-7] +3°5] 25 — 6] —2: — 3 
21-0 +45 + 74 +7°3 1-4} 16 + OI +115 +7°2 | +3°4] 20 —16 —4I -10 
22:0 t- 50 + O5 6-9 Iv! 13 +71 +134 | +7°4] +3°3] 28 -15 — 39 —10 
23°0 |] +74 | +137 | +0°5 1-0} 12 PIS tf Piss P77 1 + ST 33 ri- th=—2 
24°0 +55 + OF L 6-9 1-5 | 15 +62 fF 115 + 8-o 2-8] 30 — 7 —18 — 5 
Meanj| -+ 32 76 6-9 1-2 43 | +103 8-2 2-9 = oe > 
























































A further possible explanation for the difference between the FK3 stars and 
the others may be found in the method of formation of N30. ‘Tables I-IV 
show that the values Az;,. and Aa,s9) are of the same general character, both 
according to Dec. and to R.A. But the values are smaller for Ax(g9) than for 
Aa». As the compilation of N30 was based on differences from the G.C., 
systematic differences which might be present in G.C. between the standard and 
the other stars could therefore reappear in the catalogue N30. 

It is possible that these different causes were working together when N30 
was compiled At present it can only be pointed out that N30 does not possess 
a fully homogenous system, especially in the two regions mentioned above. 

11. In the following Tables V-XI the reduction of the N30 system to 
FK3 is given in the sense (FK3—N3o), in the same manner as the reduction 
(FK3—G.C.) was formerly given (11). 

The tables contain the values Ax,, Ad, and then the Ax, Ad, for zones of 
30° breadth for the epochs 1g00 and 1930 to 1970. ‘The values of Ax, in 
the north and south polar caps are multiplied by cos5; the values Az, 
are on the contrary not reduced to the equator. For smaller zones in Dec. the 
values of Nowacki for Ax,, AS, and the proper motions must be used (4). Correc- 
tions for magnitude equation are not given, because the values are absolutely 
insignificant. 
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TABLE V 
FK3—N30 
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TABLE VI 
FK3—N30 (++ 90° to - 60°) 
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TaBLe [X 
FK3—N30 (0° to — 30°) 
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TABLE X 
FK3—N30 (—30° to —60 ) 
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TABLE XI 
FK3— N30 (—60° to go ) 
4a, cos 6 Unit : 08-001 4b, Unit : 0”-o1 
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PHOTOGRAPHIC OBSERVATIONS OF MONOCHROMATIC 
MAGNITUDE AT SIX WAVE-LENGTHS 


R. v. d. R. Woolley, S. C. B. Gascoigne and A. de Vaucouleurs 


(Received 1954 March 9) 


Summary 


Relative monochromatic magnitudes for 85 bright southern stars have 
been determined by photographic photometry at the following six wave- 
lengths: AA 4050, 4250, 4550, 4945, 5980, 6360. The relative energy 
distribution is examined as a function of spectral type. Comparisons are 
made with the results of other observers and it appears that while the standard 
of accuracy achieved for colours is good, there is evidence for a possible 
error of about 2 per cent in the magnitude scale. 





1. Introduction and observational.—Observations of monochromatic magni- 
tude, by Hogg and Hall (1), and of colour gradient, by Gascoigne (2), have been 
made on Mt Stromlo in recent years and have already been described. Both 
these sets of observations were relative. It had been planned for some time to 
make them absolute, by fixing their zero-points with respect to an accurately 
calibrated terrestrial source. As a convenient way of doing this, the slitless 
spectrograph with which both programmes (1) and (2) had been carried out was 
removed from the Reynolds reflector and mounted separately on the polar axis 
belonging to the Farnham refractor. Here it became effectively an objective 
prism telescope of four inches aperture, and could be directed equally well at a 
star or at a standard lamp set up in the Observatory grounds. Sorne observations 
of relative magnitude and gradient, originally of a preliminary nature, were made 
with this equipment. As these showed features of interest, this part of the 
programme was expanded to include all the southern stars brighter than about 
the third magnitude, and forms the subject matter of the present paper. The 
comparisons with the standard lamp will be described in a sequel. 

A smaller angle prism was substituted for the original 45° prism in the spectro- 
graph. This gave a dispersion of 290 A/mm at Hy, adequate for the purpose in 
hand, and which, with the plates used (mostly Kodak O-F), enabled a spectrum 
of a third magnitude star to be obtained in 60 seconds. Spectra were widened 
$ mm by rocking a diagonal mirror as described in (2). The calibration of each 
plate was effected by spinning in turn in front of the prism a number of sectors, 
concentric with the circular telescope aperture. A set of these was made in } 
magnitude steps from o to 4 magnitudes. 

Dr Hogg kindly checked the sector attenuations, by attaching to the spectro- 
graph the photomultiplier he had used in (1), and measuring with it the apparent 
brightness of the standard lamp (a suitable steady light source) as seen through 
each of the sectors in turn. ‘The values determined photoelectrically agreed in 
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all cases with those computed geometrically to within a few thousandths of a 
magnitude. 

To calibrate a plate at least one and often two bright stars were given a con- 
siderable range of sectors. Three exposures were usually made on remaining 
stars, if necessary through different sectors to enable suitable densities to be 
recorded at each wave-length. Each plate accommodated about 40 spectra, 
and on a good night four or more plates could be secured. Atmospheric absorp- 
tion was determined each night, by observing stars as they rose or set through a 
considerable range of air mass. In all, 130 plates were secured on 66 nights, 
and of these, 110 were measured. The remaining plates were omitted, usually 
because of unsatisfactory absorption on the night on which they were taken. 

Mean measurements of absorption are shown in Table I, which includes 
also the standard deviations on individual nights and earlier determinations. 
The agreement between the present values and those determined previously on 
Mount Stromlo (1, g) is noteworthy. 

TABLE | 
Determinations of absorption (unit : 0-01 mag/sec 2) 
A4050 A4250 A4550 A4Q45 AS5Q80 Ab6360 
Mean absorption 37 32 26 21 15 12°6 
Standard deviation +43 +40 +34 +36 2°8 2°4 
Hogg and Hall (1) 25°4 (east) 
28-1 (west) 
Sun and Sirius (9) A4510 A4g45 As5550 Ab6209 
26 21 17 13 


Plates were measured in the projection photometer described in (2), with a 
slit of projected width 0-08 mm, and at the following wave-lengths: AA 4050, 
4250, 4550, 4945, 5980, 6360. ‘The bandwidth accepted thus ranged from 
about 20 angstroms in the blue to about 50 in the red. ‘The regions were chosen 
to avoid absorption lines as far as possible, but this cannot of course be done 
completely, and our magnitudes refer to the integrated spectrum and not the 
continuum “ between the lines”. 4550 was the wave-length at which Hogg 
and Hall made their measurements. It was hoped also to have a set of magnitudes 
at about A 5450, the photovisual effective wave-length, but the plates used had 
either rapidly varying or low sensitivity in this region, and it was therefore not 
suitable for measurement. 

2. The reductions.—The observations were reduced by selecting as standards 
the thirteen stars in ‘Table II (6 Velr and « Virg were occasionally used as second- 
ary standards). ‘These were intercompared frequently and least-squares solutions 
made of their relative magnitudes at each of the six wave-lengths. The zero- 
point of 74559 Was chosen so as to make the mean magnitudes of the standard 
stars agree with Hogg and Hall’s magnitudes; the zeros of the other five magni- 
tudes were then fixed so as to make the mean gradient agree with Gascoigne’s 
gradients. Details are shownin Table III. The zero of these gradients should, 
like Gascoigne’s, be close to that of the Greenwich gradients (3). A comparison 
of the twelve stars in common shows that we average about 0-028 bluer than 
Greenwich. If we omit four stars in Orion for which there is some evidence 
that the Greenwich gradients are anomalously red, this figure is reduced to 0-015. 
Similarly a comparison of the seventy stars in common with Gascoigne shows 
that we average 0-015 bluer. While one would like to see this figure reduced, it 
can nevertheless be regarded as satisfactorily small. 
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TasBce II 
The standard stars 

my mHH mp-MyH Gy GC Gy-Gy 
« Erid 0:296 0'274 0°022 —0°225 —0'236 O'oll 
y Orio 1°429 1°434 — ‘005 —0'272 —o'285 — 013 
e C Maj 1°307 1°295 ‘O12 —0°'277 —0o'278 ‘O01 
B Cari 1°599 1°596 003 0024 O°OI7 007 
« Leon 1°224 1°238 - ‘O14 —0'092 -o'118 "026 
y Cent 2°023 2°038 — O15 —0'022 —o'o16 006 
B Cent o*410 0°402 008 —0°327 —0°332 "005 
x Ophi 2087 2096 “009 0°290 0°305 ‘O15 
e Sgtr 1°730 1°742 - O12 ool! 0°036 —- O25 
x Aqil 0820 0°834 ‘O14 0°410 0°416 — +006 
« Pavo 1°731 1°718 O13 —0'239 —0'249 ‘O10 
a Grus 1°587 1°585 “002 —o-196 0'180 ‘016 
a PscA I°117 1°108 "009 O'lig 0°107 ‘O12 


my and my, are the magnitudes at A 4550 determined on the Farnham refractor and 
by Hogg and Hall (1) respectively. Gy and Gp are gradients determined on the 
Farnham and on the Reynolds reflector (2). 


3. Results.—The results for the 85 stars * on the programme are set out in 
Table III. These consist of the six magnitudes, together with the gradients 
computed from them. Spectral types were taken from Yerkes lists, or from 
some unpublished work of Miss Woods. ‘There was some difficulty in deciding 
on a suitable measure of colour for the later type stars. As is well known, 
magnitude differences between these and early type stars do not plot linearly 
against 1/A, the later types exhibiting a deficiency in the violet. However, over the 
limited range A 4550~—A 6360 the relation appears reasonably linear (see Section 4), 
and gradients for stars later than F5 were therefore computed from these four 
wave-lengths only. 

From internal evidence the accuracy at each wave-length for the standard 
stars is better than o-o1 mag. For the programme stars it averages 0-02 or 0:03 
mag. Red stars tend to be less accurate, possibly because of their rapid variation 
of photographic density with wave-length. 

4. Relative energy distribution.—The material in Table III enables us to 
examine the variation of relative energy distribution as a function of spectral 
type. Of previous work in this field (4, 5), ours is most nearly like that of Kienle, 
Strass] and Wempe (6). ‘To reduce observational error we represented each 
spectral type by the means of the following groups of stars: 

Bo : « C Maj B: II, 8 Cent Br II, « Orio Og V, « Orio Bo I. 
Bs : « Pavo B3 IV, « Erid Bs IV, « Grus Bs V, « Leon B8 V. 
Ao : « Sgtr Bg IV, « Pegs Bg V, 8 Cari Ao III, « Cent Ao III. 
As : « PscA A3 V, « Ophi As III, « Aqil A7 V. 

Fo I : « Cari Fo la, « Leps Fo Ib, « Cari Fo Ib. 

Fo V : « Hydi Fo V, « Cire Fo V, 8 Tr Au F2 V. 

Fg 1: 6 C Maj F8 Ia, 8 Agar Go Ib. 

G1 V : » Boot Go IV, 8 Hydi Gi IV, « Cent G2 V. 

G5 Ill : 8 Corv G5 I, w Velr Gs III, 8 Ceti G8 III. 

Ko III : 8 Gemi Ko II, ¢ Cent Ko III-IV, « Scor Ka III. 





* The following bright stars were omitted : @ Cari, which is in a cluster; y Velr, a double star 
only partly resolved on our plates; « Crux, omitted inadvertently. 
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Star 


y Pegs 
B Hydi 
B Ceti 
y Erid 
y Hydi 


B Erid 
B Orio 
y Orio 
B Leps 


6 Orio 


Leps 

Orio 

Orio 
. Colm 
f Orio 


© Orio 
C Maj 

a Cari 

. C Maj 
C Maj 


> C Maj 
C Maj 
C Min 
Gemi 
Pupp 


Pupp 
Cari 

> Velr 
Cari 
Cari 


¢ Velr 
Hyda 
Leon 
Velr 
Cent 


8 Leon 
Cent 
Cruc 
Corv 
Cruc 


R. 


h 
00 
fete) 
fol) 
oI 
ol 


95 
05 
05 
°5 
05 


95 
95 
95 
95 


Photographic observations of monochromatic magnitude 


A. 


m 
II 
23 
41 
36 
57 


05 
12 
22 
26 


Spec. 


B2 V 
G1 IV 
G8 III 
Bs IV 
Fo V 


A3 V 
B8 Ia 
B2 IV 
G2 Il 


Og-5 III 


Fo Ib 
Og V 
Bo I 
B8 Ve 


Og°5 III 


Bo II 
Br II 
Fo la 
Al V 
Br Il 


F8 Ia 
Bs I 
Fs5 IV 
Ko III 
Os 


F6 II 
Ko, B 
Ao V 
Ao III 
Fo Ib 


B2 IV 
K3 III 
B8 V 

Gs Ill 
Bo III 


A3 V 
B2 Ve 
B2 1V 
B8-5 III 
M3 II 


TABLE III 


Collected results 
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Tasie IIT (cont.) 


Star R.A. Spec. A4050 A4250 A4sso A4g45 5980 A6360 G 
h m 
B Corv 12 32 G 3°32 360 3°14 2°77 2°27 ‘10 
Musc 12 34 ] 2°32 2°36 2°41 2°43 2°62 “56 
Cent I2 39 1°97. 201 2°02 = 2°05 2°02 ‘Ol 
Cruc 12 45 O91 O94 ~~ & 1°05 I°rl 1°22 "24 
Virg 13 23 068 073 076 ~~ 0°85 0°94 95 


RDS R 


Cent ce ; 19g! 1°93 1°99 2°14 2°24 24 
Cent 13 52 a%4 290 a6 8 2°37 2°49 "44 
Boot 13464 G f 3°32 3°20 2°89 2°73 2°36 IS 
Cent 14 0°32 ©'37 O41 o'si o'61 0°63 
Cent 14 Ko III-IV 3:39 3°27 2°65 2°21 1°63 1°49 


DBD st uv a 


Boot 14 K2 p 1°73 1-48 0°76 —o-48 0°70 
Cent 14 G2V 0'74 0°56 0°13 060 0°67 
Cire 14 Fo V Sr 3°27 3°18 3°12 2°83 2°89 
Lupi 14 B2 Il 2°02 2°06 2°06 . 2°28 2°24 
Libr 14 A3 met 2°78 2°75 2°72 ; 2°57 2°53 


ReReReR 


c 
= 


B Lupi 14 B2 V 2°35 240 2°37 2%5 2°53 
y Tr Au 15 Ao V Eu 2°72 2°74 2°74 : 269 2°64 
B Libr 15 B8 V 2°38 4240 82°39 . 2°49 2°42 
y Lupi 15 B3 V 2°48 2°47 2°51 . 2°79 2°66 
« CorB 15 Ao V 2°02 , 2°07 . 2°09 2°06 


Tr Au 15 F2 V 3°13 2°98 . 2°62 
Scor 15 / 2°60 ’ 2°64 2°86 
Scor 15 2°10 2°19 2°24 
Scor 16 2°98 2°89 2°73 
Scor 16 2°47 : 2°61 . 2°83 


N= Oo 


= 


p 
5 
o 
T 


COAN OF 


N 


NNN NN 


Scor 16 fer _— 2°31 ‘ 0°33 
Ophi : ‘ 2°51 ° 2°43 
Tr Au eae a 2°90 1°30 
Scor : ‘ 2°97 1°80 
Scor . 2°45 2°64 


cmaRueR 


Arae : 2°76 2°96 
Scor 1°38 ° 1°56 
Ophi ° 2°09 
Scor : 2°05 . 1°58 
Scor 2 1V 2°19 


Sgtr 1°73 
Segtr 1°89 
Segtr ‘ 2°59 
Aqil ] 0°82 
Pavo 22 ‘65 1°73 


Agar 21 ; | 3°35 
Capr 2 2°98 
Grus . 1°59 
Psc A ; . 1°12 
Pegs 2 : E 2°40 
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From the observed magnitudes we subtract those computed from the mean 
gradient for each group, on the assumption of a linear relation between m, and 
1/A. For later types the gradient was computed from the four longest wave- 
lengths only. Residuals in the sense O-C are shown in Table IV and also plotted 
in Fig. 1. The change in the character of the residuals for types later than.F5 
is well shown. It will be noticed that with these later types departures from 
linearity are more marked for giants than for dwarfs. 


—_—_—_—______ 
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44050 4250 4550 4945 5980 6360 
Wave-length 
Fic. 1.—Relative energy distribution, The vertical spacers are 0-4 mag. 
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5. Colour comparisons.—Most comparisons so far made between gradient 
and colour indices have shown a good deal of scatter (3). This may be attributed 
on the one hand to observational errors, and on the other to effects of absolute 
magnitude, interstellar reddening, unrecognized variables, binaries of different 
spectral types and peculiar objects such as Wolf-Rayet stars. If however the 
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comparison is confined to nearby, i.e. unreddened, main sequence stars, there 
should be a strictly functional (though not necessarily either linear or one-to-one) 
relation between two measures of colour, even if they refer to widely different 
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Fic. 2.—Gradients plotted against Kron’s infra-red photoelectric colours. 
Open circles are supergiants, crosses space-reddened stars, filled circles main sequence. 
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G 
Fic. 3.—Gradients plotted against Cape photoelectric colour indices Cyp. 
regions of the spectrum. ‘The same considerations should apply to comparisons 
between monochromatic and heterochromatic magnitudes. 
In Fig. 2 our gradients are plotted against infra-red photoelectric colour 


indices measured by Kron and his collaborators (7), and in Fig. 3 against 
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“* blue-yellow ” photoelectric colours measured at the Cape (8). In these figures 
high-luminosity stars are distinguished by open circles, space-reddened stars by 
crosses and main sequence stars by filled circles. The correlations are good. 
In fact mean errors of +0-02 in the gradients, and of +0-o1 in the infra-red 
colours, are large enough to account for almost all the scatter. Exceptions are 
a Erid and 8 Hydi. A similar situation exists in the comparison with the Cape. 
Here « C Maj, « Circ and perhaps 8 Hydi stand off on one side, and 8 and y Tr Au 
on the other. « C Maj and y Tr Au, with Cape colours of —o-19 and —o-20, 
and gradients of —o-o5 and 0-08 respectively, represent an extreme case of 
scatter. (The vertical scale is much more open in Fig. 2 than in Fig. 3.) The 
anomalous infra-red colour for « Erid (see Fig. 2) is puzzling, as it has high 
weight on all systems. A possible explanation could be that it has avery late-type 
companion; but no variations in light or radial velocity have been reported for 
this star, and it seems unlikely that it is a binary. 

There are 70 stars in the present programme for which gradients have also 
been determined on the Reynolds (2). The arithmetic mean difference between 
the two determinations is 0-04, with the following large residuals, all in the sense 
Farnham-—Reynolds: « C Maj +0-11, « Cari +0-10, 8 Hydi +0-19, 8 Libr —o-14, 
« Orio —o-14. Where such differences exist the evidence is usually in favour of 
the current series. For instance, it had been pointed out in (2) that « C Maj and 
a Cari were both too blue for their spectral types; these anomalies are now 
largely removed. And the present gradients for x Orio and £ Libr agree much 
better with their photoelectric colours. 

We have twelve stars in common with Greenwich (3). ‘These include four 
in Orion, which we make about 0-05 bluer inthe mean. Most modern observers 
agree in including these among the bluest stars in the sky. The mean difference 
for the other eight stars is about 0-02, which is small. 

6. Magnitude comparisons.—The indications from the preceding discussion 
are that our gradients have reached a very satisfactory order of accuracy. When, 
however, Our 4559 are plotted against the Hogg and Hall magnitudes a scale 
difference of about 2 per cent becomes evident. ‘This is investigated in Fig. 4, 
where we plot as a function of magnitude 

(a) R— Msg —0°14 G. 

(6) myn — Ma550- 

(c) BPg—C,, — magas- 
Here the R are Kron’s red magnitudes, and BPg and C,,, colours and magnitudes 
taken from the Cape Mimeograms (8). Before heterochromatic magnitudes 
can be compared with monochromatic magnitudes in this way colour effects 
must be eliminated as far as possible. We do this by interpolating between the 
monochromatic magnitudes to the effective wave-length of the heterochromatic 
system (as was done by Hall, reference 4, p. 234). This turned out to be neces- 
sary only in case (a), as the effective wave-length of BPg—C,, is very near A 4945. 
In case (a) we have worked from misogq rather than from mgg¢g9, which was nearer 
the effective wave-length of R as there was less scatter with the shorter wave- 
length. 

Fig. 4 indicates a well-defined scale error of about 2 per cent, in the sense that 
we have made the bright stars too faint. This difference between our photo- 
graphic magnitude scale and the photoelectric ones (BPg is at least partly photo- 
electric) is disappointing; the weight of evidence undoubtedly supports the 


33 
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latter, although we have been unable to detect any source of systematic error in 
the present series. It is perhaps worth noting that the mean difference 
My — M4559 1S ONly 0-02. 





T 











4 
-2 





Magnitudes 
Fic. 4.—Comparisons between the Farnham magnitudes and those measured by Kron (a), 
Hogg and Hail (b), and at the Cape (c), to test for possible scale error. Vertical spacers are 
0°2 magnitudes apart. 


A new photographic programme on the Farnham telescope with some slight 
modifications of the present equipment will be put into operation in 1954. It is 
hoped that this will throw some light on the origin of the scale error. 


Commonwealth Observatory, 
Canberra : 
1954 March 2. 
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Addendum.—The observations mentioned in Section 6 of the paper have now 
been completed and reduced. The programme consisted of two bright stars 
(a C Ma and « Car), two faint stars (8 Leo and y TrA) and a number of standard 
stars. The results are shown in Table V. The magnitudes at 1 4550A were 
found by adopting Hogg and Hall’s values for the standard stars, i.e. all except 
the four stars already mentioned. For example, the 1954 value given for 8 Cen 
was found by assuming HH values for all the other standard stars with which it was 
compared, and so on. Some stars were used a few times as standards (e.g. 
y Ori and 5 Vel) but have not been kept in Table V as the 1954 observations of 
the stars (considered as programme stars and not as standards) is of too low weight. 

Table V shows the magnitudes found in 1954, together with the weight, 
the magnitude found in 1952, and the photoelectric magnitude according to 
Hogg and Hall. The next two columns show differences. ‘The last three 
columns show the gradient found in 1954, the gradient found in 1952 and the 
differences. The 1954 gradients were found by least squares from magnitudes 
in six colours using the values shown in Table III as standards. 


TABLE V 
Results obtained in 1954 





Magnitude at A4550A Differences Gradients G 
in 

| hundredths 
(1) | (2) (3) ofa 

Photo- Photo- Photo- | magnitude 

graphic | Wt-| graphic | electric |W ___ 
1954 | 1952 Hogg & | (1)- | (1)- 

Hall | (2) | (3) 








Difference 








—o-62 35 —o'59 —o-69 ~ 03 
‘47 27 ‘43 1°49 : 1 
2 17 ‘31 1°30 = - —OI 
57 | 34 ‘60 160 = - +o° 03 
“19 33 ‘22 1°24 — —03 


708 | 20 "04 2°10 : +07 
‘75 | 68 ‘76 0°77 . ‘eT- — 04 
‘41 96 “41 0°40 + - - +02 
‘81 30 "74 2°80 + 

‘so | 27 09 2°10 t + + 02 


B72 | 5° 1°73 1°74 eT = + O04 
0°87 | 70 082 0°83 +5 t | } +02 
rtz | $3 I‘12 I‘tl —I +o: —03 
Root mean squares, unweighted +3°3+3°0 + 0°033 






































The differences photographic 1954 minus photoelectric (Hogg and Hall) 
show a little correlation with magnitude, as is shown by the following least 
squares results : 

Including all stars = +0™”-016—(0-014 + 0°007 S.€.)ig550. 

Excluding « Carinae = +0™-006—(0-007 + 0°:008 s.€. )mg550. 
The differences also show a correlation with gradient, as follows: 

Including all stars = —o™-004 +(0™.062 + 0”-023)Gyg50. 


m™ 


Excluding « Carinae = —o"”-005 +(0”-043 + 0™-030)Gj955. 
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Since « Carinae is both bright and (relative to the B stars) red, its large differ- 
ence 1954 minus photoelectric supports either a scale effect or a colour effect. 
The general run of the 1952 differences suggests very strongly that the colour 
effect is not real, and the authors feel that the scale effect is not real either. In 
our view it is the Earth’s atmosphere rather than photometric technique which 
limits the accuracy of magnitude comparisons over wide angles. 

The agreement between the gradients determined in 1054 and 1952 is highly 
satisfactory. 

Mr Gottlieb and Dr Heintz have assisted us by taking many of the 1954 
plates. 


ERRATA 
M.N., 113, No. 6: 
S. Chandrasekhar, Problems of stability in hydrodynamics and hydromagnetics 
(George Darwin Lecture). 


P. 672, the legend to Fig. 1 should read as follows :-— 


Fic. 1.—The variation of the critical Rayleigh number, R,, for the onset of convection as a 
function of T(=4Q# d* cos* $/v*) for the three cases: (a) both bounding surfaces free, (b) both 
bounding surfaces rigid and (c) one bounding surface free and the other rigid; the abscissa in case 
(c) is log (16T). 


P. 674, the legend to Fig. 2 should read as follows :— 

Fic. 2.—The variation of the critical Rayleigh number, R,, for the onset of instability as a 
function of Q for the three cases (i) both bounding surfaces free (curve marked 1), (ii) both bounding 
surfaces rigid (curve marked II) and (iii) one bounding surface free and the other rigid (curve 
marked X11): the abscissa in case (iii) is log (4Q). The points marked I, I and III on the R,-axis 
are the limiting values for the three cases considered in the absence of a field. 

P. 678, line 2 :— 
For the Fifteenth International Congress, read the Fifth International 
Congress. 


P. 721, equation (20) should read as follows : 


(n — 3) L G) ; a T,? ied 
Vv. ——— Li _ ind ib ay gil d 
“a | Tne = or {(: To =r, (" =)) (a | ane a) 
r ba f (a—af"a sina cd) 
T 0 or J0 
te 2 £2 Gi) sin 9 <A 
T, sind \ 26) |’ 


P. 728, second paragraph, the following addition should be made :— 


The dependence of A, on (z’—2z) near the edge of the star is in general more 
complicated than that given, but the conclusion concerning the feed-back currents holds. 








P. 728, equation (51) should read :— 


1°52 X 107? M mq 
3 R q> 
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